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§1.  Introduction.  The  Basic  Properties  of  HDA 

In  this  review  highly  dispersed  aerosols  (HDA)  are  defined  as  aero- 

.  C 

dispersed  systems  with  particulate  dimension1  <  1,000  A.  In  air  at 
atmospheric  pressure  this  corresponds  to  a  Knudsen  number  Kn  (ratio  of 
mean  free  path  of  gas  molecules  to  particle  radius)  greater  than  unity. 

The  most  important  property  of  HDA  with  Knudsen  number  Kn  1  is  that 
the  processes  of  pulse,  energy  and  mass  transfer  from  the  particles  to 
the  medium  and  back  is  described  in  this  case  by  the  gas  kinetics  formulas, 
i.e.,  the  particles  can  be  regarded  simply  as  gigantic  gas  molecules. 
Therefore,  the  resistance  of  the  medium  to  the  motion  of  the  particles, 
the  rate  of  evaporation  (mass  reduction)  and  the  rate  of  heat  transfer 
to  them,  as  well  as  the  rate  of  thermophoresis  are  proportional  to  the 
square  of  the  radius  of  the  particles.  The  "Knudsen  aerosol  number", 
i.e.,  the  ratio  of  the  apparent  mean  free  path  of  the  particles  them¬ 
selves  to  their  radius,  when  Kn  >  1 ,  is  also  greater  than  unity  and  there¬ 
fore  the  rate  of  coagulation  of  HDA  is  expressed  by  the  formula  cor  the 
number  of  collisions  between  gas  molecules,  and  the  coagulation  constant 
increases  with  particle  dimension. 

For  HDA  with  a  Knudsen  number  greatly  exceeding  unity  it  is  necessary 
to  add  to  the  above-mentioned  formulas  correction  factors  that  increase 
as  Kn  diminishes.  For  small  Kn,  however,  i.e.,  in  coarsely  dispersed 
aerosols,  all  above-mentioned  processes  are  subordinate  to  completely 
different  laws,  which  are  derived  from  the  equations  of  hydrodynamics, 
thermal  conductance  and  diffusion  in  continua. 

The  vapor  pressure  of  HDA  particles  greatly  exceeds  the  equilibrium 
pressure  above  a  flat  surface  of  the  materials  of  the  particles. 

5  In  this  article  the  dimension  is  always  particle  diameter  (d) . 
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Diffusion  of  light  by  HDA  obeys  the  Rayleigh  equation,  and  the  intensity 
of  the  diffused  light  is  so  small  that  the  analysis  of  HDA  by  macro- 
optical  methods  is  fruitless;  under  ordinary  conditions  of  illumination 
and  observation  HDA  particles  are  invisible  under  the  ultramicroscope. 

Brownian  motion  of  HDA  particles  is  so  vigorous  and  their  inertia 
and  sedimentation  rate  under  the  influence  of  gravity  so  small  that  the 
setting  of  HDA  on  walls,  obstacles,  etc.  takes  place  exclusively  as  a 
result  of  diffusion.  The  charges  taken  on  by  particles  with  radius 
<0.03  micron  in  a  bipolar  ionized  atmosphere  do  not  exceed  one  elemental 
charge,  and  the  fraction  of  the  charged  particles  decreases  with  diminish¬ 
ing  dimension. 

In  every  sense  mentioned  herein  HDA  sharply  differs  from  aerosols 
with  the  dimension  >  1,000  A,  which  up  to  now  have  been  the  basic  target 
of  aerosol  investigations. 

We  mention  also  that  with  consideration  of  the  above  statement  con¬ 
cerning  the  charges  of  HDA  particles  their  electrical  mobility  exceeds 
-4  2 

10  cm  /v*sec;  thus  the  charged  part  of  atmospheric  HDA  belongs,  in  the 
terminology  of  aerophysicists,  to  "medium'1  and  "heavy"  atmospheric  ions. 
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CHAPTER  1 


FORMATION  OF  HOA  AND  PRODUCTION  METHODS 
§2.  Formation  of  Condensation  HDA 

The  investigation  of  the  theory  of  the  formation  of  HDA  should  begin 
with  the  following  question:  In  exactly  what  cases  during  the  condensa 
tion  of  vapor  are  extremely  fine  particles  formed?  The  average  mass  of 
the  particles  is  equal  to  the  ratio  of  the  mass  of  the  condensing  vapor 
to  the  number  of  particles  formed.  It  follows  from  any  theory  of 
spontaneous  condensation  that  the  rate  of  nucleation  (i.e.,  formation 
of  seeds  of  the  new  phase)  increases  more  rapidly  than  the  first  order 
of  supersaturation,  whereas  the  rate  of  condensation  growth  of  particles 
increases  approximately  in  proportion  to  supersaturation.  Therefore, 
the  average  particle  dimension  in  the  forming  aerosol  in  the  absence  of 
coagulation  should  be  less,  the  greater  supersaturation  is,  i.e.,  the 
greater  the  extent  of  supercooling  of  the  vapor  (in  the  case  of  instant* 
aneous  temperature  drop),  iu  awl— .1  conditions  cooling  takes  place  at 
a  certain  finite  rate  and  the  formation  of  the  aerosol,  including 
nucleation  and  condensation  growth  of  the  seeds,  takes  place  during  the 
entire  cooling  period.  The  faster  the  vapor  cools  at  the  given  initial 
and  final  temperatures  the  smaller  the  amount  of  vapor  that  is  used  up 
in  the  growth  of  particles  formed  in  the  early  cooling  stage  at  a  com* 
paratively  low  supersaturation  and  the  larger  the  number  of  particles 
formed. 

As  regards  the  influence  of  the  nature  of  the  substance  on  the 
dimension  of  the  aerosol  particles,  most  important  here  is  the  slope  of 
the  curves  of  the  vapor  pressure  temperature  dependence.  Given  the  same 
cooling  conditions  of  the  vapor,  substances  with  steeper  curves  produce 
more  highly  dispersed  aerosols.  This  steepness  for  a  given  substance 
increases  as  temperature  decreases.  Therefore,  if  the  vapor  of  a  substance 
that  is  saturated  at  temperatures  T^  and  T,  (T^  >  T^)  is  instantaneously 

cooled  by  AT  degrees,  a  more  highly  supersaturated  and  more  highly  dis* 
persed  aerosol  is  formed  in  the  second  case.  Here  both  high  super¬ 
saturation  and  a  reduction  in  the  equilibrium  concentration  of  the  vapor 
in  the  forming  aerosol  arc  important.  Meanwhile,  if  this  concentration 
is  not  very  low  the  HDA  is  extremely  unstable:  within  the  HDA  the  process 
known  as  "devouring"  of  the  smaller  particles  by  the  larger  takes  place 
at  a  rapid  rate.  This  process  is  attributed  both  to  isothermic  fractiona¬ 
tion  of  the  vapor  from  fine  particles  to  large  and  to  the  effect  of 
local  temperature  fluctuations  in  the  aerosol:  even  with  an  extremely 
small  temperature  rise  the  finer  particles  evaporate  completely,  and  with 
subsequent  cooling  the  liberated  vapors  condensed  on  individual  larger 
particles.  Therefore,  the  HDA  can  be  made  from  substances  that  are 
volatile  at  room  temperature,  for  instance,  water,  by  simply  a  rapid 
cooling  of  the  vapor  to  the  temperature  of  liquid  nitrogen. 


-3- 


Ordinarily  the  formation  of  IfliA  takes  place  through  spontaneous 
condensation,  since  the  concentration  of  foreign  condensation  nuclei, 
including  gas  ions,  is  small  in  comparison  with  the  concent  ion  of 
the  nuclei  that  form  spontaneously  at  high  supersaturations.  Hiring 
the  formation  of  IU>A  in  plasma  burners,  on  the  other  hand,  during 
condensation  of  cesium  vapors  in  a  coronal  discharge  in  flow  from  a 
Lavai's  nozzle,  etc.,  i.c.,  at  extremely  high  ionic  concentrations, 
condensation  apparently  takes  place  on  the  gas  ions. 

We  excluded  from  the  above  examination  the  coagulation  of  parv.cles 
which  can  be  accompanied  by  coalescence  or  sticking  together.  Meanwhile, 
as  will  be  pointed  out  below,  coagulation  plays  a  fundamental  role  in 
the  formation  of  HDA  at  moderate  concentrations  of  the  vapor. 

Classical  nuclcation  theory  (1]  is  thoroughly  discussed  in  the 
literature.  Therefore,  we  will  simply  discuss  the  basic  premises  of 
the  theory.  It  is  assumed  tnat  a  supersaturated  vapor  represents  a 
mixture  of  single  molecules  and  molecular  aggregates,  the  concentration 
of  which,  according  to  general  theory  of  thermodynamic  fluctuations,  is 
an  exponential  function  of  their  free  enthalpy  of  formation  AG.  The 
dependence  of  AG  on  the  number  of  molecules  in  aggregate  g  in  a  super¬ 
saturated  vapor  has  the  maximum.  Each  aggregate  may,  with  a  certain 
probability,  get  larger  as  the  result  of  the  sticking  to  it  of  colliding 
molecules  cr  may  shrink  through  evaporation.  For  aggregates  with 
dimensions  exceeding  the  critical  value  g*  the  corresponding  maximum 
AG,  the  probability  of  growth  is  much  greater  than  the  probability  of 
evaporation.  Such  ’’super  critical”  aggregates  are  viewed  as  nuclei  of 
the  new  phase.  The  change  with  time  of  the  concentration  of  the  aggregate 
containing  g  molecules  is  expressed  through  the  equation 


-  >  t 
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where  dg  is  the  probability  of  collision  of  the  aggregate  with  a  vapor 

molecule,  a  is  the  probability  of  evaporation  of  molecules  from  the 

aggregate.  The  value  of  6  is  classical  theory  is  expressed  through 

K  i 

4^  -  1  using  the  principle  of  minute  equilibrium.  When  g*  is  sufficiently 

large  the  set  of  equations  of  the  term  (J.l)  for  g  from  1  to  g*  can  be 
replaced  by  the  differential  equation 
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which  is  essentially  the  equation  of  the  "diffusion”  of  aggregates  along 
discontinuous  axis  g.  Under  stationary  conditions  nuclcation  rate  1  can 


be  found  in  the  assumption  that  supercritical  particles  are  eliminated 
trout  the  system  and  are  replaced  bv  a  quantity  of  vapor  equal  in  mass 
through  the  equilibrium  equation  to 
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During  the  formation  of  I  IDA,  because  of  the  high  values  of  super¬ 
saturation  required  for  this  purpose,  g*  usually  does  not  exceed  10, 
lit  ere  fore ,  the  method  used  in  classical  theory  of  expressing  the  free 
enthalpy  of  nuclei,  in  which  they  are  viewed  as  droplets  of  a  liquid 
having  the  macroscopic  properties  of  the  latter,  is  obviously  unsuitable. 
Moreover,  when  g*  is  small  the  discreteness  of  the  g  axis  increases  and 
it  cannot  be  regarded  as  continuous.  Very  often  the  process  leads  to 
supersaturat ion ,  i.c.,  cooling  of  the  vapor  or  vapor-gas  mixture,  takes 
place  simultaneously  with  the  formation  of  nuclei  of  the  new  phase  and 
with  the  growth  of  the  particles  formed.  Under  such  conditions  the 
temperature,  concentration  of  the  vapor  and  supersaturation  change  con¬ 
stantly  throughout  the  process,  and  this  makes  it  impossible  to  solve 
equations  of  the  form  (2.2). 


There  is,  in  addition  to  the  classical  method,  another  approach  to 
the  problem  of  (indention,  known  as  the  "theory  of  constant  nucleus 
dimension"  ("constant -number  theory").  The  value  of  g*  in  this  series 
is  seen  to  be  small,  slightly  larger  than  2,  and  practically  independent 
of  nuclcation  conditions.  The  probability  of  evaporation  of  a  nucleus 
containing  g*  +  1  molecules  is  ignored  and  the  concentrations  of  the  sub- 
critical  aggregates  are  considered  equal  to  the  equilibrium  concentrations 
In  classical  theory,  however,  the  difference  between  stationary  and 
equilibrium  concent  rat  ions  of  small  aggregates  is  also  small.  The  in¬ 
tuitive  considerations  that  make  up  the  foundation  of  the  theory  of  con¬ 
stant  nucleus  dimension  were  first  advanced  by  l.aMcr  [2]  and  Courtney 
[3].  based  on  similar  considerations,  Christiansen  [4]  attempted  to 
evaluate  the  time  required  for  reaching  equilibrium  concentrations  of 
subcritical  aggregates  after  the  establishment  of  supersaturat ion  in  a 
monomolecular  vapor  by  solving  a  system  of  four  equations  of  the  form 
(2.1).  These  estimates,  however,  were  only  of  a  qualitative  character 
because  of  the  lack  of  definition  of  the  free  enthalpy  of  formation  of 
small  aggregates. 


A  method  for  calculating  this  value  was  advanced  by  Reed  [5]  who 
expressed  the  chemical  potentials  of  aggregates  of  the  molecules  of  inert 
gases  through  sums  with  respect  to  the  states  of  these  aggregates.  For 
this  purpose,  Reed  assumed  that  aggregates  of  each  dimension  exist  only 


in  the  most  stable  three-dimensional  configuration,  that  within  them 
exists  only  a  paired  relationship  between  the  nearest  neighbors, 
described  by  the  Leonnard-Johns  equation,  and  that  the  molecules  making 
up  the  aggregate  retain  their  rotational  degrees  of  freedom.  Reed 
calculated  the  chemical  potentials  of  aggregates  of  nitrogen  molecules 
with  g  <  8.  Calculation  of  larger  aggregates  is  complicated  by  a  larger 
volume  of  calculations  and  greater  error  related  to  lack  of  consideration 
of  various  three-dimensional  configurations  of  the  aggregates.  The  latter 
error  can  be  eliminated  by  using,  instead  of  Reed's  method,  methods  of 
calculating  reduced  group  integrals,  although  even  the  latter  methods 
for  such  calculations  require  a  much  larger  volume  of  calculations  than 
Reed's  method.  Sutugin  and  Fuks  [6]  advance  a  method  of  calculating 
nucleation  of  high  supersaturations,  based  on  the  combining  of  the  theory 
of  constant  nucleus  dimension  with  Reed's  calculations.  According  to 
this  method  the  process  is  broken  down  into  several  tens  or  hundreds  of 
time  intervals  At.  It  is  assumed  that  at  the  beginning  of  each  interval 
the  concentration  of  subcritical  aggregates  corresponds  to  values  that 
are  equilibrium  values  for  temperature  and  concentration  of  the  vapor 
which  exists  at  that  particular  moment.  The  number  of  nuclei  of  dimen¬ 
sion  g*  ♦  1  formed  during  time  At  is  determined  with  the  aid  of  an 
equation  of  the  form  (2.1),  but  without  consideration  of  the  terms, 
described  decomposition  of  these  supercritical  nuclei.  The  expenditure 
of  vapor  to  the  growth  of  particles  formed  in  the  preceding  intervals 
At  is  considered  simultaneously.  This  method  can  also  be  used  for  nucle¬ 
ation  of  associated  vapors.  In  this  case  it  is  also  necessary  to  consider 
the  probability  of  collisions  of  associates  with  each  other  (if  such 
collisions  lead  to  the  formation  of  supercritical  nuclei)  and  with  super¬ 
critical  particles. 

This  method  was  used  for  calculating  the  condensation  of  silver  vapors 
in  a  hot  turbulent  stream  of  argon  flowing  into  cold  stagnant  air  in  the 

range  of  silver  concentrations  from  5*10^  to  5*10^  atoms/cm ^  (weight 

3 

concentration  0.5-50  g/cm  ).  During  the  calculation  it  was  found  con¬ 
venient  to  use  the  value  of  the  integrals  equal  to  the  average  interval 
between  consecutive  collisions  between  a  critical  aggregate  and  vapor 
molecules,  i'his  simplified  calculation,  since  it  could  be  assumed  that 
during  time  At  all  aggregates  of  dimension  g*  pass  into  the  supercritical 
range  of  dimensions.  It  was  understood  that  the  value  of  At  was  not 
constant  during  the  calculation,  but  increased  as  the  vapor  was  used  up. 
Since  At  is  quite  small  ( 10~ 5  to  10~7  sec),  the  equilibrium  concentrations 
changed  not  more  than  by  two  during  transition  from  one  interval  to 
another,  and  on  the  other  hand,  the  value  of  At  was  large  enough  so  that 
the  concentration  of  aggregates  took  on  a  new  value,  corresponding  to  the 
beginning  of  a  new  interval.  Therefore,  it  could  be  assumed  that  in  each 
interval  the  concentration  of  critical  aggregates  was  equal  to  the 
equilibrium  concentration,  determined  from  the  conditions  that  existed 
at  the  beginning  of  the  interval.  The  dimension  of  the  critical  nuclei 
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under  the  above-stated  conditions  was  2  to  6  atoms,  depet  aing  on  temperature 
and  concentration,  so  that  the  term  "theory  of  constant  nucleus  dimension" 
should  be  regarded  as  incorrect  The  concentration  of  particles  formed 
was  quite  high--1013  to  1014  cm.  Therefore,  coagulation  took  place  in 
parallel  with  the  nucleation  and  condensation  growth  of  the  particles. 

The  coagulation  of  very  fine  particles  may  be  accompanied  by  their  sticking 
together  or  coalescence,  even  at  relatively  low  temperatures  and  leads  to 
the  formation  not  of  aggregates,  but  of  continuous  particles.  With  suf¬ 
ficient  smallness  of  intervals  At  the  changes  caused  in  any  such  interval 
by  any  of  three  simultaneous  processes--nucleation,  condensation  growth 
and  coagulation,  are  relatively  minor  and  it  can  be  assumed  that  these 
processes  take  place  independently,  preceding  each  other  in  time  and 
each  interval. 

In- the  literature  are  described  several  programs  [7-11]  for  calculat¬ 
ing  the  coagulation  of  polydispersed  aerosols,  based  on  the  solution  of 
a  system  of  linear  differential  equations  of  the  form 
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Present-day  computers  make  it  possible  to  solve  a  system  containing 
not  more  than  several  hundred  such  equations.  Therefore,  these  methods 
arc  suitable  only  for  the  calculation  of  processes  in  which  the  particles 
formed  contain  not  more  than  100  to  150  primary  particles.  During  the 
formation  of  HDA ,  primary  particles  are  single  molecules  and  the  smallest 
value  of  [illegible]  at  which  the  calculation  should  be  started  is  1. 

For  g  the  smallest  value  is  g*  +  for  g  <  g*  the  concentrations  are 

assumed  to  be  equilibrium  at  each  given  moment.  Since,  as  a  result  of  the 

°  C 

process,  particles  measuring  100  A  may  form,  i.e.,  consisting  of  10*>  and 
more  molecules,  the  above-mentioned  programs  are  not  suitable  for  the 
calculation  of  the  formation  of  aerosols  that  are  net  very  highly  dis¬ 
persed.  Rosinski  and  Snow  [12]  advanced  a  simplified  program  based  on 

the  assumption  that  all  particles  consist  of  21  molecules,  where  i  is  a 
whole  number.  The  collision  of  two  particles  of  the  class  21  leads  to 
the  formation  of  particles  of  the  class  21+*.  Further,  it  is  assumed 
conditionally  that  during  collisions  of  particles  of  classes  21  and  21 
particles  of  the  class  21+^  are  also  formed,  but  in  the  quantity 
(2*  +  2*  *)/21+*  =  0.75  of  the  number  of  particles  of  calss  21  that  partici¬ 
pated  in  the  collisions.  In  collisions,  however,  particles  of  classes  21 
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and  21  J,  j  >  1,  the  dimension  of  the  former  decreases,  but  their  number 

increases  by  a  factor  of  (21  +  21"^)/21.  Of  these  assumptions,  the  former 
leads  to  an  overstatement  of  the  rate  of  mass  transfer  along  the  axis  of 
particle  dimension,  and  the  latter  to  an  understatement;  the  errors 
attributed  to  these  assumptions  are  at  least  partially  compensated. 

This  program  was  combined  with  the  above-described  method  of  cal¬ 
culating  nucleation  and  condensation  growth  [6].  Both  the  calculation 
and  experimental  results  obtained  by  electron  microscopic  analysis  of 
particles  sampled  at  various  stages  of  the  formation  of  aerosols.  The 
determination  of  the  capacity  of  the  powders  formed  by  the  precipitation 
of  these  aerosols  for  coalescence  prompted  the  conclusion  that  the  process 
by  which  the  HDA  of  silver  is  formed  in  this  range  of  concentrations  of 
the  vapor  can  be  broken  down  into  six  stages.  1.  Nucleation  and  con¬ 
densation  growth.  2.  Nucleation,  condensation  growth  and  coagulation, 
accompanied  by  the  aglutenation  of  particles.  3.  Coagulation  with 
aglutenation  and  condensation  growth.  4.  Coagulation  with  aglutenation. 

5.  Coagulation,  accompanied  by  aglutenation  of  fine  particles  and  co¬ 
alescence  of  large.  6.  Coagulation,  leading  to  the  formation  of  ag¬ 
gregates  of  particles  held  only  by  cohesion  forces.  The  dispersion  of 
the  aerosols  obtained  in  this  work,  determined  by  measurement  of  the 
specific  area  of  the  sediments,  was  governed  almost  exclusively  by  co¬ 
agulation,  accompanied  by  aglutenation  of  the  particles.  Therefore, 
the  average  particle  dimension  was  propcrtional  to  the  cube  root  of  the 
vapor  concentration. 

At  very  large  supersaturations  of  the  vapor  the  labile  state  may  be 
reached,  i.e.,  the  thermodynamic  barrier  between  the  phases  may  vanish, 
and  in  this  case  the  formation  of  aggregate,  even  from  two  molecules, 
becomes  an  irreversible  process.  The  assumption  of  the  reaching  of  the 
labile  state  was  made  arbitrarily  d>-  Rozinski  and  Snow  [12]  for  calculation 
of  the  condensation  of  Fe20s  vapors  and  by  Stockham  [13]  for  silver  vapors. 

By  calculating  the  condensation  of  Ag  vapors  under  conditions  similar  to 
Stockham's  experimental  conditions,  Sutugin  and  Fuks  found  that  the  labile 
state  in  this  case  is  not  achieved.  However,  since  in  the  final  analysis 
the  dispersion  is  determined  by  coagulation,  accompanied  by  particle 
aglutenation,  in  view  of  the  asymptotic  character  of  coagulation  the  path 
of  condensation  in  the  early  stages  of  the  formation  of  the  aerosol  has 
little  effect  on  the  final  result.  This  governs  the  applicability  of  the 
Rozinski-Snow  method  to  those  cases  where  the  labile  state  is  not  reached, 
and  explains  the  fact  that  Stockham  obtained  reasonable  results. 

In  the  case  where  g*  =  2  the  rate  of  nucleation  is  determined  by  the 
ratio  between  the  rates  of  formation  of  trimer  and  dissociation  of  the 
dimer.  The  formation  of  the  bonded  state  during  collision  of  two  molecules 
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was  examined  .[16] ,  although  the  lifetime  of  the  vapor  aggregate  formed 
as  a  function  of  the  temperature,  reaction  potential  and  geometric 
parameters  of  collision,  could  be  found  only  by  numerical  methods.  There¬ 
fore,  the  results  [16]  for  the  calculations  of  the  process  of  aerosol 
formation  are  still  inapplicable.  It  should  be  pointed  out  that  the 
concentration  of  dimers  in  nitrogen,  calculated  by  this  method,  differ 
somewhat  from  the  concentration  calculated  by  Reed's  method  [6]. 

In  those  cases  where  g*  is  quite  large  and  supersaturation  changes 
little  during  the  relaxation  time  of  the  nucleation  process,  as  determined 
by  the  solutions  of  time-dependent  diffusion  equation  (2.2),  it  is  pos¬ 
sible  to  use  the  ordinary  stage  by  stage  calculation  method,  in  which 
the  nucleation  rate  in  each  stage  is  determined  by  stationary  solutions 
of  classical  theory.  This  method,  of  course,  is  not  applicable  in  the 
presence  of  coagulation,  nor  for  description  of  the  condensation  of 
associated  vapors.  Using  this  method  Griffin  and  Sherman  [15]  demon¬ 
strated  that  HDA  can  form  as  a  result  of  condensation  of  copper  vapors 
that  flow  through  a  Laval's  nozzle  at  high  mach  numbers.  Such  a  descrip¬ 
tion  of  consecutive  calculations  of  condensation  in  the  nozzles  using 
stationary  solutions  of  classical  theory  can  be  found  in  the  literature 
[16,  17]. 

Under  the  condition  of  slow  change  in  supersaturation  the  condensation 
process  can  be  described  with  the  aid  of  differential  equations  that 
express  changes  in  time  of  temperature,  vapor  concentration  and  nuclea¬ 
tion  rate,  and  also  of  the  equations  of  continuity  and  adiabatic  expansion 
in  the  case  of  condensation  in  nozzles.  The  solution  of  such  an  equation 
system  for  the  simplest  cases  and  the  large  number  of  assumptions  can  be 
found  analytically,  i.e.,  the  theoretical  dependence  of  the  number  of 
particles  formed  and  their  mass  on  time  can  be  found.  Solutions  of  this 
type  were  found  by  several  authors  [18-20].  The  attempt  of  Buykov  and 
Bakhanov  [21]  to  avoid  simplification  assumptions  required  the  use  of 
numerical  solution  methods  with  the  aid  of  computers. 

Very  often  the  formation  of  condensation  HDA  is  related  to  chemical 
reactions  in  the  gas  phase,  for  instance  during  the  production  of  soot 
by  the  pyrolysis  of  hydrocarbons,  production  of  metals  by  decomposition 
of  chlorides  and  carbonyls,  etc.  Three  cases  of  the  formation  of  con¬ 
densation  aerosols  as  a  result  of  chemical  processes  can  be  singled  out. 

The  first  case  is  where  the  chemical  reaction  leads  to  the  formation  of 
condensed  vapor  and  takes  place  so  rapidly  that  it  is  finished  before 
nucleation  can  begin.  In  this  case  the  calculation  of  the  condensation 
process  is  the  same  as  that  of  "physical"  condensation.  In  the  second 
case  the  reaction  takes  place  in  parallel  with  condensation,  but  the 
aerosol  particles  formed  have  no  effect  on  the  reaction.  In  this  case 
the  chemical  reaction  can  be  taken  into  account  by  the  consecutive 
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calculation  of  condensation.  For  this  purpose  it  is  necessary  to  consider 
not  only  the  expenditure  of  vapor  for  growth  and  formation  of  the  particles, 
but  also  the  formation  of  new  portions  of  vapor  as  a  result  of  the  chemical 
reaction,  and  also,  if  necessary,  the  thermal  effect  of  the  reaction.  A 
calculation  program  of  this  type  was  described  by  Dunhem  [22],  who  was  able 
in  this  way  to  establish  the  relationship  between  the  rate  of  formation  of 
H2S0^  and  the  rate  of  nucleation  in  these  vapors.  The  third  case  is  where 

the  reaction  can  be  catalyzed  on  the  surface  of  the  forming  particles. 

This  case  poses  the  greatest  difficulties  in  calculation  and  so  far  no 
mathematical  model  has  been  devised  for  such  a  process.  Nevertheless,  it 
should  be  pointed  out  that  Tesner  [23]  proposed  a  method  permitting  cal¬ 
culation  of  the  dependence  of  the  number  of  soot  particles  formed  by 
pyrolysis  on  time  if  the  distribution  of  the  particles  with  respect  to 
dimensions  in  the  final  aerosol  is  known. 

§3.  Production  of  HDA  by  Physical  Condensation  of  Vapor 

The  methods  of  producing  condensation  HDA  can  be  classified  on  the 
basis  of  the  method  of  achieving  supersaturation  of  the  vapor:  during 
"physical  condensation"  supersaturation  is  achieved  by  coo.ing  the  vapor 
or  gas-vapor  mixture,  during  "chemical"  condensation,  a  substance  possess¬ 
ing  extremely  low  vapor  tension  at  the  reaction  temperature  is  formed  as 
the  result  of  the  gas  reaction.  The  vapor  can  be  cooled  by  mixing  it  with 
a  cold  gas,  by  heat  exchange  with  a  cold  body  (walls  of  refrigerator, 
coolant,  etc.),  by  adiabatic  expansion,  or  by  a  combination  of  these 
methods . 

The  simplest  method  of  HDA  generation  by  mixing,  but  hardest  to  control, 
is  heating  of  high-boiling  solids  or  liquids  in  a  turbulent  flow  of  cold 
gas.  Mixing  of  the  vapor  of  the  heated  material  and  the  gas  and  the 
initial  nucleation  stage  take  place  here  near  the  surface  of  the  material. 

In  the  absence  of  flow  or  in  laminar  flow  of  the  gas  mixing  takes  place 
very  slowly  and,  in  accordance  with  what  was  said  in  §2,  coarsely  dispersed 
aerosols  are  generally  formed.  In  this  method  the  presence  of  more 
volatile  contaminants  which  also  form  HDA  on  heating  pose  serious  dif¬ 
ficulties,  although  in  most  cases  these  aerosols  can  be  distinguished 
from  the  basic  aerosol  on  the  basis  of  the  criterion  that  their  concentra¬ 
tion  gradually  diminishes  by  measure  of  evaporation  of  the  contaminants 
[24],  When  the  concentration  of  the  aerosol  remains  constant  over  a  long 
period  of  time  it  can  be  assumed  that  it  consists  of  the  basic  material. 

The  method  is  often  formulated  in  such  a  way  that  the  vaporized  substance 
is  applied  as  a  thin  film  on  a  high-melting  wire  heated  by  an  electric 
current.  This  is  the  way  that  HDA  is  obtained  from  various  salts.  A 
serious  deficiency  of  this  method  is  the  fact  that  in  time  the  concentra¬ 
tion  of  the  aerosol  formed  and  the  particle  dimension  continually  decrease. 
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In  the  case  of  NaCl  at  a  temperature  of  450°C  the  particle  dimension 
diminishes  for  5  to  7  houis,  reaching  values  of  100  to  300  K  [25].  This 
is  attributed  to  the  gradual  recrystalization  of  the  salt  on  heating, 
leading  to  a  reduction  in  its  rate  of  evaporation.  To  obtain  time 
constant  aerosols  it  is  possible  to  resort  to  the  aging  of  the  NaCl 
la;  _-r  for  several  hours  and  replacing  it  after  one  day  of  work.  The  HDA 
thus  obtained  is  very  polydispersed. 

HDA  can  also  be  generated  by  the  calcination  in  air  of  pure  metallic 
wires;  first  purpose  Pt  and  W  are  most  often  used.  The  concentration  of 
the  HDA  formed  is  quite  constant  in  time  and  is  more  or  less  satisfactorily 
reproducible,  thanks  to  which  these  types  of  HDA  are  often  used  as  con¬ 
densation  seeds  and  for  the  production  of  monodispersed  aerosols.  The 
mechanism  of  the  formation  of  HDA  from  the  above-mentioned  metals  consist, 
apparently,  in  their  surface  oxidation  and  evaporation  of  oxides  that  are 
considerably  more  volatile  than  the  metals  themselves.  This  is  evidenced 
by  the  fact  that  the  formation  of  oxide  HDA  begins  at  temperatures  at 
which  the  vapor  pressure  of  the  metals  is  infinitesimally  small:  from 
platinum  at  200  [28]  or  255°C  [271;  from  tungsten  at  300  or  500°C  [281; 

from  nickel  at  650°C  [27]  and  950°C  [24].  Such  great  differences  in  the 

determination  of  the  temperature  at  which  HDA  begins  to  form  are  attributed 
to  the  difficulty  in  detecting  aerosol  particles  of  nearly  atomic  size  and 
to  the  presence  of  hard  to  remove  contaminants  in  the  metals.  These, 
aersols  are  also  extremely  polydispersed  [27-29].  During  the  blowing 
of  a  platinum  wire  heated  to  300°C  with  a  stream  of  air  at  the  rate  of 
1  liter  per  minute  the  concentration  of  aerosol  reached  ~  2*105  cm3  and 

O 

the  average  particle  dimension  was  about  40  A  [28] .  A  CrO^  aerosol  with 

a  particle  dimension  of  ~15  A  [29]  was  obtained  from  a  Nichrome  wire 

heated  to  1,000°C. 

The  analogous  method  is  sublimation  of  a  substance  in  a  rarefied 
gas:  the  substance,  placed  in  a  vessel  with  a  rarefied  inert  gas  is 
calcined  either  with  a  tungsten  coil  [30-32]  or  by  the  inducation  method 
[33,  34].  This  method  was  first  posed  by  Gen,  Ziskin  and  Petrov  [30] 
for  the  production  of  aerosols  of  aluminum,  and  then  was  applied  to  a 
large  number  of  other  metals.  The  dimension  of  the  particles  generated 
by  this  method  diminishes  with  a  reduction  in  the  rate  of  evaporation  of 
the  metal  and  pressure  of  the  gas,  where  the  minimum  average  dimensions 
of  the  particles  obtained  by  this  method  is  *”200  A.  More  highly  dispersed 
aerosols  could  not  be  obtained,  since  at  low  vapor  concentrations  the 
dimension  of  the  particles  is  practically  independent  of  the  gas  pressure 
and  concentration  of  the  vapor.  Turkevich  [34]  analyzed  at  length  the 
mechanism  of  the  formation  of  aerosols  by  this  method.  By  placing  the 
apparatus  for  the  sampling  of  aerosol  particles  at  various  distances 
from  the  vapor  source  he  was  able  with  the  help  of  electron  microscopy 
to  trace  the  process  of  the  formation  of  aerosol  or,  at  least,  the  last 
stages  of  its  formation  and  arrive  at  the  conclusion  that  the  aerosol 
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is  formed  in  three  stages,  corresponding  to  the  first,  third  and  fifth 
stages  on  p.  12  of  this  review.  A  reduction  in  the  dimension  of  the 
particles  as  the  pressure  of  the  inert  gas  is  increased,  is  in  Turkevich's 
opinion,  the  result  of  accleration  of  vapor  diffusion  and  the  resultant 
increase  in  the  rate  of  its  cooling. 


In  [13]  is  described  a  flow  version  of  this  method:  silver  vapor  from 
a  crucible  heated  by  the  induction  method  is  mixed  with  rarefied  argon 
flowing  past  the  crucible.  In  this  way  it  is  possible  to  produce  silver 
aerosols  with  an  average  particle  dimension  >40  A.  The  particle  dimension 
increases  as  the  vapor  concentration  is  increased,  but  is  practically 
independent  of  the  argon  pressure.  Probably,  the  mixing  of  vapor  in 
this  case  with  the  cold  gas  takes  place  basically  as  the  result  of 
turbulence . 


The  production  of  aerosols  in  rarefied  gases  is  in  itself,  of  course, 
extremely  useful  in  the  analysis  of  the  properties  of  aerosols  at  reduced 
pressure,  since  the  problem  of  a  substantial  change  in  pressure  in  an 
aerosol  without  the  loss  of  particles,  their  coagulation,  etc.  is  extremely 
complex.  For  this  very  reason,  however,  this  method  is  not  suitable  for 
investigation  of  aerosols  under  normal  conditions. 

More  improved  HDA  generators,  in  which  comparatively  large  quantities 
of  the  substance  are  heated  ir.  a  thermostat  furnace  in  a  slow  gas  flow, 
are  now  in  use.  Here  the  gas  is  practically  saturated  with  the  vapor 
and  is  then  rapidly  mixed  with  cold  gas.  Under  these  conditions  the  con¬ 
centration  of  the  vapor  remains  constant  in  time,  which  facilitates 
constancy  of  concentration  and  dispersion  of  the  aerosol  formed.  Two 
versions  of  such  generators  are  now  in  use.  In  one  version,  mixing  is 
accomplished  in  a  free  gas  volume,  for  instance,  with  the  outflowing  of 
a  hot  turbulent  stream  of  vapor  or  vapor-gas  mixture  into  a  large  volume 
of  cold  gas.  In  this  case  all  the  vapor  is  condensed  in  the  volume.  In 
the  second  version  mixing  takes  place  in  a  T-type  diffuser  nozzle,  etc., 
and  most  of  the  vapor  by  far  is  condensed  on  the  mixer  walls. 

A  generator  of  the  first  of  these  two  types  was  used  in  the  above- 
-mentioned  work  [6].  The  dependence  of  the  dispersion  of  the  HDA  on  the 
operating  mode  of  such  a  generator  was  investigated  by  way  of  example  of 
an  aerosol  obtained  through  the  condensation  of  MoO^  vapor  in  a  hot  air 

stream  flowing  into  cold  air  at  the  rate  of  8  to  580  m/sec  [35].  In  the 

range  of  vapor  concentrations  in  the  stream  from  0.1  to  100  g/cm^  it  was 
possible  to  obtain  aerosols  with  an  average  particle  dimension  70  to 

1,000  A.  The  particle  dimension  decreases  to  some  limit  as  the  vapor 
concentration  is  decreased  and  the  linear  flow  rate  is  increased,  and  at 
some  constant  linear  flow  rate,  as  the  nozzle  diameter  is  reduced.  It 
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is  interesting  to  note  that  the  dispersion  of  aerosols  formed  during 
condensation  in  a  stream  at  the  very  same  vapor  concentrations  was  higher 
than  in  a  heat  exchange  generator  (see  below).  The  influence  of  flow 
rate  and  nozzle  diameter  was  explained  with  the  aid  of  the  equation 
presented  in  work  [6]  for  mixing  rate  in  a  submerged  stream,  showing 
that  this  rate  increases  with  the  indicated  change  in  the  above-mentioned 
parameters.  In  the  range  of  low  concentrations  (<•  1  g/ m3)  the  increase  in 
particle  dimension  as  the  way  concentration  is  increased  is  not  very 
noticeable,  but  intensifies  as  concentration  is  increased.  This  can  be 
explained  by  the  increasing  influence  of  coagulation,  accompanied  by  co¬ 
alescence  of  the  particles. 

The  formation  of  HDA  in  arc  or  plasma  burners  takes  place  in  much 
the  same  w-iy.  The  substance  is  added  either  to  the  material  of  the 
electrodes  or  in  the  form  of  a  powder  to  the  gas  stream.  The  stream  of 
incandescent  gases  is  cooled  as  a  result  of  mixing  with  cold  gas  or  with 
atomized  coolant.  Condensation  of  the  vapors  on  cooling  takes  place 
apparently  on  ions.  According  to  the  given  electron  microscopic  analyses 
[34,  36-37],  particle  dimension  in  the  aerosols  thus  obtained  is  25  to 
0 

1,000  A,  and  the  average  dimension,  b<  terms  of  mass,  may  vary  from  200  to 

O 

800  A.  By  vaporizing  metals  in  an  atmosphere  of  inert  gases  and  in  a 
reducing  atmosphere  it  is  possible  to  produce  metal  aerosols,  and  in  an 
oxidizing  atmosphere,  oxide  aerosols.  It  is  interesting  to  note  that 
particles  of  aerosols  generated  with  the  aid  of  arc  or  plasma  burners  are 
usually  amorphous.  Holmgren,  Gibson  and  Sheer  [37]  showed  that  the  dis¬ 
persion  of  aerosols  can  be  increased  by  using  induced  cooling  of  the  gases 
flowing  from  the  arc.  Selover  [38]  investigated  the  formation  of  nickel 
aerosol  during  cooling  of  argon  plasma  containing  nickel  vapor.  To 
argon,  fed  into  a  plasma  burner,  atomized  nickel  carbonyl  was  added. 

Since  the  decomposition  of  carbonyl  took  place  at  a  very  high  temperature 
it  can  be  assumed  that  it  preceded  physical  condensation.  In  the  range 
of  carbonyl  concentrations  between  0.006  and  0.025  g/  the  particle 
dimension  of  the  aerosol  obtained  varied  within  the  limits  of  15  to  150  A. 
For  the  purpose  of  rapid  cooling  the  argon  stream  was  sprayed  with  a 
coarsely  dispersed  aerosol  of  cooled  isopropyl  alcohol.  With  an  average 
particle  diameter  (determined  by  BET  (Brunauer,  Emmett  and  Teller)  method), 
of  200  X  practically  all  particles  were  included  in  the  dimensional  range 
of  30  to  300  A.  Safronov  [39]  showed  that  it  is  possible  to  produce  a 
highly  dispersed  soot  aerosol  by  condensing  carbon  vapors  formed  by  the 
decomposition  of  benzene  in  plasma. 

By  the  plasma  method  it  is  possible  to  put  even  high-blowing  sub¬ 
stances  (with  the  exception  of  thermally  unstable  materials)  into  the 
aerosol  state,  and  in  large  quantities  at  that.  Unfortunately,  the 
mechanism  of  processes  that  take  place  during  condensation  in  cooled 
plasma  jets  has  not  yet  been  sufficiently  investigated. 
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Let  us  turn  now  to  generators  in  which  vapor  and  gas  are  mixed  in 
narrow  mixers.  A  mixer  generator  is  described  [40]  that  is  designed  for 
the  production  of  highly  dispersed  aerosols  from  salts,  for  instance, 
from  NaCl,  consisting  of  a  quartz  tube  placed  in  an  electric  furnace, 
with  one  side  fitted  with  a  T-pipe  with  an  inside  diameter  of  8  ram.  In 
the  tube  are  placed  ceramic  rings,  coated  with  a  volatile  substance,  and 
through  it  is  passed  air  at  the  rate  of  0.1  to  0.5  Z/min.  In  the  T-pipe 
this  flow  is  mixed  with  a  large  volume  of  cold  air  (5  to  20  Z/min).  To 
prevent  condensation  of  the  vapors  on  the  walls  prior  to  mixing  a  heater 
is  installed  in  the  generator.  This  generator  is  capable  of  producing 
rather  monodispersed  NaCl  aerosols  with  an  average  particle  dimension  of 
14  to  100  X.  The  maximum  calculated  concentration  of  the  particles  which 
can  be  achieved  increases  with  dimension  and  varies  between  103--109/cm3. 

By  changing  flow  parameters--f low  rate  of  air  through  the  furnace  and 
through  the  T-pipe,  current  strength  in  the  furnace  windings  and  heater, 
it  is  possible  to  change  the  average  dimension  and  concentration  of  the 
aerosol  particles  in  the  above-mentioned  limits.  '  The  generator  is  very 
stable  in  operation  and  produces  satisfactorily  reproducible  results. 

A  generator  of  this  type  was  used  successfully  [41]  for  the  production 
of  radioactive  NaCl  HDA. 

The  production  of  reproducible  monodispersed  aerosols  of  dioctyl 
sebacinate  with  3n  average  particle  radius  of  30  to  300  X  is  described 
[42].  These  aerosols  were  obtained  in  a  mixer  generator  of  the  type 
KUST,  described  in  §6,  and  NaCl  with  an  average  particle  dimension  of 
34  X  was  used  as  the  seeds.  The  particle  dimension  of  the  dioctyl 
sebacinate  could  easily  be  changed  by  varying  its  vaporization  temperature, 
concentration  of  condensation  nuclei  and  ratio  of  the  flow  rates  of  the 
gas-vapor  mixture  and  cold  air. 

Cooling  of  the  vapor-gas  mixture  was  possible  also  by  heat  exchange 
with  the  walls  of  the  apparatus.  Precisely  the  same  method  was  used  in 
the  well-known  Sinkler-LaMer  monodispersed  aerosol  generator.  LaMer  [43] 
showed  that  it  is  possible  to  use  this  generator  for  the  production  of 
sulphuric  acid  HDA.  According  to  his  report  it  is  thus  possible  to 
generate  the  HTA  with  a  particle  diameter  of  20  X. 

Matiyevich  and  Kerker  et  al.,  [44-46]  modified  the  LaMer  generator 
and  used  it  for  the  production  of  aerosols  of  various  mineral  salts,  such 
as  AgCl,  NaF,  NaCl;  these  authors  analyzed  in  detail  the  influence  of  the 
gas  flow  rate  in  the  generator  and  the  evaporator  temperature  of  the 
dispersion  of  the  aerosol  formed  under  conditions  of  spontaneous  con¬ 
densation  and  discovered  that  the  average  particle  dimension  increases 
as  these  parameters  are  increased.  The  reason  for  the  increase  of  the 
particle  dimension  with  increasing  flow  rate  is  not  completely  clear; 
unfortunately  no  detailed  plan  of  the  modified  generator  has  yet  been 
published.  In  these  works  the  authors  managed  to  produce  particles  measur¬ 
ing  between  lu  and  1,000  X.  The  degree  of  polydispersion  of  the  aerosols 
depended  little  on  the  average  particle  dimension  and  the  standard  geo¬ 
metric  deviation  was  1.15  to  1.2. 
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In  a  number  of  works  heat  exchanger  generators  were  used  for  the 
production  of  NaCl  HDA  with  high  weight  concentrations  (up  to  1  g/cm3) 
for  the  purpose  of  further  precipitation  of  these  aerosols  and  production 
of  highly  dispersed  powders  [47-49].  An  increase  in  the  particle  dimen¬ 
sion  with  increasing  vapor  concentration  was  also  noticed  in  these  works. 

Due  to  heat  exchange  it  was  impossible  to  achieve  very  high  cooling  rates, 
and  consequently  high  vapor  supersaturations.  Therefore  it  is  difficult 
to  produce  aerosols  with  particle  dimensions  less  than  100  X  in  heat 
exchanger  generators.  The  formation  of  MoO^  aerosols  by  passing  hot  air 

containing  vapors  of  this  substance  at  a  rate  of  1  to  28  1/ min  through 
a  quartz  water  cooler  with  a  diameter  of  20  min  was  investigated  in  [35]. 

At  greater  concentrations  of  0.002--8  g/m3  the  average  dimension  of  the 
particles  obtained,  determined  from  a  certain  BET  method  for  specific 
area  of  powders  obtained  by  the  precipitation  of  these  aerosols,  was 
130—1,000  A.  The  investigators  were  not  able  to  observe  whether  or  not 
the  flow  rate  has  any  clear  influence  on  the  dispersion  of  the  aerosols. 

In  addition  to  heat  exchange  and  mixing  with  cold  gas,  adiabatic 
expansion  may  also  be  used  as  a  method  of  cooling.  The  condensation  of 
the  vapor  during  adiabatic  cooling  takes  place  in  the  nozzles  of  super¬ 
sonic  wind  tunnels,  where  it  is  an  undesirable  phenomenon,  the  elimina¬ 
tion  of  which  requires  implementation  of  special  measures  (air  drying, 
heating  of  nozzle  walls,  etc.).  On  the  other  hand,  the  formation  of  HDA, 
particularly  of  strongly  contaminated  HDA,  in  nozzles  is  of  considerable 
interest  for  the  creation  of  electric  colloidal  motors. 

For  all  substances  except  helium,  a  drop  in  the  equilibrium  pressure 
of  the  vapor,  attributed  to  a  reduction  in  the  temperature  during  adiabatic 
expansion,  takes  place  more  rapidly  than  a  drop  in  static  pressure.  There¬ 
fore  sufficiently  strong  isoenthropic  expansion  of  the  vapor-gas  mixture 
or  even  of  the  pure  gas,  may  lead  to  the  creation  of  high  supersaturations. 
The  growth  rate  of  supersaturation  can  be  extremely  high,  therefore  the 
formation  of  HDA  even  from  such  volatile  material  as  water  is  observed  in 
wind  tunnels.  The  lifetime  of  these  aerosols,  however,  is  quite  short, 
since  articles  evaporate  at  high  temperatures  in  the  shock  wave.  The 
dimension  of  ice  crystals  formed  in  a  laval  nozzle  at  mach  5  was  estimated 
[50-51]  on  the  basis  of  their  inertial  properties  and  rate  of  evaporation 
in  the  shock  wave.  Particle  dimension  was  aboi't  14  A,  and  they  were  about 
identical  in  magnitude,  independent  of  the  concentration  of  water  vapor. 
Condensation  in  Laval  nozzles  was  examined  [57,  53]  from  the  point  of 
view  of  the  possibility  of  producing  a  ray  of  charged  particles  suitable 
for  use  in  an  elect rorcactivc  motor.  Hg,Cl , ,  HgCl,  and  A1C1 _  vapors  were 

passed  through  a  heated  nozzle  in  a  vacuum  [57].  For  analysis  the  particles 
were  precipitated  on  a  screen  placed  in  the  nozzle.  In  this  manner  it 
was  possible  to  obtain  particles  with  an  average  dimension  of  100  to 
1,000  X.  Cadmium  vapors  flowing  through  a  nozzle  were  condensed  in  the 
field  of  a  positive  coronal  discharge  (55),  whereupon  particles  measuring 
-lOO  X  were  obtained. 
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Cooling  of  the  vapor-gas  mixture  due  to  both  adiabatic  expansion 
and  mixing  with  cold  gas  and  ray  emission  apparently  takes  place  during 
the  production  of  aerosols  by  the  flashing  wires  method  [54-56].  In 
this  method  a  high  voltage  discharge  of  considerable  power  (about  1  kg 
per  milligram  of  wire)  is  fed  to  a  thin  wire  fastened  between  massive 
electrodes.  The  wire  is  converted  instantly  into  a  plasma  "cable",  on 
cooling  of  which  an  aerosol  is  formed.  In  an  atmosphere  of  inert  gases 
of  all  metals,  and  in  air,  metal  aerosols  are  formed  only  from  the 
noble  metals.  In  the  presence  of  nitrogen  and  oxygen  nitrides  and 
oxides  arc0formed.  The  size  of  the  particles  usually  falls  between  SO 
and  1,000  A.  The  average  particle  dimension  depends  upon  the  material 
from  which  the  wire  is  made  and  diminishes  as  the  power  of  the  discharge 
is  increased.  The  fact  that  during  instantaneous  dilution  of  the  aerosol 
cloud  it  is  possible  to  produce  aerosol  consisting  of  individual  non- 
aggregated  spheres  [54]  indicates  that  condensation  growth,  and  not 
coagulation,  plays  the  basic  role  in  the  formation  of  the  aerosol. 

Since  the  calculated  concentration  in  the  cloud  formed  by  the  flash  is 
still  great  (about  10n/cm3),  coagulation  still  takes  place  in  the 
aerosol  that  is  not  yet  cooled,  and  therefore  it  is  not  accompanied  by 
the  coalescence  of  particles.  An  aerosol  generator  is  described  in  [41] 
in  which  cooling  of  the  vapor-gas  flow  also  takes  place  as  a  result  of 
two  processes- -adiabatic  expansion  and  mixing  with  a  cold  gas.  In  this 
instrument  a  hot  gas-vapor  stream  is  passed  at  a  high  rate  through  a 
nozzle  into  a  cold  gas.  It  was  possible  to  obtain  by  this  method  rather 
monodispersed  aerosols  of  selenium,  without  nuclei,  with  particle  dimen¬ 
sions  of  600  to  2,600  X. 

§4.  HDA  Production  by  Chemical  Reactions  in  Gas  Phase 

Such  processes  as  thermal  and  photochemical  decomposition  of  carbonyl 
vapors,  thermal  hydrolysis  of  the  vapors  of  metal  chlorides  and  silicon 
tetrachloride,  thermal  destruction  of  hydrocarbons,  have  long  been  used 
on  industrial  scales  for  the  production,  through  the  aerosol  state,  of 
highly  dispersed  metal  powders  and  their  oxides,  carbon  black  and  AeTosil. 
Pertaining  to  all  these  processes  !s  a  vast  amount  of  patent  literature, 
examination  of  which  extends  beyond  the  scope  of  this  review.  Therefore 
we  will  discuss  only  a  few  works,  the  authors  of  which  focus  their  atten¬ 
tion  on  the  mechanism  of  formation  of  aerosols  during  these  processes. 

The  dispersion  of  aerosols,  metals  and  their  oxides  bbtained  in  the  works 
discussed  below  was  characterized  by  a  specific  sediment  area  measured 
by  the  BET  method. 

Photolysis  of  the  vapors  of  iron  pentacarbonyl  in  the  presence  of 
oxygen  has  been  used  repeatedly  for  the  production  of  aerosols  [57-62]. 

It  has  been  learned  [61]  that  one  quant  of  light  causes  the  separation 
of  one  CC  molecule  from  the  carbonyl  complex,  after  which  the  Fe^  (CQ)g 

complex  is  formed,  which  is  oxidized  in  the  presence  of  0^  to  Fe^Oj. 
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Kogan  [63]  used  these  reactions  for  the  production  of  "molecular  condensa¬ 
tion  nuclei".  The  process  of  the  format  ion  of  Fe,,03  was  anal/zed  in  depth 

[64].  It  was  established  that  as  the  ti  .  of  illumination  of  a  mixture 
of  nitrogen  and  the  vapors  of  iron  pentacarbonyl  by  a  mercury  lamp  is 
increased  by  0.1  to  6  seconds  the  degree  of  decomposition  of  the  latter 
increases  from  10  to  100%.  After  illumination  the  nonreactive  carbonyl 
was  absorbed  by  passing  the  aerosol  through  a  layer  of  activated  char¬ 
coal,  which  reduced  the  effect  of  autocatalytic  decomposition  of  the 
carbonyl  on  the  surface  of  the  particles.  At  a  carbonyl  concentration 
of  36.2  g/cm1  and  degree  of  decomposition  10  to  100%  it  w-s  possible  to 
produce  in  this  manner  HDA  with  an  average  particle  dimension  of  22  to 

82  X. 


The  production  of  the  aerosols  of  metal  oxides  during  thermal  hydroly¬ 
sis  of  the  vapors  of  their  chlorides  was  investigated  by  a  group  of 
French  researchers  [6S,  66]  who  fed  a  hydrogen  burner  with  oxygen  sat¬ 
urated  with  these  vapors.  The  concentration  of  the  chlorides  varied 
from  3  to  12  mg H,  the  temperature  of  the  oxygen -hydrogen  flame  was  1,600 
to  3,000"K.  The  average  particle  dimension  of  the  precipitates  of  the 
HDA  varied  from  100  to  1,200  X  and  decreased  as  the  vapor  concentration 
was  reduced.  More  highly  dispersed  aerosols  were  obtained  in  the  flame 
with  the  highest  temperature  of  1,500°K  than  in  the  flame  measuring 
3,000°K,  but  these  aerosols  were  coarser  than  at  1,900“K.  At  high  temp¬ 
eratures  the  dispersion  is  reduced,  probably  because  of  the  aglutenation 
of  particles  during  coagulation;  at  exceedingly  low  temperatures  hydrolysis 
takes  place  slowly  and,  accordingly,  the  rate  of  nucleation  is  lower.  The 
coarsest  particles  were  formed  in  the  largest  flame,  and  this  is  under¬ 
standable  considering  that  the  rate  of  cooling  of  a  submerged  jet  is 
reduced  when  its  diameter  is  increased. 

Zakutinskiy  and  Blyckher  [67]  analyzed  thermal  hydrolysis  of  AlClj 

vapor  during  mixing  of  a  string  of  nitrogen  containing  this  vapor  with 
air  saturated  with  water  vapor.  Hydrolysis  took  place  at  temperatures 
between  400  and  4S0°C.  Here,  more  highly  dispersed  aerosols  were  formed 
than  during  hydrolysis  in  a  flamc--at  chloride  concentrations  of  5  to 
30  mg fl  the  average  particle  dimension  was  SO  to  4Q0  X.  Obviously  the 
low  temperatures  at  which  Zakutinskiy  and  Blyakhcr  worked  precluded  the 
aglutenation  of  the  particles  and  resulted  in  greater  dispersion  than  ir 
works  (6S,  66] . 

A  review  of  works  pert  aining  to  the  product  ion  of  carbon  black  and 
Aerosii  can  be  found  in  [t»8,  t9],  and  the  formation  of  M>A  by  radiolysis 
of  gaseous  hydrocarbons  is  discussed  in  [70]. 

Fuks  and  Oshman  [71]  described  production  of  H,SO,  HDA  by  mixing, 
ir.  a  capillary  T-pipe,  flows  of  air  containing  SO*  and  H,0.  Investigation 
of  the  dispersion  of  these  aerosols  is  described  in  £?. 


-17- 


Si5.  Product  ion  of  Dispersed  HDA 


The  atomization  of  liquids  produces  particles  of  various  dimensions 
including  extremely  fine  particles.  However,  a  real  HDA  that  does  not 
contain  coarse  particles  at  all  or  contains  only  a  few  coarse  particles 
cannot  be  obtained  by  atomization.  Likewise,  HDA  cannot  be  obtained  by 
the  atomization  of  polydispersed  powders,  since  due  to  strong  adhesion 
among  the  fine  particles  total  deaggregation  of  such  powders  is  impos¬ 
sible.  By  atomizing  extremely  dilute  solutions  and  drying  the  droplets 
of  the  latter,  however,  it  is  possible  to  produce  HDA  from  a  dissolved 
substance.  This  method  was  used  [72]  for  the  generation  of  aerosols  of 
different  salts  with  an  average  particle  dimension  of  100  to  1,000  A. 
Aerosols  with  an  average  particle  dimension  of  600  X  were  obtained  by 
the  atomization  [73]  of  salt  solutions  with  a  concentration  of  ~0.1%, 
separation  of  large  droplets  by  bubbling  and  subsequent  drying  of  the 
droplets.  In  modification  of  the  LaMer  generator,  known  as  the  Rapaport- 
Vaynshtoi.  generator,  was  used  [74]  for  the  production  of  monodispersed 
particles  of  dioctylphthalate  measuring  >400  X.  The  dioctylphthalate 
solution  obtained  by  mechanical  atomization  in  a  methanol  fog  was  passed 
through  a  heater,  whereupon  evaporation  of  both  components  took  place. 

Then  the  gas-vapor  mixture  was  passed  through  an  air  cooler,  whereupon 
the  vapors  condensed  on  the  walls,  or  remained  in  the  vapor  phase,  and 
the  dioctylphthalate  vapors  condensed  on  the  walls  and  on  the  nonvolatile 
sediments  of  the  vaporized  particles  that  served  as  the  seeds.  The  di¬ 
mension  of  the  particles  obtained  could  be  varied  by  changing  the  con¬ 
centration  of  dioctylphthalate  in  solution. 

The o product ion  of  monodispersed  aerosol  with  a  particle  dimension 
r  =  320  A  by  atomization  of  a  virus  culture  is  described  in  [75]. 

In  connection  with  the  development  of  research  on  the  creation  of 
the  electric  aerosol  motors  and  energy  converters  reports  have  appeared 
concerning  the  possibility  of  atomizing  liquid  cesium  in  a  strong  electrical 
field  in  a  vacuum,  forming  particles  measuring  >40  X  [76].  Although  these 
particles  are  suspended  not  in  the  gas  medium,  but  in  rarefied  cesium 
vapors,  the  stated  system  can  be  conditionally  classified  as  an  aerosol. 

It  is  interesting  to  note  that  electrostatic  atomization  of  liquids  at 
normal  pressure  does  not  produce  particles  measuring  <2,000  X  [77]. 
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CHAPTER  I! 


HDA  INVESTIGATION  METHODS 
§6.  Investigation  of  HDA  in  a  Suspended  State 

Due  to  differences  in  the  properties  of  aerosols  with  particle 

O 

dimensions  less  than  1,000  A  from  those  of  coarser  systems  investigation 
methods  also  differ.  HDA  investigation  methods  can  be  broken  down  into 
two  groups--investigations  in  suspended  state  and  investigations  based 
on  measurement  of  the  rate  of  sedimentation  of  the  particles  or  investi¬ 
gation  of  the  sediments. 

Since  the  minimum  dimension  of  particles  visible  under  the  ultra¬ 
microscope  under  normal  illumination  conditions  is  '-1,000  X  and  the 
intensity  of  diffusion  of  light  in  HDA,  according  to  the  Rayleigh  law, 
diminishes  very  rapidly  as  particle  dimension  decreases,  HDA  has  very 
little  optical  density  if  only  the  concentration  is  not  very  high,  and 
in  the  latter  case  they  very  rapidly  coagulate  and  become  coarsely  dis¬ 
persed.  Therefore,  in  order  to  use  the  optical  investigation  methods 
on  HDA  it  has  to  be  enlarged  by  means  of  condensation  on  vapor  particles 
of  some  other  substance.  For  instance,  in  the  adiabatic  counters  of 
Aytken  and  Shol'ts  [78],  long  used  by  geophysists  for  determination  of 
the  concentration  of  condensation  nuclei,  the  particles  are  enlarged 
by  condensation  of  water  vapor  on  them.  The  aerosol  is  admitted  into 
a  cylindrical  chamber  containing  humid  air.  Expansion  is  brought  about 
manually  with  the  aid  of  a  piston  and  is  accompanied  by  adiabatic  cool¬ 
ing  and  condensation  of  the  supersaturated  vapor.  The  water  droplets 
formed  by  condensation  on  the  particles  settle  out  on  a  glass  plate  and 
are  counted  with  the  aid  of  a  magnifying  glass. 


The  degree  of  supersaturation  created  in  such  counters  is  usually 
sufficient  for  the  detection  of  r.cnhygroscopic  particles  with  a  dimension 
exceeding  100  X. 

The  precise  determination  of  concentration  of  even  smaller  particles 
is  extremely  difficult,  even  by  increasing  the  degree  of  supersaturation, 
due  to  the  diffusion  sedimentation  of  particles  on  the  walls  of  the 
counter  chamber  and  the  sample  dosing  apparatus  during  the  taking  of 
the  sample  and  during  the  time  period  preceding  expansion.  Nolan  and 
Poliak  devised  a  photoelectric  adiabatic  counter  [79-81],  in  which  the 
fog  that  is  formed  does  not  settle  out  on  the  floor  and  the  concentra¬ 
tion  of  the  particles  is  determined  on  the  basis  of  light  diffusion  or 
extinction  in  the  fog  with  the  aid  of  a  photoelectric  sensor.  This  type 
of  counter  was  initially  calibrated  with  an  Aiken  counter,  with  the  result 
that  all  errors  contributed  by  the  latter  were  also  transmitted  to  the 
automatic  counter.  In  the  literature  [81,  82]  it  described  a  counter 
designed  for  calibration  purposes,  in  which  the  concentration  of  aggregated 
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particles  was  determined  by  means  of  stereo  microphotography  in  the 
suspended  state,  thanks  to  which  certain  deficiencies  of  the  Aytken 
counter  were  eliminated,  although  the  principal  source  of  errors-- 
diffusion  losses  of  fine  particles,  remained.  Early  modifications  of 
the  photoelectric  counter  had  insufficient  reproducibility  and  dis¬ 
played  a  tendency  toward  reading  drift  in  time  [83],  but  eventually 
these  deficiencies  were  eliminated  [84] .  It  is  noteworthy  that  calibra¬ 
tion  of  the  counter  depends  on  temperature  and  pressure  [85,  88].  The 
optical  density  of  the  fog  [85]  is  related  to  its  calculated  concentra¬ 
tion  n  by  the  empirical  relation  E  98  n0-1|6M0,51+,  where  M  is  the  mass  of 
the  vapor  condensed  per  unit  volume.  Other  instruments  based  on  this 
principle  are  described  by  Vertsar  [86]  and  Rich  [87]. 

A  number  of  authors  have  described  automatic  adiabatic  counters  in 
which  expansion  takes  place  at  a  rather  high  frequency--up  to  5  times 
per  second  [89-90].  Here,  the  residence  time  of  the  aerosol  in  the 
chamber  and  connecting  lines  is  reduced  and  diffusion  losses  should 
also  have  been  decreased.  However,  such  a  high  frequency  of  expansions 
inevitably  leads  to  turbulization  of  the  aerosol  flow  when  fed  into  the 
chamber,  which  can  increase  the  losses.  In  concluding  this  review  of 
adiabatic  counters  it  should  be  pointed  out  that  they  are  hardly  useful 

in  investigations  of  particles  measuring  less  than  100  X  due  to  the 
above-mentioned  losses. 

We  know  of  flow  instruments  designed  for  counting  with  aggregation, 
in  which  diffusion  chambers  are  used  for  the  purpose  of  creating  super¬ 
saturation  [91-93].  With  these,  however,  it  is  possible  to  create  only 
low  supersaturation,  sufficient  for  detecting  only  rather  large  particles, 
approximately  with  d  >  400  X,  and  the  flow  rates  of  the  aerosol  in  these 
chambers  are  quite  low. 

LaMer  developed  the  original  flow  method  of  enlarging  liquid  parti¬ 
cles  [43,  94-95].  Fogs  of  dibutylphthalate  and  sulphuric  acid  were  passed 
over  water  and  toluene,  respectively.  After  a  very  short  time  the  com¬ 
position  of  the  droplets  goes  into  equilibrium  with  the  vapor  phase  and 
the  droplets  enlarge  correspondingly  and  are  investigated  by  optical 
methods  so  that  in  addition  to  the  calculated  concentration  it  is  also 
possible  to  find  the  dimensions  of  the  original  particles  with  the  aid 
of  simple  calculations. 

In  [96]  particles  were  enlarged  by  means  of  condensation  of  water 
vapor  on  them  during  mixing  of  flows  of  HDA  and  a  hot  vapor-gas  mixture. 

The  method  of  flow  enlargment  of  HDA  proposed  by  Kogan  and  Bumasheva 
[97]  turned  out  to  be  exceptionally  successful.  In  the  KUST  enlarger 
which  they  designed  (see  Figure  1)  a  flow  of  cold  HDA  is  mixed  in  a 
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conical  diffusor  with  a  flow  of  hot  gas  containing  the  vapors  of  a  high- 
-boiling  liquid,  for  instance  dibutylphthalate.  The  flows  enter  the 
diffusor  coaxially--the  hot  vapor-gas  mixture  is  fed  to  the  center  and 
the  enlarged  HDA  is  admitted  on  the  periphery.  During  turbulent  mixing 
of  the  flows  supers at urat ion  occurs,  the  degree  of  which  can  be  control lvd 
by  the  flow  ratio,  absolute  flow  rates,  temperature  in  the  vaporizer  of 
the  instrument,  width  of  the  annular  gap  through  which  the  HDA  enters 
the  diffusor,  and  also  by  the  choice  of  the  working  liquid.  As  a  result 
of  enlargment  a  stable  monodispersed  aerosol  is  obtained  with  an  average 
particle  dimension  of  a  few  tenths  of  a  micron,  which  can  be  investigated 
by  ordinary  optical  methods.  The  remarkable  feature  of  the  instrument  is 
the  independence  of  the  dimension  of  the  enlarged  particles  on  the  counted 
HDA  concentration.  This  is  explained  by  the  fact  that  much  of  the  vapor 
is  condensed  not  into  volume,  but  on  the  walls  of  the  diffusor,  and  the 
actual  supersaturation  of  the  vapor  in  the  instrument,  and  consequently, 
the  dimension  of  the  particles  formed,  are  determined  basically  by  con¬ 
densation  on  this  surface.  In  an  improved  version  of  the  instrument, 
the  KUST-4  [63]  the  upper  concentration  limit  at  which  it  has  no  effect 
on  the  dimension  of  the  enlarged  particles  is  ~105/cm3.  Owing  to  this 
property  of  the  instrument  the  relative  HDA  concentrations  can  be  deter¬ 
mined  with  the  aid  of  a  nephelometer  without  preliminary  calibration. 

Kogan  managed  with  this  instrument  to  detect  even  individual  molecules  of 
certain  substances  [63]. 


Figure  1.  Diagram  of  KUST 
I  ns t rumen t . 
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Measurement  of  light  diffusion  by  enlarged  and  nonenlarged  aerosols 
makes  it  possible  to  determine  the  average  particle  dimension  of  l‘DA 
exceeding  400  X  if  its  counted  concentration  is  sufficiently  high  to 
ensure  measurable  light  diffusion  [97]. 

The  question  of  the  magnitude  of  particle  losses  during  enlargment 
in  the  KlIST  instrument  was  discussed  [98,  47].  For  this  purpose  a  flow 
enlarger  was  used,  in  which  HDA  was  admitted  into  the  diffusor  axially, 
rather  than  peripherally  as  in  the  KUST  instrument,  while  the  hot  gas- 
-vapor  mixture  was  admitted  along  the  walls,  with  the  result  that  dif¬ 
fusion  losses  were  eliminated.  By  comparing  the  concentrations  of  the 
enlarged  aerosols  obtained  in  both  instruments  it  was  shown  that  the 
losses  of  NaCl  particles  measuring  50,  90  and  200  X  in  the  KUST  instru¬ 
ment  were  33,  17  and  10%.  Unfortunately  the  modified  enlarger  produced 
less  nwnodispersed  aerosols  and  the  particle  dimension  depended  somewhat 
on  their  concentration. 

Aytken  [99]  first  measured  the  activity  spectrum  of  condensation 
nuclei  and  Nolan  [100]  found  that  supersaturation  required  for  detection 
of  particles  in  his  counter  could  be  used  for  determination  of  their 
dimensions.  Later  on,  in  [101],  this  method  was  used  for  evaluating  the 
particle  dimensions  of  NaCl.  The  "chemical"  diffusion  chamber,  in  which 
one  of  the  walls  is  moistened  with  water  and  the  other  with  an  HC1  solution 
[101]  can  be  used  for  creating  low  vapor  supersaturation.  With  this 
chamber  it  is  possible  to  achieve  supersaturations  in  the  0.01-1%  range 
with  an  accuracy  of  0.01%.  It  was  shown  [103,  72]  that  the  relative 
humidity  at  which  particles  of  various  mineral  salts  measuring  200  X 
begin  to  enlarge  is  close  to  that  calculated  with  the  aid  of  the  widely 
accepted  thermodynamics  relations.  By  this  method  it  is  obviously  pos¬ 
sible  to  determine  the  average  dimension  or  even  distribution  of  particles 
with  respect  to  dimensions  in  aerosols  consisting  of  one  specific  sub¬ 
stance  with  a  known  wettability  and  hygroscopic  properties.  However, 
distribution  with  respect  to  dimensions  cannot  be  found  from  the  activity 
spectrum  of  aerosol  particles  of  unknown  composition  [104].  The  theory 
of  detecting  aerosol  particles  in  a  diffusion  chamber  was  developed  by 
Buykov  [105,  106]. 

The  dimension  of  uncharged  particles  is  measured  in  [107]  on  the 
basis  of  their  capture  coefficient  of  gaseous  ions.  This  coefficient 
can  be  found  from  measurements  of  the  reduction  of  ion  concentration  when 
the  investigated  aerosol  is  admitted  into  a  bipolar  ionic  atmosphere. 

Here,  however,  it  is  necessary  to  determine  by  independent  method  the 
counting  concentration  of  the  aerosol.  Considering  evei  that  the  depend¬ 
ence  of  the  capture  coefficient  on  particle  dimension  is  not  yet  well 
defined,  the  possibilities  of  this  method  are  unquestionable. 

A  flame  counter  was  used  [108,  109]  for  determination  of  th^  counting 
concentrations  of  aerosols  from  substances  whose  vapors  produce  strong 
illumination  at  high  temperature,  for  instance,  the  salts  of  sodium  or 
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lithium.  In  this  device  the  aerosol  enters  the  colorless  flame  of  a 
hydrogen  burner.  The  passing  of  the  aerosol  particles  through  the  flame 
is  accompanied  by  flares,  which  are  recorded  by  a  photomultiplier.  A 
modification  of  this  instrument,  described  in  [110],  permits  the  record¬ 
ing  of  NaCl  particles  larger  than  400  X.  The  upper  limit  of  the  counting 
concentration  determined  by  this  method  is  limited  by  the  error  resulting 
from  coincidence,  and  is  therefore  not  very  high. 

Recently  the  narrow  angle  diffusion  of  x-rays,  by  means  of  which 
it  is  possible  to  determine  the  distribution  of  particles  with  respect 
to  dimensions,  has  been  used  for  investigation  of  highly  dispersed 
suspensions.  It  is  logical  to  expect  that  this  method  will  be  used  in 
the  near  future  for  concentrated  HDA,  for  instance,  for  investigation 
of  the  process  of  their  formation  under  stationary  conditions. 

§7.  HDA  Investigation  Methods  Based  on  Articulate  Sedimentation 

Of  the  various  methods  of  this  group,  we  will  first  discuss  methods 
of  determining  the  "electrical  mobility"  of  particles,  i.e.,  their  rate 
of  motion  in  an  electric  field  with  a  voltage  of  1  V/cm.  The  classical 
method  of  determining  electrical  mobility  [concentration)  of  gaseous 
ions  is  the  passing  of  ionized  gas  through  a  condensor  and  measurement 
of  the  current  flowing  through  the  latter  as  a  function  of  the  applied 
voltage,  and  was  first  used  by  Lanzheven  [111]  to  investigation  of  heavy 
ions,  i.e.,  charged  particles  of  atmospheric  aerosols.  Since  in  HDA  with 
d  <  600  X  there  are  practically  no  particles  with  a  charge  greater  than  1 
elemental  charge  (see  p.  62),  it  is  possible  to  calculate  from  the  electri¬ 
cal  mobility  spectrum  of  particles  in  such  aerosols  the  distribution  of 
charged  particles  with  respect  to  dimensions.  Since  the  percent  of  charged 
particles,  generally  speaking,  increases  with  their  dimension  it  is  impos¬ 
sible  by  this  method  to  find  the  distribution. of  all  particles  in  HDA 
with  respect  to  dimensions.  Thus,  this  method,  strictly  speaking,  is 
applicable  only  to  monodispersed  aerosols.  Nevertheless,  due  to  its 
simplicity  it  is  used  extensively  in  investigation  of  atmospheric  aerosols. 
In  [112]  this  method  was  used  for  analyzing  HDA  of  sulphuric  acid.  Israel 
[113]  contributed  greatly  to  the  development  of  this  method. 

In  a  similar  method  HDA  samples  are  placed  in  the  chamber  of  a 
Nolan-Pollak  counter,  in  which  there  are  flat  electrodes  with  a  potential 
difference  applied  to  them  [114].  The  concentration  of  particles  that 
remain  in  the  condensor  after  various  intervals  of  time,  i.e.,  the  kinetics 
of  diminuition  of  aerosol  concentration  is  measured,  after  which  particle 
distribution  with  respect  to  electrical  mobilities  and  the  percent  of 
uncharged  particles  are  calculated.  Here  it  is  necessary  to  disregard 
the  diffusion  sedimentation  of  uncharged  particles  on  the  condensor  walls 
which  in  the  case  of  very  fine  particles  can  lead  to  considerable  errors. 
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After  determining  the  percent  of  uncharged  particles  in  a  more  or 
less  monodispersed  HDA  under  conditions  of  a  charge  in  a  bipolar  ionic 
atmosphere  it  is  possible  to  calculate  particle  dimension.  Hie  Boltzmann 
equation  can  be  used  in  the  first  approximation  for  this  purpose:  for 
the  fraction  of  uncharged  particles  is  obtained  the  expression 

«*7n|Tn* + "«i  +" . . .  “  i  +  i  exp }  +  i exp \(^}TSJS]'  (7.1) 


where  nQ  is  the  concentration  of  uncharged  particles,  n^  is  the  concentra¬ 
tion  of  particles  with  ±1  elemental  charge,  etc.  According  to  the  latest 
data  [115]  an  exponential  formula  analogous  to  the  Boltzmann  equation 
properly  describes  the  distribution  of  charges  in  HDA  with  a  particle 
radius  >200  X.  For  more  highly  dispersed  aerosols  the  concentration  of 
charged  particles  under  stationary  conditions  can  be  calculated  only 
approximately  (see  §12) .  This  method,  advanced  by  Wright  [116]  and 
improved  by  Rich  [117]  can  be  used  for  particles  of  any  dimension,  but 
is  not  valid  for  polydispersed  aerosols.  In  the  latter  case  a  certain 
equivalent  radius  a  can  be  calculated  from  formula  (7.1).  As  shown  in 

[118],  the  ratio  of  a  to  the  arithmetic  mean  radius  is  greater  the  higher 

© 

the  polydispersion  of  the  aeroso1  and  can  be  used  as  a  measure  of  the 
latter. 

When  using  methods  based  on  measurement  of  the  stationary  distribution 
of  electrical  charges  it  must  be  borne  in  mind  that  a  certain  time  is 
required  to  achieve  such  a  state,  in  particular,  under  conditions  of 
natural  ionization  tens  of  hours  may  be  necessary  for  this  purpose  [119]. 

The  diffusion  method  of  determining  the  dimensions  of  HDA  particles 
first  used  by  Townsend  for  measurement  of  the  mobility  of  gaseous  ions, 
i;i  extremely  fruitful.  This  method  was  first  used  on  aerosols  by  Nolan 
[120]  and  Radushkavich  [121].  The  method  is  based  on  measurement  of 
the  diffusion  sedimentation  of  particles  on  the  walls  during  laminar  flow 
of  the  aerosol  through  a  channel  of  circular  or  rectangular  cross-section. 
The  passage  of  the  aerosol,  i.e.,  the  ratio  n/nQ  of  concentrations  at  the 

outlet  and  inlet  of  the  channel  can  be  precisely  calculated  theoretically 
and  the  diffusion  coefficient  D  of  particles,  related  to  their  mobility  B 
(rate  of  motion  under  the  influence  of  a  constant  force  of  1  dine)  by  the 
relation  D  ■  kTB  (see  §10  concerning  the  dependence  of  B  on  particle 
dimension)  can  be  determined  from  passage  measurements.  The  formulas  for 
diffusion  sedimentation  (and  for  the  mathematically  equivalent  problem- 
heat  transfer  in  a  channel)  were  derived  by  several  authors  in  the  form  of 
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series.  In  a  cylindrical  channel,  in  the  case 
(when  on  the  particle  concentration  profile  in 
there  is  a  plateau  with  initial  concentration) 
the  formula  [122] 


of  nonsteady  diffusion 
the  channel  cross-section 
passage  is  expressed  by 


n/n0  - 1  -  2,56  + 1,2|»  +  0,17V»  4-... 


(7.2) 


2—  — 
where  y  =  Ds/R  v,  x  is  the  length  of  the  channel,  R  is  its  radius,  v  is  the 

mean  linear  flow  rate.  For  the  steady  diffusion  state  (y  >  0.05)  [173] 


nln  o  —  0,819  exp  (—  3,657jt)  +  0,097  exp(—  22,3jt)  + 
+  0,032  exp  (—5,7^). 


In  a  plane-parallel  channel  (in  whose  cross-section  length  greatly  exceeds 
width  2h)  passage  in  both  states  is  quite  accurately  expressed  by  the 
Oe  Marcus  equation  [124] 


nln o  -  0,9149  exp  (-  1 ,885ft)  +  0,0592  (-  22,33ft)  + 

+  0,0258  exp  (-151 ,8ft),  (7*4) 


2_ 

in  which  u  »  Dx/h  v  . 

The  numerical  integration  of  the  equation  of  diffusion  in  a  plane- 
-parallel  channel,  carried  out  in  [125]  on  a  computer,  yielded  a  result 
that  agrees  satisfactorily  with  the  analytical  solution  of  De  Marcus. 

In  practice  an  aerosol  is  often  passed  through  a  system  of  parallel 
channels— an  apparatus  for  which  Rodebash  [126]  offered  the  name  "diffusion 
battery".  In  this  case  the  expression  for  u  should  contain  in  the  num¬ 
erator  a  factor  equal  to  the  number  of  channels.  The  -e  of  sedimenta¬ 
tion  of  HDA  particles  is  so  low  that  it  is  not  necessary  when  working 
with  HDA  to  place  the  battery  in  the  vertical  position.  The  above  formulas 
are  valid  for  sedimentation  of  monodispersed  aerosols,  but  actually  all 
aerosols  are  polydispeTsed  and  when  they  flow  through  a  channel  it  is 
primarily  the  finest,  -*s t  mobile  particles  that  settle  out.  Therefore, 
the  apparent  diffusion  coefficient  of  a  polydispersed  aerosol,  which  can 
be  calculated  with  the  aid  of  the  above  equations,  increases  as  the 
characteristic  parameter  y  ■  u/D  increases.  Hence,  diffusion  measurements 
are  valuable  only  when  they  make  it  possible  to  find  the  diffusion  coef¬ 
ficient  that  corresponds  to  a  certain  average  particle  dimension  for  the 
entire  dimension  distribution  of  the  particles.  At  the  present  time. 
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there  are  four  methods  used  for  processing  diffusion  experimental  data1, 
about  which  it  is  possible  to  find  the  parameters  of  particle  dimension 
distribution.  Poliak  and  Metniks  [127]  worked  out  the  so-called  exhaustion 
method,  which  consists  in  the  following.  If  we  assume  that  an  aerosol 
is  a  mixture  of  several  monodispersed  fractions  characterised  by  in¬ 
creasing  diffusion  coefficients  Dj,  D^,  ...,  and  fractions  of  calculat¬ 

ing  concentrations  p^,  ps,  ...,  p.s ,  we  can  show  that 

exp  (-  KD'/v)  -  S  ft  exP  (“  KDfik  (7.5) 

i 


where  D'  is  the  empirical  diffusion  coefficient  corresponding  to  flow 
rate  v;  K  is  a  coefficient  that  includes  the  structural  parameters  of 
the  battery  in  accordance  with  equation  (7.3)  or  (7.4).  When  v  ♦  •,  by 

expanding  the  coefficients  in  (7.S)  into  a  series  we  obtain  D'  ■ 

“1£  i*®«  *  to  v  -*•  •  corresponds  the  Man  weighted  value  D  of  the 

i 

mixture  of  the  fractions.  When  v  -*■  0  only  the  term  corresponding  to  the 
fraction  with  the  smallest  remains  in  formula  (7.5).  Therefore,  by 

extrapolating  the  experimental  function  n/n^  »  f(y)  to  v  ♦  0  we  can 
determine  Dj  and  pj.  After  subtracting  px  from  the  total  number  of  parti¬ 
cles  we  can  make  the  conversion,  determining  how  the  function  f(y)  looked 
for  the  aerosol^  from  which  was  excluded  the  first  fraction,  then  extrap¬ 
late  again  te  v  «  0,  find  02  and  p2,  etc. 

This  Mthod  has  beat  used  in  works  of  the  Irish  school  [128,  129) . 
Equation  (7.5)  was  derived  with  consideration  only  of  the  first  term  in 
the  right-hand  side  of  equation  (7.4);  more  awkward  expressions  are 
obtained  with  consideration  of  the  second  term.  Nolen  and  Scott  [130] 


iStechkina  [163]  has  just  offered  the  fifth  method  of  determining  the 
parameters  of  articulate  dimension  distribution,  based  on  the  asymptotic 
expansion  of  the  integrals  in  (7.6)  with  respect  to  the  small  parameter 
8  .  Using  the  analytical  expressions  obtained  in  this  work  it  is  pos¬ 
sible  ^  find  o  and  8  froa  measurements  of  D  for  two  values  of  v. 

8 
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delved  into  the  problem  of  the  influence  of  the  number  of  fractions 
removed  on  the  accuracy  of  tht  method  and  found  that  the  correct  rep¬ 
resentation  of  the  real  dimension  distribution  can  be  reached  by  making 
the  aerosol  of  at  least  nine  fractions.  Together  with  this,  this  method 
yields  low  accuracy,  despite  cumbersome  calculations,  due  to  the  need 
for  repeated  extrapolation.  Fuchs,  Stechkina  and  Starosel'skiy  [131] 
calculated  on  a  computer  the  values  of  the  integral  diffusive  passage 
as  a  function  of  n/nfl  ■  f(y)  for  aerosols  of  the  logarithmically  normal 

distribution  with  respect  to  dimensions  g(a)  from  the  formula  (see  7.4) 


a  n»  —  0.91 49  ^  g  (a)  exp  |—  D (a)  y j  da  + 
i 

oo 

—  0,059^  g(o)  exp  |  -D'(a)y]da  -f  . .. 

1 


(7.6) 


With  systems  of  curves  f(y)  corresponding  to  a  certain  mean  geometrical 
particle  radius  a^  in  the  range  between  10  and  1,000  A  and  to  the 

logarithms  of  standard  geometrical  deviation  from  0  to  0.4  were  plotted. 

The  Mi llikena- Knuds en  equation  was  used  for  particle  mobility  (see  $9). 

It  should  be  pointed  out  that  logarithmically  normal  distribution  was 
observed  in  aerosols  of  most  varied  origin.  The  curves  n/nQ  •  f(y) 

borrowed  from  [131]  are  illustrated  in  Figure  2.  For  practical  use  of 
this  "method  of  best  approximation"  the  calculation  curves,  constructed 
in  the  same  scale  on  tracing  paper,  are  superimposed  on  the  experimental 
data  plotted  on  coordinates  (n/nQ ,  log  y),  and  are  moved  along  the  abscissa 

until  the  best  matching  of  the  experimental  points  with  one  of  the  cal¬ 
culation  curves  is  achieved.  In  this  way  8^  is  determined.  V  is  found 

from  the  amount  of  displacement  along  the  log  y  axis  required  for  matching 
with  the  family  of  curves  corresponding  to  the  given  a^. 

From  the  theoretical  point  of  view,  Tuai  used  the  most  logical  approach 
[132],  but  it  turned  out  to  be  unsuccessful.  By  transforming  equations 
(7.3)  or  (7.4)  for  the  case  of  a  polydispersed  aerosol  and  deriving  an 
expression  analogous  to  (7.7),  he  showed  that  the  experimental  dependence 

n/i»Q  ■  f(y)  can  be  equated  to  a  converging  series  of  functions  of  the 

form  g(a)A.e‘*iy  ,  which  are  the  Laplace  transforms  of  the  desired  dis- 

*  ^  ,  V 

tribution  fraction  with  respect  to  dimensions  g(a).  Therefore,  the 
latter  can  be  found  with  the  aid  of  the  inverse  Laplace  transform. 
Theoretically,  this  method  can  be  used  for  aerosols  with  any  unknown 
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dimension  distribution.  As  shown  by  calculations  carried  out  by  Tumi  him¬ 
self  with  the  aid  of  a  computer,  however,  if  for  a  given  g(a)  f(y)  is 
calculated  and  then  if  conversion  is  made  back  to  g(e),  the  latter  will 
differ  substantially  from  the  first  g(a) .  The  error  associated  with 
this  circumstance  increases  considerably  if  the  function  f(y)  is  found 
inaccurately  as  the  result  of  imperfect  measurements.  Only  with  an 
accuracy  in  the  determination  of  passage  in  the  entire  range  of  measure¬ 
ments,  corresponding  to  errors  <0.5%  can  good  results  be  obtained  with 
the  aid  of  this  method.  Since  such  accuracy  is  usually  unattainable, 
the  Tumi  method.  Obviously,  is  suitable  only  for  estimation  of  the  width 
of  dimension  distribution  in  extremely  polydispersed  aerosols,  for 
instance  in  atmospheric  condensation  nuclei. 

Nolen  and  Scott  [130]  derived  an  approximate  relation  for  deter¬ 
mination  of  particle  distribution  parameters  in  terms  of  the  diffusion 
coefficients 


(7.7) 


Here  0  is  the  value  that_  corresponds  to  a,  D*  is  the  value  found  from 
experiment  at  the  rate  v,  a  is  standard  deviation  of  D.  By  determining 
D  for  two  values  of  v  it  is  possible  to  find  F  and  o  from  which  the  con¬ 
version  can  be  made  to  the  average  dimension  and  dispersion  of  particle 
dimensions.  Metniks  [133]  recently  compared  by  the  calculation  method 
three  methods  for  determining  the  distribution  parameters  of  particles 
with  respect  to  dimensions— exhaustion  method,  best  approximation  method 
and  the  Nolen-Scott  method.  Metniks  considered  the  Laplace  inverse 
transform  method  [132]  practically  impossible.  The  values  of  n/nQ  ■  f(y) 

found  by  the  calculation  method  from  various  particle  dimension  distri¬ 
butions  were  processed  by  all  three  methods.  Here,  the  calculations  were 
formed  for  distributions  in  which  the  number  of  particles  was  proportional 

a"1,  a'3  and  e'*,  and  also  for  normal,  logarithmically  normal,  mid  a 
certain  bimodal  distribution.  It  was  explained  that  the  curves  cal¬ 
culated  in  [131]  satisfactorily  approximated  the  function  f(y)  for  all 
these  distributions.  The  exhaustion  method  also  produced  good  results* 
but  the  Nolen-Scott  method  turned  out  to  be  suitable  only  in  the  ease 
of  normal  and  logarithmically  normal  distribution,  and  even  then  its 
accuracy  was  less  than  that  of  the  other  two  methods.  The  advantages 
of  the  best  approximation  method  over  the  exhaustion  method  are  shown 
in  the  work  of  Nolen  and  Scovt  [130]. 
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Figure  2.  Diffusive  Sedimentation  of  Polydispersed  Aerosols  in  Plane 
Parallel  Channel.  Curves ,  N,  a.  A:  1 ,  10;  1 1 ,  16;  1 1 1 ,  25;  IV,  *0; 
V,  63;  VI,  100;  VII,  160;  VIII,  250;  IX,  *00;  X,  630;  XI,  1,000; 
Symbols:  logB^  o-ai:  a  -cfc  ■  -**  *-A« 


Hie  solutions  of  equation  (7.2)  and  (7.4)  of  diffusion  from  a 
laminar  floe  were  found  on  the  assumption  that  a  parabolic  flow  profile 
is  established  in  the  vwxy  beginning  of  the  channel.  Meanwhile,  this 
profile  is  established,  for  instance,  in  a  circular  tube  only  at.  a 
distance  of  0.0S  R.Re  from  the  inlet,  where  Re  is  the  Reynolds  number. 
Thorough  examination  shows  that  diffusion  sedimentation  in  the  intake 
part  is  much  greater  than  in  the  case  of  a  parabolic  flow  profile. 

Since  the  diffusion  parameter  corresponding  to  the  end  of  the  intake 
part  is  W|n  **0.1  D/v,  where  v  is  the  kinematic  viscosity  of  the  gas, 

the  reduction  in  passage  through  the  battery  caused  by  this  condition 
increases  as  the  diffusion  coefficient  increases.  This  effect  is  of 
practical  importance  only  in  terms  of  extremely  highly  dispersed  aerosols 
with  a  particle  radius  <10  A.  Thomas  (134]  attempted  to  eliminate  the 
error  caused  by  the  intake  effect  by  measuring  the  passage  ratio  n/nQ 

in  two  batteries  having  lengths  of  47  mid  S  cm  and  by  assuming  this  ratio 
to  be  equal  to  passage  through  a  battery  42  cm  in  length.  However,  this 
approach  can  produce  considerable  error  in  the  region  of  small  u  due  to 
the  fact  that  passage  can  only  approximately  be  considered  an  exponential 
function  of  length,  then  the  influence  of  the  second  mid  third  terns  in 
equations  (7.3)  and  (7.4)  is  great  passage  through  batteries  of  length 
and  t2  is  not  oqual  to  passage  through  a  battery  of  length  l j  ♦ 


This  error  was  recognized  by  Thomas  himself  [136].  Nevertheless,  Tumi 
[137],  in  i  ier  to  find  the  points  corresponding  to  various  y,  passed 
an  aerosol  Ji  rough  the  same  battery  a  different  number  of  times.  This 
method  involves  the  same  error,  since  passage  through  a  battery  of 
length  L  is  different  from  the  square  of  passage  through  a  battery  of 
length  L/2. 

An  empirical  correction  for  the  intake  effect  was  used  in  [40]. 

This  correction  is  suitable  for  investigation  of  aerosols  with  a  small 
dimension  distribution  width.  In  the  case  of  a  monodispersed  aerosol 
the  values  of  D  calculated  experimentally  for  various  y  should  coincide. 
For  practical  purposes,  however,  D  depends  on  y  even  in  the  case  of  a 
monodispersed  aerosol,  and  this  dependence  is  attributed  to  more  in¬ 
tensive  sedimentation  of  particles  in  the  intake  section  in  comparison 
with  sedimentation  in  the  established  flow.  If  y  is  changed  not  as  the 
result  of  a  change  in  t\  e  velocity  cf  flow  through  the  battery,  but  by 
using  a  battery  of  difterent  length,  then  the  relative  effect  of  the 
intake  effect  will  diminish  as  the  battery  length  is  increased  and  the 
absolute  value  of  the  effect  will  be  the  same  for  measurements  for  all 
y.  Therefore,  the  intake  effect  can  be  evaluated  by  comparing  passage 
measured  in  batteries  of  various  lengths.  As  seen  in  Figure  2  the  limbs 
of  the  calculation  curves  for  log  $  =1.0  and  1.26  for  values  n/n.  > 

g  0 

>  0.45  coincide.  This  permits  evaluation  of  the  magnitude  of  the  intake 
effect  even  for  aerosols  possessing  some  polydispersion,  but  in  this 
case  it  is  necessary  to  compare  passages  that  exceed  0.45.  It  was  ex¬ 
plained  in  [40]  that  this  correction  is  notable  only  for  extremely  highly 
dispersed  aerosols;  nevertheless,  even  for  aerosols  with  a  =  7-10  X  the 
correction  factor  for  passage  through  very  short  batteries  (through 
which  passages  were  measured  in  the  range  of  0.85-0.9)  did  not  exceed 
i .05-1 .02 1 . 

Experimental  determination  of  the  function  f(y)  is  made  by  measuring 
the  concentrations  of  consolidated  aerosol  before  and  after  the  battery. 
In  the  works  of  Nolen  and  Poliak  and  coworkers  their  adiabatic  counter 
was  used  for  this  purpose,  and  in  the  works  of  the  authors  of  the  present 
review  [40,  42,  98J  a  KUST  flow  enlarger  [97]  (see  pp,  20-21)  was  used  in 
combination  with  a  photoelectric  mephelometer.  Sometimes  the  radiography 


equation  of  the  diffusion  of  particles  to  the  vails  of  a  circular 
channel  during  the  process  of  establishing  a  paraboiic  profile  was 
recently  solved  by  calculation  methods  [163].  The  additional  sediments 
tion  cf  particles  in  the  intake  section  was  with  a  Schmidt  number 
Sc  *  v/D  -  10  and  5%  for  Sc  =  0.9.  This  agrees  satisfactorily  with  the 
* mpirical  correction  introduced  in  [40]. 


of  filters  through  which  is  passed  an  aerosol,  passed  and  not  passed 
through  a  battery,  is  used  for  determination  of  n/n^  during  work  with 

radioactive  aerosols  [138,  139].  In  [139]  the  measurements  of  the  total 
radioactivity  of  the  filters  were  compared  with  the  concentrations  of 
the  particles  measured  by  an  adiabatic  counter  before  and  after  the 
battery,  and  it  was  shown  that  the  radiographic  method  yields  two  to 
three  times  lower  values  of  D.  This  is  understandable,  since  the  radio¬ 
activity  of  activated  decay  products  of  the  emanations  of  aerosol  parti¬ 
cles  increases  with  dimension.  The  distribution  of  activity  in  the 
sediment  on  the  walls  of  the  diffusion  channel  was  measured  [140]  in 
order  to  determine  the  diffusion  coefficient  of  the  decay  products  of 
radon.  The  diffusion  coefficients  were  calculated  from  the  measurement 
results  with  the  aid  of  Berezhniy's  solution  [141]  of  the  equation  that 
describes  sedimentation  of  products  formed  inside  the  diffusion  channel 
itself. 

When  designing  diffusion  batteries  it  is  necessary  to  bear  in  mind 
that  it  is  impossible  to  create  a  "universal"  battery  suitable  for  the 
entire  range  of  dimensions  of  HDA,  since  their  diffusion  coefficient 
varies  fron.  about  10" 1  to  about  10“ 5  cm2/sec.  To  find  the  required 
values  of  u  in  the  entire  dimension  range  it  would  be  necessary  to  vary 
by  four  orders  the  flow  rate  through  the  batter,  which  is  inconceivable, 
since  at  high  velocities  the  hydraulic  resistance  increases  greatly  and 
at  small  velocities  the  air  attributed  to  convection  currents  increases. 
Calculations  and  personal  experience  of  the  authors  of  this  review  in¬ 
dicate  th^t  it  is  advisable  to  use  slit  batteries  for  particles  larger 
than  200  A,  and  for  dimensions  of  the  order  of  10  X- -batteries  made  of 
cylindrical  channels.  A  set  of  batteries  of  different  length,  consist¬ 
ing  of  eight  tubes,  each  one  6  mm  in  diameter,  was  used  in  [198]  for 
investigation  of  aerosols  in  the  particulate  dimension  range  of  14-60  X, 
and  for  particles  measuring  60  to  200  X,  a  set  of  60-channel  batteries 
consisting  of  tubes  measuring  1.4  mm  in  diameter.  Generally  speaking 
it  is  not  advisable  to  use  a  large  number  of  small  diameter  tubes  for 
batteries,  but  rather  a  few  wide  tubes,  since  in  this  case  there  is  a 
greater  possibility  of  expansion  of  the  range  of  values  y  of  the  battery 
due  to  change  in  the  flow  velocity.  The  use  of  a  set  of  batteries  of 
different  lengths  for  determination  of  the  function  f(y)  makes  it  pos¬ 
sible  to  exclude  error  in  the  measurement  of  flow  velocity,  since  all 
measurements  are  made  at  the  same  velocity,  and  facilitates  the  intro¬ 
duction  of  the  correction  factor  for  the  intake  effect. 

As  regards  errors  that  may  occur  with  the  diffusion  method,  instead 
of  the  above-mentioned  intake  effect  we  should  point  out  the  acceleration 
of  sedimentation  of  charged  particles  on  the  channel  walls  due  to  mirror 
forces.  Calculations  [142,  143]  indicate  that  this  effect  is  negligibly 
small  if  the  particles  carry  no  more  than  one  elemental  charge,  but  can 
be  considerable  in  a  case  of  large  charges.  In  checking  the  diffusion 
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method  with  the  aid  of  a  monocispersed  aerosol  obtained  from  polystyrene 
latex,.  Migo  and  Uiffen  [144]  showed  that  the  correct  results  can  be  ob¬ 
tained  only  when  grounded  metallic  batteries  are  used.  Latex  particles 
undoubtedly  could  carry  large  charges.  It  was  revealed  in  [40]  that  if 
the  investigated  aerosol  had  a  higher  temperature  than  the  battery  wall, 
precipitation  of  small  u  was  greater  than  in  isothermic  conditions, 
obviously  as  a  result  of  thermophoresis  and  convective  transfer.  When 
working  with  exceedingly  highly  dispersed  aerosols  it  is  important  to 
keep  the  losses  in  the  connections  constant  during  measurements  in  bat¬ 
teries  of  different  lengths.  For  this  purpose  the  stopcocks  were 
omitted  in  [40]  and  the  difference  in  the  length  of  the  battery  was 
compensated  by  moving  the  KUST  instrument  on  a  dolly.  The  latter  measure 
is  necessitated  by  the  fact  that  according  to  the  observations  of  the 
authors  and  also  of  Nolen  and  Kaffel  [145]  the  sedimentation  of  particles 
takes  place  much  more  vigorously  in  curved  rubber  tubes  than  in  straight 
ones  due  to  the  electrical  charges  created  by  the  bending  of  the  tubes. 

In  addition  to  Migo  and  Uiffen,  Thomr  [135]  also  checked  the  dif¬ 
fusion  method,  comparing  the  results  of  the  optical  and  diffusion  deter¬ 
minations  of  particle  dimensions  in  the  1,000--3,000  X.  range  and  obtaining 
satisfactory  agreement.  The  slight  discrepancies  noted  can  be  attributed 
completely  to  incorrect  compensation  of  the  intake  effect.  For  coarser 
particles  the  diffusion  method  is  unsatisfactory  due  to  an  increase  in 
sedimentation,  which  cannot  be  eliminated,  even  by  using  vertical  batteries. 

Essentially,  the  diffusion  method  can  be  used  not  only  in  channels, 
but  also  in  any  system  of  bodies  that  permit  theoretical  calculations  of 
the  diffusive  sedimentation  of  particles  from  a  flow,  for  instance,  in 
a  model  of  fibrous  filters  consisting  of  a  system  of  parallel  cylindrical, 
randomly  arranged  fibers  oriented  perpendicular  to  the  flow.  The  passage 
of  HDA  through  such  a  model  is  expressed  by  the  formula 

«  «o  exp (— 25t// /“/?,), 

where  t,  is  the  fraction  of  the  volume  filled  by  the  fibers;  H  is  the 
thickness  of  the  model,  is  fiber  radius,  n  is  the  coefficient  of  capture 

of  the  particles  by  the  fibers,  expressed  by  the  theoretical  formula  [146] 

-  2,9/\  'Pc  ".t>24Pe~i. 


where  Pe  =  '2.RX  D.  A’*=  1/2  In  J-  + ;4--i,.75  ,  V  is  the  average  flow  velocity 

in  the  model.  Formula  (7.9)  is  in  satisfactory  agreement  with  experimental 
data  [147]  on  the  passage  of  monodispersed  HDA  with  a  =  18-80  A  thorugh 
model  filters  (see  Figure  3).  The  diffusion  coefficients  of  the  particle 
were  determined  in  these  experiments  with  the  aid  of  diffusion  batteries. 
The  advantage  of  this  version  of  the  diffusion  method  is  the  absence  of 
input  effects. 
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Figure  3.  Capture  Coefficients 
of  Diffusive  Sedimentation 
of  Particles  in  Model  Filter 
with  Checkerboard  Distribu¬ 
tion  of  Cylinders: 

Theoretical  curves 

1— * ■0.01:  1I-E  ■«.<»;  1)1— {■0,135; 

IV— E  ■0,27. 

Mean  particle  radius 

:  a.  X:  /—IS;  2-18:  J-60;  4—60; 

5-70;  6—41;  7-70;  *■  55 


In  [148]  the  diffusion  method  was  combined  with  the  determination 
of  the  number  of  charged  particles:  for  the  latter  purpose  the  diffusion 
battery  was  used  as  a  plane-parallel  condenser.  Measurement  of  the  per¬ 
cent  of  charged  particles  in  stationary  conditions  makes  it  possible  to 
determine  the  "equivalent"  dimension  of  aerosol  particles  (see  p.  29), 
from  diffusion  measurements  it  is  possible  to  determine  some  average 
particle  dimension  and  the  ratio  between  these  two  average  values  can 
characterize  qualitatively  polydispersion.  The  general  idea  of  this 
method  is  valid,  but  since  the  values  of  D,  measured  by  the  diffusion 
method,  depend  on  y,  the  result  thus  obtained  will  depend  both  on  the 
parameters  of  the  aerosol  and  on  the  measurement  conditions. 

In  some  works  [149,  150]  the  so-called  "static  diffusion  method" 
was  used.  In  this  method  the  change  in  the  concentration  of  an  aerosol 
at  rest  in  a  cylindrical  or  spherical  vessel  in  time  is  traced.  For 
particles  with  a  dimension  of  about  10  A  this  method  is  unsuitable,  since 
such  particles  settle  out  in  time  not  exceeding  a  few  seconds,  and  for 
coarser  particles  it  becomes  less  accurate  than  the  dynamic  method  due 
to  errors  attributed  to  convection. 


Labeyri  [151]  offers  a  method  for  investigation  of  aerosols,  in¬ 
cluding  highly  dispersed  aerosols,  based  on  their  radioactivation  by 
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adding  radon.  After  settling  out  of  the  decay  products  on  aerosol  parti¬ 
cles  these  latter  can  be  trapped  in  a  filter,  after  which  the  filter  is 
subjected  to  radiographic  analysis  and  the  number  of  particles  counted. 

In  conclusion,  we  will  examine  HDA  investigation  methods  based  on 
their  preliminary  sedimentation.  The  method  of  investigating  highly 
dispersed  sediments  does  not  depend  on  the  method  by  which  they  are 
obtained.  Therefore,  it  is  not  useful  here  to  examine  in  detail,  for 
instance,  the  electron  microscopic  methods,  since  they  are  disucssed 
in  an  abundance  of  special  literature.  We  will  confine  ourselves  here 
only  to  examination  of  specific  moments  that  occur  during  investigation 
of  aerosol  sediments.  The  chief  condition  here  is  the  representation 
of  the  sample,  i.e.,  the  method  of  particle  sedimentation  need  not  be 
selective.  Usually  thermal  precipitators,  designs  of  which  insure 
sampling  of  representative  samples  were  described  in  [152]  and  [153], are 
employed  for  sampling.  The  most  representative  samples  are  obtained 
with  thermal  precipitators  with  a  moving  wire  [154].  Sometimes  electro¬ 
precipitators  are  also  used  for  sedimentation.  Since  classification  of 
particles  takes  place  to  some  extent  in  these  instruments  during  the 
motion  of  aerosol  along  the  sedimentation  electrode,  it  is  necessary 
to  analyze  samples  taken  at  various  points  of  the  electrode.  The 
electroprecipitators  are  used  almost  exclusively  for  investigation  of 
aerosol  sediments  by  means  of  the  BET  method,  x-ray  diffraction  analysis 
and  other  methods . 

An  ingenious  electrostatic  precipitator  design  is  described  in  [155]. 
A  jet  of  charged  HDA  in  coronal  discharge  is  admitted  to  a  laminar  flow 
between  plane-parallel  or  coaxial  electrodes.  Classification  during 
sedimentation  is  a  result  of  the  fact  that  the  rate  of  displacement  of 
the  particles  to  sedimentation  electrode  is  proportional  to  the  charge 
and  inversely  proportional  to  the  resistance  coefficients  of  a  particle. 

Electron  microscopy  is  an  extremely  important  HDA  investigation 
method,  since  it  permits  direct  determination  of  the  dimensions  and  shapes 
of  aerosol  particles  and  also  of  the  count  concentration  of  HDA.  However, 
the  use  of  electron  microscopy  for  analysis  of  HDA  involves  certain  dif¬ 
ficulties  which  increase  by  measure  of  reduction  of  particle  dimension, 
caused  by  the  fusing  and  evaporation  of  particles  by  the  electron  beam, 
and  also  by  reduction  in  image  contrast  due  to  increased  particle  trans¬ 
parency.  The  known  electron  microscopic  methods--shading  of  particles 
by  vacuum  sublimation  of  heavy  metals,  shading  of  replicas,  sandwich 
method  [156,  157],  are  of  limited  benefit,  since  during  analysis  of 
extremely  fine  particles  the  thickness  of  the  shaded  layer  becomes 
commensurate  with  the  particle  dimension.  Therefore,  it  is  not  surpris¬ 
ing  that  although  many  authors  have  analyzed  sediments  of  particles  with 

a  dimension  approaching  the  limit  resolution  of  the  electron  microscope 
[32,  44,  45],  as  illustrations  for  their  articles  they  prgfer  to  use 
photographs  of  particles  measuring  between  100  and  1,000  A,  It  is 


-34- 


I 


pointed  out  in  [47]  that  it  is  impossible  to  take  sharp  photographs  of 
NaCl  particles  measuring  less  than  200  A,  but  in  [158]  it  is  not  explained 
that  Pt  and  Ni  particles  measuring  less  than  100  a  are  either  trans¬ 
parent  to  electrons,  or  evaporate  rapidly  under  their  influence.  The 
authors  of  the  present  review  are  not  successful  with  the  aid  of  the 
electron  microscope,  with  a  resolution  limit  of  15  A,  to  investigate 
the  structure  of  MoO^  aggregates  consisting  of  individual  particles 

measuring  about  100  A  (as  was  done  by  the  BET  method) .  The  problem  of 
microscopy  of  liquid  particles  should  be  considered  unsolved.  In  those 
cases  where  it  is  possible  to  gather  weighable  quantities  of  aerosol 
sediment,  the  average  particle  dimension  can  be  calculated  from  the 
specific  area  of  the  sediment,  determined  by  any  of  the  common  methods. 

If  after  coagulation  the  particles  did  not  stick  together  or  unite,  but 
loose  aggregates  were  formed,  the  dimension  of  the  former,  noncoagulated 
aerosol  can  be  determined  from  the  specific  surface  of  the  sediment, 
although  the  formation  of  solid  aggregates  and  particularly  their  ag¬ 
glomeration  can  lead  to  overestimation  of  particle  dimension  calculated 
in  this  way.  In  [65,  66]  the  results  of  the  BET  method  were  compared 
with  the  electron  microscopic  analyses  and  it  was  established  that  in 
determination  of  the  particle  dimension  of  A1203,  Ti02,  Zr02,  Fe203  and 

SiO,,  sediments  obtained  by  the  interaction  of  the  vapors  of  the  chlorides 

of  these  metals  with  water  vapor  at  1, 600-2, 000°C,  both  methods  yield 
similar  results. 

The  results  of  both  methods  also  agree  satisfactorily  [37]  with 
respect  to  sediments  of  aerosols  of  oxides  of  the  metals  obtained  in  an 
electric  arc.  On  the  other  hand,  the  reduction  has  beer  noted  [47,  48] 
in  the  specific  area  of  the  sediments  of  NaCl  and  nickel  HD A,  which  can 
be  a  serious  source  of  errors,  particularly  if  there  is  a  lengthy  time 
lapse  between  sediment  production  and  its  analysis. 

Membrane  filters  are  also  used  fcr  HDA  'sis  [159-161].  Cartright 
[162]  for  electron  microscopy,  used  the  metht  A  producing  relics  from 
particles  precipitated  on  membrane  filters.  For  analysis  of  radioactive 
aerosols  radiography  of  filters,  including  membrane  filters,  is  often 
employed.  Sometimes  membrane  filters  are  used  for  detection  of  highly 
dispersed  particles  possessing  icing  activity.  For  this  purpose  the 
filter,  after  the  settling  out  on  it  of  these  particles,  is  immersed 
in  a  super-cooled  sugar  solution  and  ice  crystals  form  on  the  active 
particles.  Here,  preliminary  determination  of  the  icing  background  of 
the  filter  is  required,  since  each  filter  contains  two  to  three  centers 
that  possess  activity  of  this  type. 

In  concluding  our  discussion  of  HDA  investigation  methods  we  might 
point  out  that  the  present  methods  insure  reliable  determination  of  all 
parameters  of  HDA  consisting  of  solid  particles. 
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CHAPTER  I  I  I 


HDA  PROPERTIES 


d8.  Transfer  Process  in  HOA.  General  Comments 


The  processes  of  mass,  energy,  pulse  and  charge  transfer  from  a 
medium  to  particles  and  back  are  of  fundamental  importance  in  aerosol 
physics.  The  character  of  these  processes  is  governed  by  the  Knudsen 
number  Kn  *  Z'a.  When  Kn  -*  0  these  processes  are  described  accurately 
!>y  the  equations  of  thermoconductance,  diffusion  and  hydrodynamics  of 
continua.  When  Kn  >  «•  particles  do  not  violate  Maxwellian  distribution 
of  the  velocities  of  gas  molecules  moving  toward  the  particles,  and  for 
molecules  moving  away  from  the  particles  the  distribution  function  is 
relatively  simple  to  find.  Therefore,  transfer  processes  in  this  "free 
molecular"  mode  are  expressed  through  simple  formulas,  derived  from 
molecular  kinetic  theory  of  gases. 

Otherwise,  we  are  dealing  with  transfer  at  intermediate  Knudsen 
numbers:  here  the  theory  is  fraught  with  tremendous  mathematical  dif¬ 
ficulties.  For  quite  a  long  time  the  "stitching"  method  was  used  to  solve 
this  problem:  it  was  assumed  that  transfer  in  the  immediate  vicinity  of  a 
particle  occurs  in  the  free  molecular  mode  and  starting  at  some  distance 
of  the  order  of  l  from  the  surface  of  the  particle,  according  to  the 
equations  of  transfer  in  continua.  Although  in  certain  cases  a  stitching 
method  produced  rather  good  agreement  with  experimental  data  it  essentially 
cannot  be  accurate.  Nevertheless,  this  method  alone  is  still  used  for 
charge  transfer. 


The  rigorous  calculations  of  the  rate  of  transfer  for  intermediate 
Knudsen  numbers  can  be  based  only  on  the  solution  of  the  Boltzmann  kine¬ 
matic  equation,  and  this,  in  turn,  is  possible  only  after  some  simplifi¬ 
cation,  (.linearization)  of  the  latter.  The  simplest  linearization  method 
advanced  independently  by  Welander  and  Bkhatnagar  and  coworkers,  and 
described  in  the  ensuing  section,  permitted  the  successful  solution  of 
several  problems  of  theoretical  physics.  Brock  and  Willis  used  this 
method  in  conjunction  with  "Knudsen  iteration"  for  transfer  phenomena  in 
aerosols.  Knudsen  iteration  consists  in  expansion  of  distribution  func¬ 
tion  f  into  a  scries  by  powers  Kn-1 


where  f()  is  the  Maxwell  distribution  function.  By  substituting  this 

expression  into  the  linearized  Boltzmann  equation  the  above-stated  authors 
were  able  to  solve  it,  retaining  only  the  first  two  terms  in  the  expansion. 


and  thereby  found  the  function  fj,  i.e.,  the  first  approximation  to  the 

distribution  function,  and  hence  determined  the  transfer  rate  or  flow  of 
the  corresponding  value  to  a  spherical  particle.  The  formulas  thus  de- 
rived  or  flows,  have  the  form 


=  #t/(l  +  XKn-1) 


uhere  4>k  is  flow  in  the  free  molecular  mode  and  X  is  sum  constant.  It 
is  noteworthy  that  formulas  of  this  form  can  be  accurate  only  when 
Kn'1  <  1. 

Considerable  progress  has  been  made  in  recent  years  in  the  theory 
of  transfer  processes  at  intermediate  Knudsen  numbers,  and  this  progress 
has  been  made  by  researchers  having  very  little  to  do  with  HDA.  Sakhin, 
working  on  the  theory  of  transfer  and  absorption  of  neutrons  in  nuclear 
reactors,  made  an  accurate  calculation  of  the  rate  of  mass  transfer  to 
a  sphere  for  all  Knudsen  numbers  using  the  so-called  integral  equation. 
Since  transfer  processes  at  high  and  intermediate  Knudsen  numbers  are 
of  great  importance  in  the  aerospace  industry,  this  problem  attracted 
the  attention  of  many  theoreticians.  Using  the  variation  method  developed 
by  Cherchin'yanin  Pagani,  they  were  able  to  make  an  accurate  calculation 
of  the  heat  flux  to  a  spherical  particle  and  the  resistance  of  the  medium 
to  its  motion  for  all  Knudsen  numbers,  and  the  calculation  results  agree 
very  satisfactorily  with  the  experimental  data. 

Naturally,  in  works  related  to  the  aerospace  industry  a  great  deal 
of  attention  is  devoted  to  transfer  processes  at  high  particle  velocities 
in  relation  to  the  medium  and  high  temperature  and  concentration  gradients. 
In  the  present  review  we  will  examine  only  the  case  of  low  values,  in 
particular,  low  velocities  in  comparison  with  the  thermal  velocity  of 
molecules. 

We  should  mention,  finally,  that  in  the  theory  of  the  processes  of 
transfer  to  particles  it  is  inevitably  necessary  to  proceed  on  the  basis 
of  the  concept  of  mean  free  path  of  gas  molecules.  At  the  same  time 
this  concept  is  generally  not  used  in  modem  theory  of  transfer  processes 
in  a  gas  volume,  chiefly  because  it  is  somewhat  lacking  in  definition. 

In  essence,  in  the  theory  of  transfer  to  particles  l  should  be  regarded 
not  as  the  mean  distance  traveled  by  molecules  between  consecutive  col¬ 
lisions,  but  as  the  mean  effective  free  path  l  which  may  differ  for 

the  various  transfer  phenomena.  Thus,  l ^  for  mass  transfer  is  the 

average  distance  that  separates  a  molecule  from  the  initial  position  at 
the  moment  when  its  velocity  vector  becomes  independent  of  the  initial 
velocity  vector.  In  heat  transfer  l  ^  is  the  distance  at  which  the 
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molecule  acquires  the  average  temperature  of  the  surrounding  gas,  and 
in  moment  transfer  is  the  distance  at  which  the  molecule  acquires  the 
average  (i.e.,  hydrodynamic)  gas  velocity.  Without  going  into  detailed 
analysis  of  this  question,  it  may  be  stated,  for  instance,  that  in  view 
of  the  incomplete  transfer  of  energy  corresponding  to  the  inner  degrees 
of  freedom  of  molecules  in  collisions,  ^e££  for  heat  transfer  may  be 

greater  than  for  moment  transfer.  Up  to  now  this  important  fact  has  been 
taken  into  account  only  in  the  work  of  Lord  and  Harbor  [165]. 

59.  Hass  Transfer  in  HDA:  Evaporation  and  Condensation  Growth  of  Particles 

Evaporation  and  condensation  growth  of  particles  play  an  extremely 
important  role  in  the  formation  and  life  of  aerosols  and  have  been  the 
object  of  numerous  theoretical  and  experimental  investigations.  There¬ 
fore,  we  will  examine  these  processes  at  some  length.  When  Kn  ■*-  0  they 
are  accurately  described  by  the  division  equation.  Ordinarily  an  extremely 
idealized  process--evaporation  (or  growth)  of  a  single  spherical  particle 
at  rest  in  an  infinite  volume,  not  interacting  with  a  gas  particle,  is 
examined,  The  problem  is  especially  simplified  if  the  partial  pressure 
of  the  vapor  is  substantially  less  than  that  of  the  gas.  With  this  as¬ 
sumption,  collisions  among  vapor  molecules  can  be  disregarded  and  the 
process  can  be  assumed  to  be  quasistationary,  since  the  relaxation  time 
of  the  vapor  concentration  distribution  around  the  particle  is  much  less 
than  the  total  particle  evaporation  time  [166].  In  the  ensuing  discussion 
we  will  consider  this  case  alone. 

At  a  sufficiently  great  distance  from  a  particle  vapor  transfer  takes 
place  by  diffusion.  The  difference  of  equation  of  stationary  vapor  dif¬ 
fusion  from  a  sphere  or  to  it  has  the  form 

P  2p“,/’  (9.1) 

where  n  is  vapor  concentration  at  distance  0  from  the  center  of  the  sphere; 

1  is  the  constant  found  from  the  boundary  condition  on  the  surface  of  the 
sphere . 

We  will  examine  the  two  simplest  cases --evaporation  of  particle  in 
a  medium  net  containing  vapor  and  condensation  on  a  "black",  i.e.,  one 
that  absorbs  all  vapor  molecules  that  strike  it,  sphere.  Since  in  these 
cases  n  is  proportional  to  the  concentration  n  of  the  saturated  vapor 

3 

of  the  particle  or  to  the  concentration  of  the  vapor,  respectively, 
for  0  -  «,  by  integrating  (9.1)  we  obtain  in  the  former  case 
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n  —  —  S/ij/p, 


(9.2) 


and  in  the  latter, 

n  n *  0  “  &/*>)•  •  (9.3) 

where  6  is  a  constant,  positive  for  condensation  and  negative  for  evapora¬ 
tion.  Diffusive  flow  to  the  sphere  (and  from  it)  is  expressed  by  the 
formulas 

d^^DSfloo  and  G>—4s Din,.  f0  ... 


If  the  vapor  concentration  at  infinity  during  evaporation  is  not  equal  to 
zero,  but  to  n#,  then  it  is  necessary  in  all  corresponding  formulas  to 

replace  n  by  n  -  n  and  n  by  n  -  n  .  The  analogous  substitution  must 

oo  S  5  oo 

be  made  for  condensation  on  a  particle  with  a  vapor  pressure  not  equal  to 
zero,  but  to  n  . 

3 

All  these  equations  are  valid  for  any  Knudsen  number,  but  only  at  a 
great  distance  from  the  sphere.  When  Xn  0,  however,  vapor  transfer  is 
of  a  diffusive  character  all  the  way  to  the  very  surface  of  the  sphere, 
the  boundary  condition  for  evaporation  will  be  np-a  •  ns,  and  for  con¬ 
densation  on  a  black  sphere  n  *  0.  Thus,  in  the  case  6  •  ±a  and  (9.4) 
also  converts  to  the  familiar  Maxwell  formulas 

d>,  —  4eDon00  or  (9.5) 

We  will  note  that  when  Kn  -*■  0  the  mass  transfer  rate  is  independent  of 
the  evaporation  (condensation)  coefficient  a  .  Otherwise,  according  to 

recently  published  data  (167],  c  is  close  to  unity  for  all  liquids. 

V 

When  Kn  ■*  •  the  rate  of  evaporation  of  a  particle  (when  n  ■  0)  and 
condensation  on  it  (ns  *  0)  is  expressed  by  the  Knudsen  formula  for 

evaporation  in  a  vacuum 

or  (9.6) 


where  Vy  is  the  average  velocity  of  thermal  motion  of  vapor  molecules. 

The  indices  c  and  k  denote  transfer  in  the  continuous  and  in  the  free 
molecular  modes.  We  will  note  that 
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(9.7) 


®f/Oj  » -j  Kna,1 


As  already  pointed  out,  the  first  attempts  to  solve  the  problem  of 
the  evaporation  of  particles  at  intermediate  Knuds en  numbers  amounted  to 
the  stitching  together  of  the  free  molecular  and  diffusion  flows.  In  a 
number  of  works  (168-170]  the  flows  were  stitched  together  directly  on 
the  surface  of  the  particle,  i.e.,  it  was  assumed  that  the  diffusion 
equation  was  valid  ail  the  way  to  the  very  surface.  The  formula  obtained 
by  this  method 


d»  ---  U  Jh'n  3af)  r  =  (1  f  3Kn~,a,  4) 


(9.8) 


gives  the  correct  functional  dependence  of  *  on  Kn,  but  considerably 
overstated  coefficients  for  Kn  and  Kn"1. 

In  the  "boundary  sphere"  method  [166]  the  stitching  was  made  on  the 
surface  of  a  sphere  surrounding  a  particle,  consisting  of  the  surface  of 
the  particle  at  a  distance  8 l.  by  which,  on  the  average,  evaporating 
vapor  molecules,  having  traveled  distance  l,  are  separated  from  the 
particle  (when  Kn  •*  0  8  «  1/2  and  when  Kn  -*•  •  8  ■  1).  The  formula  for 
this  is 


-j  Ka/a,+  l/(t+PKa) 


When  Kn  <  1  and  a  «  1 
c 


nhr* 


>+i*» 


(9.9) 


which  is  quite  close  to  the  accurate  expression,  although  when  Kn  ►  1  an 
incorrect  functional  dependence  is  obtained: 


<S>  -  0,1(1  +  3Kn-*/4V 


(9.10) 


thus,  we  see  that  the  stitching  methods  do  not  produce  the  proper  results 
he  will  also  note  that  the  "universal"  interpolation  formula  of  Sherman 
(1?1),  which  he  proposed  for  several  transfer  processes 


0  -  *»)  •  C  -r  ♦,). 


(9.11) 


cotncido  with  formula  (9,7)  when  *  l. 
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The  rigorous  solution  of  the  prol ,em  of  particle  evaporation  at 
intermediate  Kn  numbers  is  possible  only  on  the  basis  of  the  Boltzmann 
kinetic  equation.  Under  the  above  assumption  of  low  vapor  concentration 
in  comparison  with  the  concentration  of  the  gas,  the  distribution  func¬ 
tion  of  the  molecules  of  the  latter  will  be  only  slightly  perturbed  by 
collisions  by  vapor  molecules  and  it  can  be  assumed  to  satisfactorily 
approximate  the  Maxwellian  concentration  all  the  way  to  the  very  surface 
of  the  particle.  However,  for  vapor  molecules  the  Maxwellian  distri¬ 
bution  occurs  only  at  a  great  distance,  compared  to  l,  from  the  surface 
of  the  particle.  In  the  case  of  a  "black"  sphere,  at  its  very  surface 
the  projections  of  the  velocities  of  all  vapor  molecules  onto  the  normal 
to  the  surface  are  directed  toward  the  latter,  i.e.,  here  the  distribu¬ 
tion  function  is  described  by  a  single  "half"  of  Maxwellian  distribution 
directed  toward  the  sphere.  By  measure  of  distance  from  the  sphere  this 
asymmetry  of  distribution  gradually  vanishes.  Vapor  transfer  at  a  great 
distance  from  the  sphere  is  purely  diffusive,  i.e.,  determined  by  the 
vapor  concentration  gradient.  By  measures  of  convergence  upon  the  sphere 
the  diffusive  vapor  flow  acquires  the  flow  attributed  to  the  asymmetry 
of  the  distribution  function  and  only  this  latter  flow  remains  at  the 
surface  of  the  sphere. 

The  chief  difficulty  in  the  solution  of  the  kinetic  equation  is 
the  collision  integrals  it  contains.  These  integrals  determine  the 
change  in  the  number  of  molecules  in  a  given  elemental  volume  and  which 
possess  a  given  velocity  vector,  due  to  their  collisions  with  other 
molecules  that  happen  to  be  in  this  volume.  This  difficulty  is  attri¬ 
buted  principally  to  the  fact  that  beneath  the  integrand  are  the  products 
of  the  distribution  functions.  The  collision  integral  was  replaced  in 
1172,  173]  ty  the  simple  expression  v(fg  -  f)  dt,  where  dT  is  an  element 

of  space,  f  is  the  distribution  function  in  this  space,  (Q  is  the  Maxwellian 

distribution  function,  v  is  thennumber  of  molecular  collisions  per  unit 
time  in  space  (in  our  case,  collisions  among  molecules  of  the  vapor  and 
gas).  This  method  was  used  successfully  for  the  solution  of  several 
physical  problems,  in  particular,  by  tfelander  for  calculating  heat 
transfer  from  a  heated  cylinder  of  radius  commensurate  with  l. 

In  the  problem  at  hand,  this  method  takes  on  a  graphic  physical 
meaning:  instead  of  examining  in  detail  the  changes  in  the  distribution 
function  of  the  vapor  molecules,  caused  by  (heir  collisions  with  the  gas 
molecules  (according  to  what  was  said  above,  eolUaiont  between  vapor 
molecules  can  be  disregarded),  we  will  assume  that  during  these  col¬ 
lisions  the  vapor  molecules  acquire  Maxwellian  (isotropic)  velocity 
distribution,  which  gas  molecules  possess.  It  should  be  stressed  that 
this  assumption  is  valid  only  when  the  mass  ratio  of  the  vapor  and  gas 
molecules  is  m^m  *  0.  Otherwise  there  may  be  some  persistence  in  the 
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velocities  of  the  vapor  molecules  and  their  velocity  distribution  may  not 
be  completely  isotropic  after  collisions. 

Brock  [174,  175]  used  this  method  of  linearization  of  the  kinetic 
equation  in  combination  with  the  Knudsen  iteration  for  the  theory  of  mass 
transfer  for  Kn  >  1 .  Brock  gives  only  the  final  formula  without  the 
calculations.  The  correction  factor  of  this  formula  contains  the  mean 

free  path  of  gas  molecules  l  .  Since  all  formulas  presented  in  this 

S 

section  contain  the  mean  free  path  l y  of  vapor  molecules,  which  appear 

as  an  impurity  in  the  gas,  we  transformed  Brock's  formula  with  the  aid 
of  the  we li -known  gas  kinetics  formulas 


lK  -  iy2::iy$~\ 


!,  -  ir/?  dl,)  1  .5. 


mjm. 


(9.12) 

(9.13) 


where  the  indices  v  and  g  pertain  to  the  vapor  and  gas;  d  denotes  molecular 

diameter  and  d  =  0.5(d  d  ) .  Here  the  Brock  formula  acquires  the 
vg  v  gJ  n 

following  form: 


«I*  ==  ’  —3.0,807  V'2  6Kn-  J, 


(9.14) 


where  0  is  a  factor  related  to  persistance  of  the  vapor  molecules  and  is 
equal  to  0.295,  0.150  and  0.115  for  m  /m  =  0,  1  and  9,  respectively. 

As  is  clear  from  the  foregoin°,  and  as  Brock  himself  points  out,  *’  3se 

calculations,  strictly  spea'  ,  are  applicable  only  in  the  case  m  /m  H 

2 

►  0,  when  formula  (9.14)  can  be  represented  in  the  for”,  (when  -  1): 


4>  -U,l2Kn  ')  ««!>*  (1  ,  0.42Kn-’). 


(9.15) 


As  we  will  iearn  below,  this  formula  yields  the  more  or  less  correct 
value  of  ?  when  Kn  >  1 . 

We  can  find  the  lower  boundary  for  the  coefficient  for  Kn from 
elementary  considerations.  We  calculate  the  probability  that  the  evapo¬ 
rated  molecule  will  fly  off  to  distance  l  along  the  normal  to  the  particle 
surface,  there  to  collide  with  a  gas  molecule  and  return  to  the  particle 
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without  experiencing  new  collisions.  This  probability,  obviously,  is 
proportional  to  the  ratio  of  the  solid  angle,  at  which  the  particle  is 
seen  at  the  given  distance,  to  4  .  Taking  into  account  the  distribution 
of  free  paths  of  the  molecules  on  their  way  to  and  from  the  particle, 
we  find  for  the  mean  probability  of  return  of  the  molecules  the  ex¬ 
pression 


(9.16) 


When  Kn"1  is  small,  Wj^  *  0.285  Kn"1.  If  we  were  to  consider  the 

fact  that  the  molecules  evaporate  from  the  surface  not  normally,  but  at 

various  angles,  then  we  would  find  a  somewhat  higher  value  for  W1 .  The 
total  probability  of  return  of  the  molecules  to  particle  W,  is  obviously 
greater  than  W^.  Since  the  rate  of  evaporation  of  a  particle  with 

consideration  of  returns  is  $,  =  (1  -  W),  obviously  the  true  coefficient 
-1  K 

for  ;.n  should  in  any  case  be  greater  than  0.285. 

When  the  diffusion  of  vapor  molecules  is  nonisotropic  and  their 

velocities  persist  after  collisions,  the  probability  of  their  return 

to  the  particle  after  the  first  collision  is  obviously  less,  and  this 

is  increasingly  true  as  persistence  increases,  i.e.,  the  larger  the 

ratio  m  /m  .  Thus,  in  formula  (9.16)  l  should  be  defined  as  an  "effective" 
v  g 

mean  free  p-*'  of  molecules  (see  §8).  Formula  (3.15),  obviously, 

should  be  handled  in  a  like  manner  until  someone  manages  to  develop  an 
accurate  transfer  theory  that  takes  into  account  velocity  persistence. 

As  has  been  shown  [176], 

l  «  l  (m  /m  +  11/(1  ♦  S), 

eff  vv  v'  g  ’’  K  J  ’  (9.17) 


where  l is  the  mean  distance  traveled  by  vapor  molecules  between  two 


collisions;  B  is  the  number,  equal  to  0,  when  m  /m 

v  g 

monotonically  with  lAy/nig  and, reaches  the  value  1/3  when  mv/rog 
follows  from  (9.17)  and  (9.13)  that 


0,  which  increases 
».  It 


fell  =  +  1  ^0  t  ?)'«/**;•• 


(9.18) 
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After  replacing  in  the  Meyer  formula  for  ti.e  diffusion  coefficient 
of  a  vapor  in  a  gas 

Dv  —  Vvlv/3  (9.19) 

l  with  i  we  obtain  the  Dzhin  [176] 
v  eff  L  1 

D  =  V„Vl+m  yjnu  (9.20) 

v  3  (1  +  P)  Mgdlt 


which  differs  from  the  familiar  Stefan-Maxwell  formula  only  by  a  factor 
(1  +  6)  close  to  unity.  Thus,  in  formula  (9.15)  it  is  easier  to  use  the 
value  l  calculated  from  the  coefficient  of  mutual  diffusion  of  the  vapor 

and  gas  according  to  formula  (9.19). 

The  most  accurate  solution  of  the  problem  of  mass  transfer  in  aerosols 
at  intermediate  Knudsen  numbers  was  found  in  the  theory  of  neutron  trans¬ 
fer.  One  of  the  most  important  problems  of  this  theory  is  the  so-called 
"Miln's  spherical  problem",  i.e.,  calculation  of  a  flux,  to  a  black  sphere, 
of  neutrons,  diffused  isotropically  by  heavy  atoms  [177],  is  completely 
equivalent  to  the  examined  problem  for  the  case  my/mg  0* 

We  will  first  assume  that  all  vapor  molecules  possess  the  same 
absolute  velocity,  i.e.,  the  distribution  function  contains  only  a  single 
velocity  vector.  In  this  case  the  collision  integrals  were  considerably 
simplified  and  several  attempts  were  made  to  solve  the  kinetic  equation. 

The  most  complete  and  accurate  results  were  obtained  by  Sakhni  [178], 
who  based  his  work  on  the  so-called  integral  equation  and  used  in  his 
calculations  both  analytical  and  mathematical  methods. 

The  problem  at  hand  is  characterized  by  the  fact  that  a  "diffusion" 
approximation  is  often  found  instead  of  its  accurate  solution,  i.e., 
vapor  concentration  distribution  around  a  sphere,  which  satisfies  dif¬ 
fusion  equations  (9. 1-9. 3),  and  for  large  p  it  coincides  with  the  accurate 
distribution,  and  therefore  gives  the  accurate  magnitude  of  flux  to  a 
sphere  according  to  formula  (9.4).  Usually,  instead  of  the  constant  6 
in  this  formula,  the  "linearly  extrapolated  length",  i.e.,  the  expression 


is  calculated.  The  geometric  meaning  of  this  expression  is  clearly  shown 
in  Figure  4,  which  depicts  concentration  distribution  according  to  the 
accurate  solution  and  according  to  diffusion  approximation  for  Kn  *  0.5, 
by  Sakhni's  calculations.  As  seen  in  Figure  4,  at  a  distance  >  l  from 
the  surface  of  the  sphere  the  concentration  according  to  diffusion  ap¬ 
proximation  is  quite  close  to  its  accurate  value,  and  here  mass  transfer 
is  almost  exclusively  the  result  of  diffusion;  for  shorter  distances, 
however,  the  diffusion  equation  is  not  valid.  The  same  is  also  true 
of  other  Knudsen  numbers. 


Figure  k.  Concentration  Distribution 
Near  Absorbing  Sphere. 

By  comparing  formulas  (9.3)  to  (9.21),  we  find  that  6 
and  formula  (9.4)  acquires  the  form 

®-rnnGr  “me? 


■  a/(l  +  AKn) 


(9.22) 


The  values  of  A  as  calculated  by  Sakhni  are  presented  in  Table  1. 
According  to  the  author  himself,  their  accuracy  is  about  0.3%. 


TABLE  1. 


VALUES 

OF  C0EFFI 

Cl  ENT  A  ACCORDING 

Kn-‘ 

A 

Kn-‘ 

A 

v 

1,333 

1.30 

0,997 

0.1 

1 ,  CfO 

1 ,5o 

0.2 

I 

2.o  ; 

[  O.'vo 

0,5 

1.155  ! 

2,5  1 

0>'»2 

0  7  j 

1  I.H'-t  , 

3  0 

0,'l.» 

1.0  1 

I.U-tt  | 

PC  | 

0  710 

Kn  <  1  we  obtain  from  (9.22) 


<1*  =.  <\>t  (1  n,71UKn). 


(9.23) 


For  Kn  >  1  X  must  be  expressed  in  accordance  with  the  table  in  the 
form  \  -  1.333  -  0.37  Kn'1  and  D  is  replaced  by  vl/ 3.  Consequently,  we 
obtain 


<1>  -=«!»*/(!  +0,472Kn-1), 


(9.24) 


which  is  close  to  Brock's  formula  (9.15). 

Since  it  is  not  convenient  to  use  the  table  in  calculations  of 
evaporation  kinetics,  it  can  be  replaced  by  the  interpolation  formula 


1, 333  +  0.71  Kn-* 


1  + Kir* 


(9.25) 


The  maximum  deviation  of  values  calculated  according  to  this  formula  from 
those  presented  in  Table  1  is  2  to  6%.  When  the  evaporation  time  of  a 
particle  is  calculated  using  formula  (9.25)  we  obtain  an  integral  Jhat 
is  easily  taken  as  a  rational  function  of  a. 

Since  the  average  absolute  velocity  vy  molecules  does  not  figure 

in  the  derivation  of  formula  (9.24),  or  in  calculation  of  X,  the  above 
limitation- -assumption  of  equal  absolute  velocity  of  all  molecules,  is 
discarded:  the  velocity  distribution  can  be  anything  desired,  including 
Maxwellian,  and  only  the  retention  of  vy  during  collisions  and  diffusion 

isotropy,  i.e.,  the  condition  "*  0>  are  required.  Since  l  expends 

on  v  ,  the  average  should  be  used  in  (9.22). 


The  theory  of  neutron  transfer  also  takes  into  account  the  case  of 
persistence  of  their  velocities,  although  Miln's  spherical  problem  has 
no  solution  that  takes  into  account  persistence,  and  ..  r  lack  of  a  better 
approach,  persistence  should  be  taken  into  account,  as  indicated  above 
on  p.  49. 

Since  the  condensation  coefficient  a  for  certain  solids  may  be 

c 

less  than  1,  the  case  of  a  "grey"  sphere,  that  absorbs  only  the  fraction 
of  vapor  molecules  colliding  with  it,  must  be  examined.  It  is  obvious 
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that  a  grey  sphere  is  equivalent  to  one  that  possesses  some  natural 
vapor  tension  and  a£  *  1.  Through  n  and  n+  we  will  denote  the  con¬ 
centrations  of  vapor  molecules  at  the  surface  of  the  sphere,  moving  to 
and  from  it.  Obviously, 

n+  =  «_(1  —  af)  (9.26) 


The  flux  of  molecules  to  and  from  the  sphere  can  be  expressed  through 
the  formulas 


*  2 na2uyn"  -  2% a2vyn+  +  2 it a2vy(n_  -  n+),  {9.21) 

•  2ira2vyM+.  (9.28) 

However,  such  fluxes  will  also  occur  in  the  case  of  a  sphere  that  possesses 
vapor  tension  2n  ,  when  a  =  1  and  for  the  rest  cant  flux  to  the  sphere, 

equal  to 


tl>K,  ■«  2-fl  :v.  (n_  —  n+). 


(9.29) 


Thus  these  two  spheres  are  equivalent,  i.e.,  the  vapor  concentration 
distributions  around  them  are  identical.  As  follows  from  what  was  stated 
on  p.  50,  the  resultant  flux  may  also  be  described  in  the  form 

Ogt  wm  AtmD  (nx  —  2//_ )  (1  r  /J  a).  (9  J0) 

By  excluding  n+  from  formulas  (9.26)  and  (9.29)  and  setting  $gr  and 
V f  equal  we  can  find  the  expression  for  n  and,  hence,  the  formula  for  the 
flux  to  a  grey  sphere 


Or 

»  (9.31) 


which  Smirnov  derived  by  a  different  method  [179). 

We  will  return  to  the  problem  of  the  effect  of  capillary  forces  on 
the  evaporation  of  flying  particles  in  $15. 
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Proceeding  now  to  comparison  of  the  theory  outlined  here  with  the 
sparse  existing  experimental  data,  we  notice  that  direct  measurements 
of  the  rate  of  evaporation  of  HDA  particles  are  extremely  difficult  to 
make,  and  therefore  it  is  necessary  to  work  with  low  gas  pressures  (to 
increase  l) .  All  that  can  be  said  of  the  known  tests  of  Bradley  and 
coworkers  [180]  and  Monchik  and  kayse  [181]  is  that  these  tests  do  not 
contradict  the  above  theory,  although  their  accuracy  is  too  low  for 
solution  cf  the  remaining  unsolved  problem  of  the  influence  of  the 
persistence  of  vapor  molecules  on  the  mass  transfer  rate. 

Calculation  of  the  growth  or  evaporation  of  aerosol  particles  under 
actual  conditions  requires  that  the  heat  emissions  of  phase  transition, 
influence  of  particle  motion  and  other  factors  be  taken  into  account. 

Several  papers  have  been  published  in  connection  with  examination  of 
these  questions.  Brock  calculated  the  condensation  flux  to  a  particle 
moving  in  the  free  molecular  mode  [182],  and  in  another  of  his  works 
[1731,  calculated  the  liberation  of  the  heat  of  condensation  during 
free  molecular  growth.  Using  the  Belander-Bkhatnagar  model,  he  also 
studied  the  problem  of  nonstationary  particle  growth  rate  for  intermediate 
Knudsen  numbers.  Smirnov  [184]  found  expressions  for  growth  rate  with 
evaporation  rate  of  particles  with  consideration  of  the  heat  of  phase 
transition  and  a  number  of  secondary  factors  based  on  the  results  of 
neutron  transfer  theory.  Kang  Sang-wook  [18S]  used  the  boundary  sphere 
method  not  only  for  calculating  the  growth  rate  of  particles,  but  also 
for  calculating  the  rate  of  dissipation  of  the  heat  of  condensation. 

The  theoretical  results  obtained  from  his  formula  were  compared 
with  Weatherstone's  tests  on  the  condensation  of  mercury  vapors  in  a 
jet  [186]  although  this  comparison  did  not  produce  any  conclusions 
concerning  the  validity  of  the  formula  derived  in  (185]. 

Noteworthy  also  are  methods  of  calculating  the  rate  of  evaporation 
of  a  collection  of  aerosol  particles,  proposed  by  Okuyama  and  Zung  [187] 
and  rate  of  premolecular  evaporation  of  water  droplets  sprayed  into  a 
greatly  rarefied  plasma  [188]. 

In  concluding,  we  will  note  briefly  the  adsorption  ofgases  and 
vapors  on  HDA  particles.  Obviously,  the  equations  for  the  vapor  con¬ 
densing  on  particles  may  also  be  used  for  describing  adsorption  kinetics. 

The  only  difference  should  be  the  dependence  of  the  probability  of  ad¬ 
hesion  on  the  amount  of  material  sorted,  i.e.,  on  the  degree  of  coating 
of  the  surface  by  the  adsorbate.  The  over-all  process  of  the  adsorption 
of  vapors  on  acro-ol  particles  can  be  calculated  if  the  desorption 
kinetics  of  the  material  are  known. 

the  boundary  sphere  method  was  used  in  [189]  for  evaluation  of  the 
rate  of  adsorption  of  iodine  vapors  by  atmospheric  condensation  nuclei 
and  it  was  found  that  iodine  should  be  sorbed  almos'  entirely  within 
50  seconds.  In  practice,  however,  due  to  simultaneous  desorption,  rte  ( 
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fraction  of  iodine  bonded  to  the  particles  did  not  exceed  75%  of  the 
total  amount  of  iodine  in  the  atmosphere,  and  this  value  was  reached 
only  after  2  1/2  to  3  hours. 

As  regards  the  adsorption  of  the  decay  products  of  emanations,  they 
almost  entirely  ionized  [190].  Therefore  the  process  of  their  capture 
should  be  described  by  the  equations  for  the  rate  of  charging  of  HDA 
particles  by  gaseous  ions  (see  §12). 

§10.  Pulse  Transfer  in  HOA:  Resistance  of  Medium  to  Particle  Motion 

This  process  occupies  a  special  place  in  the  theory  of  transfer  in 
aerosols,  since  only  in  this  case  do  there  exist  accurate  and  completely 
reliable  experimental  data,  furnished  by  Milligan  [191],  for  oil  droplets, 
and  therefore  accurate  checking  of  the  theoretical  calculations  is  pos¬ 
sible. 

At  sufficiently  small  Reynolds  numbers  the  resistance  of  a  medium 
is  expressed,  for  Kn  -*•  0  by  the  Stokes  formula 

Ft  =  6*rflU,  (10.1) 


where  U  is  particle  velocity,  n  is  gas  viscosity. 

For  Kn  -*>  °°  Epstein  [192]  derived  the  following  formula: 

rk  (10.2) 


where  6  is  a  coefficient  whose  value  depends  on  the  mechanism  by  which 
molecules  are  reflected  from  a  particle.  In  the  case  of  total  thermal 
accommodation  and  diffusion  scattering  of  molecules  6  •  1  ♦  ir/8,  and  in 
the  case  of  mirror  reflection  0  ■  i.  By  substituting  in  (10.2)  the  known 

expression  for  gas  viscosity  n  3  0.499  n  m  v  l,  we  may  write  (10.2)  in 

S  8  & 

the  fora 


Fk  ■  (.4  -f  Q)l, 


(10.3) 


where  the  value  of  A  ♦  Q  has,  in  the  two  above-mentioned  cases,  the  values 
1.616  and  2.250.  From  (10.1)  and  (10.3)  follows 


Fe  Ft  *  (A  t  Q)  Kn. 


(10.4) 
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Mil liken  expressed  the  results  of  his  tests  with  oil  droplets  through 
the  empirical  equation 

+  *■//>]-'.  (10.5) 

in  which  A  =  1.234,  Q  =  0.414  and  b  =  0.876.  Obviously,  in  I/a  ■  Kn  09 
(10.5)  is  converted  to  (10.3).  By  comparing  the  theoretical  values  of 
A  +  Q  with  the  experimental  value  1.65,  we  may  make  the  following  formal 
conclusion:  that  a  small  part  of  the  molecules  is  reflected  from  the  -sur¬ 
face  of  the  droplets  in  the  mirror  fashion  and  that  the  rest  of  the 
molecules  experience  complete  accommodation  and  are  diffused. 

The  first  attempt  at  a  theoretical  calculation  of  the  resistance  of 
a  medium  at  intermediate  Knudsen  numbers  was  undertaken  by  Fuchs  and 
Stechkinu  [193]  using  the  boundary  sphere  method.  On  the  basis  of  a  number 
of  incompletely  substantiated  assumptions,  these  authors  derived  the 
following  formula  for  resistance  to  motion  of  oil  droplets: 

FlF*  ”  [ i  -t  oWn"  +  1 ,67Kn]  ‘  (10.6) 


The  values  of  F/F  calculated  from  this  formula  differ  by  not  more 

c 

than  2%  from  Milliken’s  data. 

By  breaking  down  the  collision  integral  in  the  kinematic  equation  in 
terms  of  Hermit^  polynomials  and  using  the  Knudsen  iteration,  Liu  and 
coworkers  [194]  derived  the  formula 


F:F„ .«  1 0,298Kn->. 


(10.7) 


By  linearizing  the  kinetic  equation  as  indicated  in  §9,  Willis  [195] 
using  the  Knudsen  iteration,  derived  the  expression 


FFk.*  '1  -  0.36GKn-».  (10.8) 

Chercignani  and  Pagani  [196]  solved  transfer  problems  in  the  range 
of  intermediate  Knudsen  numbers  by  a  new  method.  After  the  linearization 
of  the  kinetic  equation  the  authors  constructed  a  functional  from  the 
functions  that  express  the  density  of  the  gas,  velocity  of  the  molecules 
and  temperature,  and  used  the  variation  method  to  find  the  stationary 
value  of  the  function.  Here,  the  above-specified  physical  values  acquire 
values  that  correspond  to  the  solutions  of  the  linearized  equation,  but 
the  fluxes  of  pulse,  heat,  etc.  arc  calculated  directly  from  the  stationary 
values  of  the  functional.  Using  this  method  the  authors  were  able  to 
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solve  an  entire  series  of  transfer  problems  in  the  systems  with  different 
geometry. 

Recently,  Cherciyani  and  Fagani  [197]  published  a  work  in  which  the 

above  method  was  used  for  the  problem  of  the  resistance  of  a  medium  to 

the  motion  of  spheroidal  particles.  Their  calculation  results  for 

spherical  particles  are  depicted  in  Figure  5  as  the  dependence  of  F/F.  on 

-1  K 
Kn  .  The  calculations  quite  satisfactorily  approximate  the  experimental 

formula  (10.5);  the  slight  difference  between  them,  which  never  exceeds 

2%,  can  be  attributed  to  the  fact  that  in  the  expression  for  F^  the  authors 

used  for  the  coefficient  8  the  value  1  +  tt/8,  i.e.,  they  assumed  that  all 
gas  molecules  experience  tdtal  accommodation  on  the  surface  of  particles, 
whereas  in  actuality  a  small  portion  of  them  is  reflected  in  mirror  fashion. 
The  works  of  Cherciyani  and  Pagani  doubtless  have  made  an  extremely  im¬ 
portant  contribution  to  the  theory  of  transfer  phenomena  at  intermediate 
Xnudsen  numbers. 


i 


Figure  5«  Kinetics  of  Transfer 
Processes  at  Intermediate 
Knudsen  Numbers. 


It  is  clear  in  Figure  5  that  Willises'  formula  (10.8).  as  expected 
gives  the  correct  values  of  F/F^  only  approximately  to  Kn"1  m  0.5,  while 

formula  (10.7)  is  completely  inaccurate.  As  regards  the  "universal" 
Sherman  formula  (see  equation  (9.11)), 
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(10.9) 


F  t\  ---  ,  l  v-. 

*  I  rr  n.  r  ( 

which  can  be  represented  through  formula  (10.4)  as 


P  f  1  1 

*  i  -  kir*  (.t  4*1"*  | -r  n.itO<»Kn ' 


which  then  yields  somewhat  understated  (up  to  7%)  values  of  F/F^ .  It  is 

extremely  fortunate  that  the  values  of  calculatet^  according  to 

Saklmi's  data  for  mass  transfer  (from  Table  1)  with  the  aid  of  formula 

(9.7),  as  functions  of  Kn ' 1 ,  quite  satisfactorily  rely  on  the  curve  for 
pulse  transfer,  as  seen  in  Figure  5.  At  the  same  time  the  curve  for  heat 
transfer  lies  somewlut  higher. 

Brock  [198]  derived  quite  a  complex  expression  for  the  resistance  of 
a  medium  to  the  motion  of  growing  or  evaporating  particles  at  high  Knudsen 
numbers,  which  undoubtedly  will  be  of  benefit  for  calculation  of  the 
formation  of  aerosols  in  rocket  nozzles  and  influence  of  the  particles 
formed  on  the  flux  properties. 

511.  Heat  Transfer  in  HOA:  Heating  and  Cooling  of  Particles 
0 

The  importance  of  the  process  of  heat  transfer  to  particles  from  its 
surroundings  lies  primarily  in  the  fact  that  the  temperature  of  an 
evaporating  or  growing,  due  to  condensation,  particle  diffors  from  the 
temperature  of  the  surrounding  environment,  and  this  fact,  in  the  case  of 
materials  with  notable  vapor  tension,  has  a  considerable  effect  on  the 
kinetics  of  evaporation  and  growth  of  particles. 

For  Kn  »  0  the  rate  of  heat  transfer  to  a  particle  is  expressed  by 
a  formula,  derived  in  classical  thermal  conductance  theory,  analogous  to 
(9.5), 


Qe  ■=»  4sa/A7\ 


(11.1) 


where  \  is  the  heat  transfer  coefficient  of  the  medium;  AT  is  the  tempera¬ 
ture  difference  between  the  medium  and  the  particle.  We  will  confine  our 
examination  to  the  simplest  case  where  AT  is  small  compared  to  the  absolute 
particle  temperature  T. 

To  derive  the  expression  for  heat  flux  to  a  particle  when  Kn  •*  •  we 
will  introduce,  as  in  the  derivation  of  formula  (9,31),  the  indices  ♦  and 
-  to  denote  molecules  moving  away  from  and  toward  a  particle.  From  the 
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or 


condition  of  material  balance  it  follows  that  n+  v  +  =  n"  v 

g  g  g  g 


(11.2) 


The  flux  of  molecules  to  a  particle  is  (insert  a) 
and  the  mean  thermal  energy  of  one  molecule  is  c^T" ,  where  *  3/2k  for 

monatomic  and  5/2k  for  diatomic  gases.  However,  the  energy  flux  to  a 
particle  cannot  be  found  by  multiplying  these  two  values,  since  the  greater 
the  kinetic  energy  of  molecules,  the  relatively  more  often  they  collide 
with  the  particles.  The  correct  expression  for  energy  flux  [  1 99]  is  of 
the  form 


4 


/  - 


(11.3) 


The  analogous  expression  (in  which  the  -  is  replaced  by  the  +)  is  also 
found  for  the  energy  flux  carried  away  by  the  molecules,  and  the  resultant 
flux  is  equal  to  the  difference  of  these  expressions: 


Q<*  - 


i  T  )  !  •'  —  «*  (f  r 


(11.4) 


We  will  introduce  the  coefficient  of  temperature  accommodation  of 
molecules  a.j.  =  A'T/AT,  where  A'T  is  the  difference  between  the  average 

temperature  of  molecules  departing_ from+particles  and  the  temperature  of 
the  environment;  obviously  AT'  »  T"  *  T  .  Assuming  A'T/T  to  be  a  small 

value,  we  find  from  (11.4),  using  condition  (11.2),  the  formula 


Q*  “■  ~aHlu^)  {cp  ~  l)nsT'  '  A  T  ^-a:n^/e{cvn-^)ir^T.  (11.5) 


We  will  now  find  the  value  of  the  ratio  Q^/Q^ .  From  the  known  formula 
for  the  thermal  conductance  of  gases  x  =  u ncv/m^ ,  where  n  is  the  viscosity 

of  the  gas  and  the  coefficient  c  is  equal  to  2.5  for  monatomic  and  1.9 

for  diatomic  gases,  and  from  the  formula  n  *  O.Sn  m  v  l,  we  obtain 

g  g  g 
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X  =  0,5  tc0ntvgl. 


(11.6) 


hence 

0,/Q*  =  2icvKn  I  [cv  +  |)  ar.  ( n 

Thus,  for  monatomic  gases  Qc/Qk  =  3.75  Kna,^ ,  and  for  diatomic  gases 
Qc/Qk  =3.17  Kna'1. 

For  an  intermediate  Knudsen  number  Springer  and  Psai  [200]  derived, 
using  the  boundary  sphere  method,  the  formula 


(11.8) 


where  B  -  3.75  for  monatomic  and  3.17  for  diatomic  gases.  When  Kn  -*■  « 
this  formula  acquires  the  form  Q/Qc  =  (B  Kna^1)"1.  Considering  the 

expressions  for  Qc/Qk»  we  see  that  Q/Qk  -*■  1,  in  this  case  as  should  be 
expected. 

Using  the  moment  method,  Lees  [201]  found  that  for  diatomic  gases 
the  formula 


Q  Qc  ~  (1  -f  2,5  Kn)'1. 


(11.9) 


which  is  clearly  unsuitable  for  large  Knudsen  numbers. 

Sherman's  universal  formula  has  in  this  case  the  form 

^  :r=  l  -  Qk-Qc  I  +  parKn-‘  (11.10) 

[see  formula  (11.7)  concerning  the  value  of  B] . 

The  above-described  Knudsen  iteration  method  was  used  by  Brock  [202] 
for  derivation  of  the  formula 

q  Qk  l  —  ar5Kn-’,  (11.11) 
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where  B  =  0.153  for  monatomic  and  0.239  for  diatomic  gases.  Formulas  of 
this  type,  as  we  have  already  stated,  are  valid  only  for  large  Knudsen 
numbers . 

Cherciyani  and  Pagani  [203]  used  the  variation  method  mentioned  in 
the  preceding  section  (assuming  a,j,  =  1)  to  calculate  the  values  Q/Q^  for 

monatomic  gases  in  the  range  of  Knudsen  numbers  from  0.0001  to  50. 

Figure  4  shows  the  curve  plotted  according  to  these  data. 

Up  until  just  recently  the  only  experimental  determination  of  the 
rate  of  heat  transfer  to  a  .sphere  at  low  Reynolds  numbers  and  at  inter¬ 
mediate  Knudsen  numbers  was  made  by  Takao  [204]  who  measured  the  kinetics 
of  cooling  of  the  heat  glass  sphere  1.5  cm  in  diameter,  placed  in  a  vessel 
measuring  about  1  m,  from  which  the  air  was  exhausted  to  the  desired 
residual  pressure.  All  the  above-cited  authors  compared  the  formulas 
which  they  derived  (Takao  himself  also  derived  an  extremely  complex  formula 
with  Takao' s  data  and  they  all  found  excellent  coincidence.  It  should  be 
pointed  out  that  Takao  expressed  his  results  in  the  form  of  the  depend¬ 
ence  of  Q/Q  on  Knudsen  number.  Generally  speaking,  however,  Takao 's 
c 

data  are  sufficiently  accurate  only  for  small  Knudsen  numbers  and  are  er¬ 
roneous  for  Kn  >  1,  as  indicated  by  the  strong  scattering  of  experimental 
points  in  this  range.  For  Kn  >  1,  apparently,  recently  published  results 
of  Harbor  [205],  also  in  error,  are  more  accurate. 

Figure  5  shows,  in  addition  to  Cherciyani 's  curve,  Brocks  curve 
(11.11),  which  nearly  coincides  with  it  at  large  Knudsen  numbers,  and 
the  experimental  data  of  Harbor  and  Takao.  Unfortunately,  the  values 
of  in  the  test  of  Takao  and  Harbor  are  unknown,  which  makes  it  dif¬ 
ficult  to  compare  the  theory  with  the  experiment. 


As  already  mentioned  in  §8,  it  is  necessary  in  heat  transfer  cal¬ 
culations  to  use  the  mean  free  path  of  molecules  specifically  for  this 
process,  i.e.,  I rather  than  the  generally  accepted  value  for  moment 

transfer  l .  According  to  the  calculations  of  Lord  and  Harbor  [165]  for 

diatomic  molecules,  *  1.6.  If,  however,  Chercyani's  curve  is 

related  to  rather  than  to  it  must  be  displaced  to  the  left  in  such 

a  way  as  to  reduce  the  abscissa  by  a  factor  of  1.6.  Then  this  curve  will 
nearly  coincide  with  the  analogous  curves  for  mass  transfer  and  pulse 
transfer  and  we  will  have  practically  the  same  curve  for  all  three  processes. 
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612.  Charge  Transfer  in  HDA:  Charging  and  Discharging  of  Particles 
by  Gas  Ions 

One  of  the  most  important  physical  processes  in  aerosols  is  charge 
transfer  to  particles  by  ions  captured  by  them.  It  is  noteworthy  that 
until  just  recently  we  did  not  know  the  exact  nature  of  a  gas  ion.  The 
whole  series  of  facts  indicates  that  primary  ions,  formed  when  an  electron 
separates  from  or  combines  with  molecules  of  a  gas,  are  united  by  the 
electrical  forces  of  neutral  molecules  that  possess. constant  dipole 
moment,  forming  aggregates.  The  electrical  mobilities  of  atmospheric 
"bright"  ions  lie  in  the  range  from  about  0.5  to  2.5  cm2*B_1 ‘sec-1 . 

Hence  the  number  of  molecules  in  ionic  aggregates  also  varies  substantially. 

Since  the  expressions  for  the  electrical  mobility  and  coefficient 
of  diffusion  of  ions  [199] 


B 


mg  +  mi  \»<» 
mtmikT  )  ’ 


(12.1) 

(12.2) 


contain  both  mean  free  path  of  ions  l.  and  their  mass  m.  (m  is  the  mass 
of  neutral  gas  molecules) ,  we  can  determine  l i  only  by  assigning  a  certain 
value  to  iik  .  If  we  assume  itk  =  m^  and  take  for  mobility  the  mean  value 
B  =  1.4  cm2*B-1*sec”1,  then  in  air  of  20°  and  at  normal  pressure  we  will 
have  l ^  =  1.3*10'6cm,  and  for  the  mean  thermal  velocity  v^  =  4.62*104  cm/sec. 


If,  however,  like  Young  [206]  we  assume  that  itu/m^  =  6-18,  we  find 

for  l.  values  of  about  1.8*10-6cm  and  for  v. ,  1.0-2.0*104  cm/sec.  This 
i  i 

uncertainty  in  the  values  of  l ^  and  v^  greatly  impede  comparison  with 

experimental  data  in  charge  transfer  theory  in  aerosols.  We  will  note 
also  that  the  mode  of  charge  transfer  to  particles  is  determined  in  this 
case  by  the  "Knudsen  ion  number"' Kn^,  the  value  of  which  is  several  times 

less  than  the  "molecular"  number  Kn,  so  that  Kn^  =  0.13-0.18  corresponds 
to  the  upper  boundary  of  HDA  particle  radii  which  we  have  selected,  namely 
10~5  cm. 


We  will  first  examine  the  charging  of  particles  by  ions  of  the  same 
sign  in  the  absence  of  an  external  electrical  field.  The  differential 
equation  of  stationary  diffusion  of  ions  to  a  particle  in  electrical  field 
created  by  it  has  a  form 


-56- 


i  12.  3) 


where  E  is  the  voltage  of  the  field.  When  Krr  *  0  we  have  the  boundary 

conditions  n  =  n*  for  a  ♦  *•  and  n  =  0  for  .  =  a  and  we  arrive  at  the 
formula 


where 


(12.4) 


['o' 

:>P»  (p*-o‘)  ’ 


(12.5) 


--the  potential  of  electrical  forces  acting  on  an  ion,  e  is  the  elemental 
charge;  i  is  the  number  of  elemental  charges  on  a  particle,  positive  for 
like  ionic  and  particulate  charges  and  negative  for  unlike.  The  first 
term  in  (12.5)  corresponds  to  the  Coulomb  force  and  the  second  to  mirror 
force.  We  assume  here  that  the  particle  is  a  conductor  of  electricity,  j 

which  is  valid  in  most  cases  for  extremely  fine  particles.  It  is  assumed  ? 

in  the  derivation  of  (12.4)  that  the  process  of  diffusion  of  ions  is 
quasistationary.  Under  what  conditions  this  is  true  will  be  explained 
below. 

To  answer  the  question  whether  mirror  forces  can  be  disregarded  in 
the  theory  of  charging  of  HDA  particles  we  will  calculate  their  effect 
during  diffusion  of  ions  to  an  uncharged  particle  (i  =  0).  The  magnitude 
of  the  integral  in  (12.4)  in  this  case  is  0.97,  0.83  and  0.62  for  a  = 

=  10"4,  10" 5 ,  and  10~6  cm,  respectively.  Thus,  mirror  forces  in  HDA  can¬ 
not  be  disregarded. 

Proceeding  to  the  case  Kn^  -»  »,  we  will  note  that  the  mechanism  of 

charging  of  particles  by  gaseous  ions  consists  in  the  fact  that  ions  in 

collision  with  particles  give  up  their  charge  and  fly  off  in  the  form  of 

neutral  molecules.  Gentry  and  Brock  [207]  instead  of  disregarding  mirror 

forces,  make  an  incorrect  assumption  with  Boltzmann  distribution  of  ion 

concentrations  around  the  particles.  Since  this  distribution  is  governed 

by  molecular  collisions,  and  since  the  latter  do  not  exist  near  a  particle 

for  Kn.  -+<»,  the  formula 
i 
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4»,  --  ra'r.rt* exp  (—U‘  akT) 


(12.6) 


derived  in  the  above-cited  work  is  incorrect.  Liu,  Whitby  and  Yu  (208} 
made  the  same  mistake. 

The  correct  solution  of  the  problem  is  found  in  the  works  of  Keefe, 
Nolen  et  al.,  [209-211],  who  solved  the  problem  of  the  free  trajectory  of 
an  ion  travelling  in  an  electrical  field  created  by  a  particle.  By  dis¬ 
regarding  mirror  forces,  these  authors  found  for  the  ionic  flux  to  a 
particle  of  opposite  charge,  in  consideration  of  Maxwellian  velocity 
distribution  of  the  ions,  the  expression 


‘I\.  -ii-  ak'l T, 


(12.7) 


which  coincides  with  (12.6)  only  when  ie2/akT  <  1.  When  mirror  forces 
are  taken  into  account  the  ionic  flux  to  an  uncharged  particle  is 


•I**--  r.Q:v,nx(\  —  ]  r.t  'lakT). 


(12.8) 


In  the  case  of  charged  particles  the  problem  can  be  solved  only  by  the 

mathematical  method  [211].  The  values  of  the  "electrostatic"  factor  6 

(inclosed  in  braces  in  (12.7)  and  (12.8))  are  presented  for  this  case  in 

Figure  6  as  functions  of  particle  radius.  The  curves  in  this  figure  were 

plotted  on  the  basis  of  the  mean  square  velocities  of  the  ions.  As  Keefe 

et  al.,  show,  the  curves  consider  the  velocity  distribution  of  the  ions 

and  can  he  satisfactorily  approximated  by  displacing  the  curves  shown  in 

the  figure  to  the  right  so  that  the  abscissas  are  1.5  times  greater.  The 

authors  assume  in  their  calculations  that  ionic  aggregates  consist  of 

eleven  water  molecules,  i.e.,  they  assumed  m./m  *=  7. 

i  g 

The  problem  of  the  charging  of  particles  at  intermediate  Knudsen 
numbers  Km  is  fraught  with  enormous  difficulties  and  remains  far  from 

solved  even  today.  Proceeding  from  formula  (12.6),  Gentry  and  Brock  [207] 
found  by  the  Knudsen  iteration  method  an  expression  for  in  a  form 

with  which  we  are  already  familiar.  Lie,  Whitby  rnd  Yu  went  even  farther 
[208]:  they  assumed  not  only  the  Boltzmann  distribution  of  ion  concentra¬ 
tions  near  a  particle,  but  also  the  Maxwellian  (isotropic)  distribution 
of  their  velocities  and  arrived  at  the  con-lusion  that  (12.6)  is  valid 
for  all  Knudsen  numbers! 
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0,6  1 


i 


Figure  6.  Electrostatic  Factor 
During  Diffusion  of  Ion  to 
Charged  Particle. 


All  other  authors  solved  this  problem  by  the  boundary  sphere  method. 
The  diffusive  flux  to  it  was  determined  on  the  basis  of  formula  (12.4) 
in  which  a  is  substituted  by  the  radio  of  the  boundary  sphere  a  +  SZ^. 

The  free  molecular  flux  to  a  particle  from  the  boundary  sphere  was  cal¬ 
culated  by  various  methods.  Brickard  [217]  and  Siksna  [213]  did  not  take 
into  account  the  curvature  of  the  trajectory  of  the  ions  in  the  electrical 
field  cf  the  particle.  Keefe  and  coworkers  [211]  used  in  their  calcula¬ 
tions  the  results  outlined  above,  which  they  obtained  for  Kn.  -*•  “,  i.e., 

they  assumed  that  the  ions  fly  toward  a  particle  as  though  from  infinity. 

In  his  calculations  Fuchs  [214]  assumed  that  the  initial  point  of  the 

trajectories  of  the  ions  lies  on  the  boundary  sphere.  Natanscn  [215] 

selected  as  his  boundary  sphere,  a  sphere  such  that  an  ion  experiencing 

a  collision  within  it  cannot  be  ejected  from  it.  All  above  cited  authors 

based  their  calculations  on  extremely  varied  values  of  Z.  and  v.  and  it 

ii 

is  therefore  difficult  to  compare  their  results.  Since  the  boundary 
sphere  method  generally  does  not  produce  accurate  ■esul1'",  it  is  hardly 
worthwhile  to  consider  the  question  as  to  which  of  the  above-mentioned 
methods  is  more  accurate. 

In  all  the  foregoing  we  assumed  the  process  of  charging  and  dis¬ 
charging  of  particles  to  be  quasi-stationary.  Meanwhile,  during  each 
event  of  recharging  of  particles  the  distribution  of  ion  concentrations 
around  them  must  change.  It  is  obvious  that  in  order  that  the  discharging 
process  be  considered  quasi-stationary  it  is  essential  that  the  time  of 
relaxation  of  the  ionic  atmosphere  around  the  particles  be  much  less  than 
the  average  time  that  elapses  between  recharging  events.  As  analysis  of 
the  problem  indicates  [214],  the  condition  a^n^  <  1  when  a  >  10*5  cm  must 


-59- 


be  satisfied  for  this  purpose,  When  a  »  10*f  c*  the  more  rigid  condition 

a:nt  <  0.01  must  be  satisfied.  In  practice,  these  conditions  may  be  re¬ 
garded  as  satisfied,  both  in  the  laboratory  and  in  the  atmosphere. 

Under  the  above  stated  conditions,  as  indicated  by  calculations, 
the  electrical  field  created  by  the  nonuniform  distribution  of  ion 
concentration  near  a  charged  particle,  i.8.,  the  screening  effect,  may 
also  be  disregarded. 

In  the  case  where  particles  become  charged  in  a  bipolar-ionized 
environment,  it  can  be  assumed  that  ions  of  both  charges  diffuse  toward 
the  particle  independently  of  each  other.  A  stationary  distribution, 
ia  which  the  number  of  particles  with  i  charges  of  a  given  sign,  which 
acquire  yet  another  charge  of  the  same  sign  per  unit  of  time,  is  equal 
to  the  number  of  particles  with  i  +  1  charges,  acquiring  one  charge  of 
the  opposite  sign,  will  eventually  be  established  in  an  aerosol,  regard¬ 
less  of  the  initial  distribution  of  the  particles  with  respect  to  charges. 
For  simplicity  we  will  assume  that  an  aerosol  is  monodispersed  and  that 
the  concentration,  mobility  and  other  properties  of  ions  of  both  signs 

coincide.  By  determining,  on  the  ba.iis  of  the  above  formulas,  the  fluxes 
of  like  and  unlike  ions  to  particles  with  i  charges  ^  and  $!,  we  can 

formulate  a  series  of  equations 


2Vo  ■  nNi; 


4  N 
11 


♦  .K. 

i  t 


®v*-lNi+l’ 


(12.9) 


when  N.  is  stationary  concentration  of  particles  with  i  charges,  and  hence 

we  can  find  the  distribution  of  particles  in  terms  of  charge  (the  coef- 
ficient  2  in  the  first  equation  is  a  result  of  the  fact  that  an  uncharged 
particle  can  capture  an  ion  of  either  sign).  We  will  note  that  this  dis¬ 
tribution  is  independent  of  ion  concentration,  al 'hough  the  rate  of  its 
establishment  is  proportional  to  n°°. 

For  large  particles  (a  >  10'4  ■  ),  when  formula  (12.4)  is  applicable 
and  mirror  forces  can  be  disregarded,  this  formula  acquires  one  of  the 
following  forms: 


<|*  -\r.Dnn  ;  4>, 
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Dunjf  j 
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(12.10) 


where  X^  = | ic2/akl| . 


It  follows  from  (12.10)  that 


expt-A^ . 
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By  aultinlvins  equation  (12.9)  we  obtain 


*  A'o  exp  (—  *7i:  2akT)  (exp  Ui/2)  -  exp  (—  /.,/2)|V 


(12.11) 


As  shown  by  calculations  the  factor  [exp(X../2)--exp(-X^/2]/x^  for 

particles  with  a  >  10'4  cm  and  with  not  very  large  charges  fthe  over¬ 
whelming  majority  of  particles  fits  this  description)  is  quite  close  to 
1.  If  we  discard  it  we  obtain 


Nt  -  N0  exp  (-i2t-l2ab.T), 


(12.12) 


i.e.}  particle  distribution  with  respect  to  charges  is  expressed  by  the 
Boltzmann  eqir.cion. 

We  will  note  that  formula  (12.12)  is  extremely  accurate  when  a  >  10‘4 

cm,  but  as  particle  dimension  decreases  it  becomes  increasingly  inaccurate. 
In  connection  with  the  method  described  in  §4  for  determining  particle 

dimension  according  to  the  fraction  of  the  charged  particles  in  the 

stationary  state,  we  point  out  that  for  a  =  10'5  cm  this  fraction, 
calculated  in  accordance  with  the  Boltzmann  formula,  is  ~2%  (relative) 

less  than  the  more  accurately  calculated  value,  ~10%  for  3*10~6  cm  and 

about  double  that  for  10'6  cm. 

Keefe  et  al.,  [212]  used  the  Boltzmann  particle  distribution  by 
charge  a  priori,  and  various  explanation?  are  given  for  departures  from 
it  for  small  a.  Meanwhile,  the  Boltzmann  equation  is  valid  only  for 
equilibrium  states;  if  ions,  having  combined  with  particles,  could  break 
away  from  them  in  the  form  of  ions,  i.e.,  if  there  existed  "detailed 
equilibrium",  the  Boltzmavn  formula  unquestionably  would  be  applicable. 

In  reality,  however,  we  don't  have  equilibrium  here,  but  rather  a 
stationary  state. 

As  we  already  pointed  out  in  the  introduction,  the  number  of  elemental 
charges  acquired  by  HDA  particles  in  a  bipolar  ionic  atmosphere  is  very 
small.  By  way  of  illustration  we  present  a  table,  compiled  by  Fuchs  [214], 
of  HDA  particle  distribution  by  cl. irge  in  the  stationary  state  (rounded 
off  to  the  nearest  hundred) . 
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TABLE  2.  STATIONARY  01  ST*.ltUT ION  OF  HOA 
PARTICLES  BY  CHARi'.C  IN  BIPOLAR 
IONIZED  ATMOSPHERE. 


j  (cm) 

!0“7 

r-~ - 

3-I0"7 

I0“6 

3* 10“6 

10“5 

Fraction  of  uncharged 
particles 

0.99 

0.95 

0.76 

0.51 

0.29 

Fraction  of  particles 
with  1  c 

l 

0.01 

i 

0.^5 

0.24 

0.45 

0.44 

2  £ 

- 

- 

- 

0.04 

0.20 

3  c 

- 

- 

- 

- 

0.06 

4  e 

- 

- 

- 

- 

0.01 

1 _ 

In  the  presence  of  an  external  electrical  field  the  charge  is 
superposed  on  the  diffusion  charge  of  the  particles  by  the  flux  of  ions 
moving  in  the  field.  In  the  case  of  unipolar  ionization  this  flux  is 
expressed  by  the  Poten’ye  formula 


$/:  -  3i:£na8fla1<l~H/3£aI)\ 


(12.13) 


where  E  is  the  field  voltage.  For  comparison  of  this  flux  with  diffusion 
flux  we  will  examine  the  charge  of  a  neutral  particle  with  a  *  10" 5  cm. 

In  this  case  we  can  express  diffusion  flux  without  large  error  by  formula 
(12.10)  for  <t  .  This  formula  and  formula  (12.13)  it  is  easy  to  find  that 


<I'c  /  «I>0  -  0.75i£a  kT  *  0.9-10 *Ea 


(12.14) 


or  for  a  *  10  ^  cm,  expressing  E  in  v/^-a, 


< t»a  •--=  310*4£. 


(12.15) 


This  ratio  is  even  smaller  for  particles  of  smaller  dimension  and  for 
charged  particles.  Thus,  the  rate  of  "drift”  charging  of  HDA  is  compared 
with  the  d;ffusion  rate  only  in  sufficiently  strong  fields  (in  excess  of 
3,000  v/cm),  for  instance,  in  clectro-fi Iters ,  and  it  can  be  disregarded 
in  the  atmosphere. 

Proceeding  to  experimental  data  on  the  kinetics  of  the  charging  of 
aerosol  particles  by  gas  ions,  we  will  point  out  that  most  of  the  measure¬ 
ments  were  made  on  comparatively  large  particles  at  atmospheric  pressure, 
i . c . ,  at  small  Kn^  numbers.  Flanagan  [216]  measured  the  reduction  of  the 

stationary  concentration  of  gas  ions  in  air  u.  Ter  the  influence  of  the 
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. polydxspersed  aerosol,  formed  by  the  heating  of  Nicbrome  with  a  from 
0.01  to  0.13  m,  i.e.  for  Kn^  >  1.  Considering  that  the  particles  thereby 

did  not  acquire  more  than  one  elemental  charge  (which  is  true  only  for  th*. 
lower  end  ot  this  range),  Flanagan  was  able  from  these  measurements  to 
determine  the  sum  of  values  and  : ^ ,  specifically  :Q  +  : ^ ,  *  0.58  a. 

For  Kn.  >  1  this  dependence  correlates  thoroughly  with  the  theory. 

Liu,  Wiitney  ;*nd  Yu  [216]  measured  the  kinetics  of  the  charging  of 
monodispersed  aerosols  of  dioctylphthaiate  with  a  =  0.32  and  0.67  u  by 
unipolar  ions  at  pressures  of  0.03  to  1  atm,  i.e.,  Kn.  <  1  and  found 

that  the  curve  of  particle  charging  as  a  function  of  time  agrees  satis¬ 
factorily  with  formula  (12.6)  which  docs  not  contain  Kn^,  and  tbe  authors 

assumed  this  formula  to  be  valid  for  all  Kn^.  For  this  purpose,  of  course, 

they  had  to  assume  that  the  mass  of  «!»c  ions  was  15  times  greater  than 
the  mass  of  the  neutral  ir  molecules.  As  already  pointed  out,  formula 
(12.6)  is  invalid  at  such  Km  numbers,  and  it  is  impossible  to  explain 

how  such  agreement  was  achieved.  Otherwise,  it  would  have  to  be  said 
that  the  charge  tests  conducted  for  low  pressure  cannot  accurately  simulate 
the  process  of  charging  of  HDA  particles  due  to  the  effect  of  mirror  forces, 
the  magnitude  of  which  depends  not  on  Kni  but  on  particle  dimension. 

A  number  of  researchers  measured  the  ratio  of  charged  and  uncharged 
aerosol  particles  in  the  stationary  state  in  a  bipolar  ionic  atmosphere 
from  which  was  calculated  the  relation 


-  2N>iNx 


These  studies,  were  made  with  poiydispersed  aerosols  and  can  hardly 
be  used  for  checK.ng  the  theoretical  calculations.  Fortunately,  in  the 
tests  of  Poliak  and  Metniks  (217]  with  aerosols  obtained  by  the  heating 
of  Nichrooc,  for  which  a  *  (1  -  4)*10'1  cm,  it  was  found  that  the  experi¬ 
mental  values  of  the  ratio  N^/N^,  in  accordance  with  theory,  are  higher 

than  those  calculated  according  to  the  Boltzmann  formula  for  a  <  1.35* 

•  10’^  cm,  but  lower  for  larger  particles.  As  pointed  out  b..  the  above- 
cited  authors,  the  reason  for  this  is  the  polydispersion  of  the  aerosols, 
and  the  larger  the  particles  the  more  it  distorts  the  results. 
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:>I3.  Thermophoresis  and  D i  f f us iophores i s  In  HDA 

In  our  discussion  of  these  complex  phenomena  we  will  confine  ourselves 
only  to  the  range  Kn  >  1.  * 

A  particle  located  in  heterogeneous  gaseous  medium  is  acted  upon  by 
a  force,  known  as  thermophoretic  if  it  is  caused  by  the  presence  of  a 
temperature  gradient  in  the  mixture,  and  diffusiophoretic  in  the  case  of 
a  composition  gradient.  In  these  cases,  as  in  the  case  of  the  problem 
concerning  the  resistance  of  a  medium  to  particle  motion,  the  force  acting 
on  area  S  of  a  particle  is  defined  by  the  integral 

F».—  5,iS2S5«,|»,V,l/<«r 

*4» 

where  f+  and  f  are  the  distribution  functions  of  gas  molecules,  the  pro¬ 
jections  of  whose  velocities  onto  the  direction  of  the  concentration 
gradient  are  positive  or  negative,  respectively,  is  the  velocity  of 

gas  molecules.  The  temperature  gradient  or  concentration  gradient  appears 
in  this+equation  in  implicit  form  through  the  coordinate-dependent  func¬ 
tions  f  and  f* .  The  causes  of  motion  of  aerosol  particles  under  influence 
of  temperature  gradient  differ  for  Kn  >  1  and  Kn  <  1.  When  Kn  >  1  this 
phenomenon  is  quite  analogous  to  thermodiffusion.  In  the  case  of  Kn  <  1 
the  cause  of  their  motion  is  the  existence  of  a  temperature  gradient  within 
the  particle  itself  and  the  tangential  flow  of  gas  along  the  unevenly  heated 
surface  of  the  particle. 

In  contrast  to  other  transfer  phenomena,  even  in  the  Kn  >  1  range,  the 
zero  approximation  for  f  cannot  be  used  in  the  calculation  of  thermo¬ 
phoretic  force,  since  the  gas  is  in  a  heterogeneous  state. 

The  theory  of  gaseous  thermodiffusion  was  used  for  solving  the  problem 
of  thermophoresis  in  tne  free  molecular  regime  by  Einstein  [219],  Clusuis 
[221],  and  in  a  more  rigorous  form,  by  Bakanov  and  Deryagin  [223],  and  also 
by  Waldmann,  [222].  In  the  two  latter  works  the  Chapman-Enskog  theory  was 
used  in  the  Lorentz  problem,  i.e.,  it  was  assumed  that  the  mass  of  the 
aerosol  particles  is  great  compared  to  the  mass  of  gas  molecules  and  that 
the  role  of  the  collisions  of  molecules  with  eath  other  is  negligibly  small 
compared  to  their  collisions  with  the  surface  of  particles.  Assuming  that 
part  of  the  molecules  1  -  <J>  is  reflected  in  mirror  fashion  and  the  ethers 
arc  diffused,  i.e.,  with  Maxwellian  distribution,  corresponding  to  the 
temperature  of  the  surface  of  a  particle,  the  thermophoretic  force  will 
not  depend  on  $  and  is  expressed  by  the  formula 

Ft*  =  —  (32,15)  (a-'y.jVr/o,,  (13.1) 
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where  AT  is  the  temperature  gradient. 


x  is 


the  thermal  conductivity  of 


the  gas.  In  the  case  of  a  multi-atomic  gas  the  value  of  the  translation 
component  of  thermal  conductivity,  Xg^  =  15n/4m^,  where  n  is  viscosity, 

should  be  substituted  into  (12.1).  Mason  and  Chapman  [224]  assumed  that 
all  molecules  are  reflected  elastically,  some  of  them  in  mirror  fashion 
and  portion  <j>*  in  such  a  way  as  to  retain  absolute  velocity,  but  in  random 
direction.  They  obtained  a  value  of  1  +  4n/9  greater  than  according  . 


(13.1).  Mohchik,  Inn  and  Mason  [227],  using  a  model  of  inealstic  col¬ 
lisions  developed  by  Wang  Chang  and  Uhlenbeck  [226],  strictly  substantiated 
the  validity  of  formula  (13..1)  for  multi-atomic  gases. 


The  theory  of  thermophoresis  in  the  range  of  intermediate  Knudsen 
numbers  was  developed  by  Brock  [227]  by  the  Bkhatnagar  method  using 
Knudsen  iteration.  The  thermophoretic  force  for  a  single  component 
multi-atomic  gas,  according  to  Brock,  is 


Ft  —  FTi;  exp  {1  —  {•),06+  0,09a,l+,i.28(I  —  2/JKn'. 


(13.2) 


where  F_.  is  given  by  equation  (13.1),  x„  and  x„  are  the  thermal  conductiv-  . 


g 


ity  of  gas  and  particle,  a  is  the  coefficient  of  mechanical  accommodation. 

m 

This  equation  is  equivalent  to  the  emp  leal  formula  derived  by  Schmitt 
[228] 


Ft  =  Frk  exp  (~  b  Kn), 


(13.3) 


assuming  that  for  a  monatomic  gas 


b  —  0,1)6+0,37  am— 0,09am/a/2/r , 


(13.4) 


and  for  a  multi-atomic  gas 


b 


-  0,06+0,37am  —  0,28ama1ya/2y,, 


(13.5) 


where  Xg*  is  the  total  thermal  conductivity  of  a  multi-atomic  gas,  is 
the  coefficient  of  thermal  accommodation. 
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The  values  of  the  constant  b  according  to  [229]  are  listed  in 
Table  3. 


TABLE  3-  VALUES  OF  CONSTANT  b  IN  EQUATION  (13-3) 
FROM  VARIOUS  EXPERIMENTAL  DATA. 


Dwyer  [230]  recently  developed  the  theory  of  thermophoresis  in 
thirteen-moment  approximation  which,  according  to  him,  is  valid  for  the 
entire  range  of  Knudsen  numbers.  Such  an  assertion  is  unfounded,  since 
a  thirteen -moment  approximation  was  developed  for  the  slip-mode,  i.e., 
for  Kn  <  1.  The  use  of  this  approximation  for  describing  coagulation 
kinetics  [231]  yielded  results  substantially  different  from  the  exact 
theory,  even  for  Kn  1.  In  the  range  Kn  -*  00  Dwyers'  result  does  not 
convert  to  the  accurate  formula  (13.1),  but  predicts  a  monotonic  dim- 
inuition  of  with  increasing  Kn.  Dwyer's  conclusion  that  should 

have  a  negative  value  for  values  of  the  coefficient  of  mechanical  accom¬ 
modation  0.9  <  ac  <  1.0  contradicts  all  known  experimental  data  [232]. 

Gardner  [233]  examined  the  problem  of  thermophoresis  of  an  evaporat¬ 
ing  particle  located  in  its  own  rarefied  vapor  in  the  absence  of  a  carrier 
gas  and  proved  that  the  thermophoretic  force  should  be  about  the  same  in 
this  case  as  in  the  absence  of  evaporation  of  the  particle. 

Experimental  investigation  of  the  thermophores i~  of  particles  of 
silicone  oil  in  Ar,  N2>  C02  and  Hn,  partially  encompassing  the  range 

Kn  >  1,  was  undertaken  by  Schmitt  [228].  The  tests  were  carried  out 
through  observation  of  the  displacement  of  a  single  droplet  in  a  Milliken 
condenser  with  plates  thermostatized  at  different  temperatures.  The 
particle  radii  varied  from  0.7  to  1.2  u  and  pressure  was  10  to  760  mmHg 
so  that  the  range  of  Knudsen  numbers  between  0.05  and  3.3  was  encompassed. 
For  Kn  35  3.3  F^  is  strictly  proportional  to  a2,  which  qualitatively  proves 

the  theory.  In  the  range  Kn  >  0.2  the  experimental  data  are  approximated 
by  formula  (13.5),  and  for  all  gases  except  H2>  b  *  0.38.  For  Kn  88  3  the 
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ratio  of  theoretical  values  of  to  experimental  values  in  Ar,  N2» 

C02  and  H2  are  1.05,  1.21,  1.22  and  1.38,  respectively,  and  experimental 

accuracy  was  lowest  in  the  latter  case.  An  analogous  method  was  used  by 
Schad  and  Cadli  [234]  for  sodium  chloride  and  mercury  particles  with 
particle  radii  of  0.22-1.15  and  0.096-0.145  ,  respectively,  and  here  the 

range  0.2  <  Kn  <  4.0  was  encompassed.  Measurements  showed  to  be 

independent  of  the  thermal  conductivity  of  the  particles,  i.e.,  the 
measurements  confirmed  the  conclusion  of  the  theory,  although  the  measured 
value  was  closer  to  the  Haywood-Einstein  theory  than  to  the  Deryagin- 
-Waldman  calculations. 

Two  partial  cases  should  be  singled  cut  for  examination  of  diffusio- 
phoresis--l)  equimolecular  counterdiffusion  of  two  components  and  2) 
diffusion  of  one  component  through  a  stagnant  second  component.  The 
theory  of  diffusiophoresis  at  large  Knudsen  numbers  was  developed  by 
Waldman  [235]  and  Deryagin  and  Yalamov  [236].  In  case  1  the  force  acting 
on  a  particle  is: 

F d*  -  (8/3)  a\  ‘(2r.kT )'  -'Du  (?*,)„  | mj?  (1  +  %«/8)  - 

-  mil2  (1  +  (13.6) 

where  DJ2  is  the  coefficient  of  mutual  diffusion,  Xj  is  the  mol  fraction 
of  component  1,  (Ax^)^  is  the  concentration  gradient  of  component  1  at 
infinity  from  the  particle.  For  case  2: 


Fett  —  (8/3)  a'-ng  (2~mg^kTy  7  (1  ■+•  (pi”/®)  DJ2 


(13.7) 


Diffusive  phoresis  in  the  transition  mode  in  a  two-component  mixture 
was  investigated  by  Brock  [237]  using  Bkhatnagar's  model  so  as  to  limit 
the  application  of  the  theory  to  molecules  with  similar  dimensions  and 
mass.  For  case  1  and  >  x^,  assuming  that  *  <t>2  =  1:  • 


=  F, 


dk 


'(f’mjjf.'df,  —  0,311  (2m4,m„n  (mgl +mi2)]1  2  ] 
'  2 - 1  ’ 


1  '■> 


Kn 


(13.8) 


where  dj  and  d2  are  the  effective  diameters  of  the  molecules  d^  =  (d^  + 

+  d2)/2. 
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(13.9) 


In  case  2  and  <  x2: 

Fd  =  F<tk  H  • ' '*071  (2mg,  (mK,  -  w^l^'  -’Kn'1) 


The  values  of  F^  in  equations  (13.8)  and  (13.9)  are  taken  from  (13.6)  and 
(13.7),  respectively. 

Waldman  and  Schmitt  [238]  investigated  diffusiophoresis  for  large 
and  intermediate  Knudsen  numbers.  They  measured  droplets  of  silicone 
oil  at  20°  for  various  binary  gaseous  mixtures.  For  case  1  they  found 
agreement  with  formula  (13.6)  in  the  following  mixtures:  C02/CjHg, 

N2/C2H2,  N2/C2H4,  N2/C2H6,  tf2/02,  N2/Ar.  Brock  [237]  compared  the  theo¬ 
retical  values  of  the  ratio  v,/D.~(V  x.)  ,  where  v,  is  the  rate  of  dis¬ 
placement  of  the  particle  in  the  composition  gradient  field,  for  Kn  ■+  ® 
with  the -results  of  Waldman  and  Schmitt  and  found  that  there  was  agreement 
between  nitrogen-02  mixtures  and  the  above  hydrocarbons.  No  agreement 

was  found  for  N2/Ar  and  CC>2/C3Hg  mixtures  between  the  theory  and  experiment. 

this  is  explained  by  the  great  difference  between  the  diameters  of  the 
components,  which  makes  it  impossible  to  use  the  simplified  relaxation 
model  of  the  kinetic  equation  of  a  binary  mixture.  For  case  2  Brock  found 
satisfactory  correspondence  for  asymptotic  values  of  the  ratio  v^/D^ 

(VXj),  obtained  from  theory  and  experiment  for  Kn  -+  »  in  H20/N2  mixtures. 

By  way  of  concluding  our  examination  of  thermal-  and  diffusiophoresis  we 
might  point  out  that  the  agreement  between  the  corresponding  theories  in 
the  transition  range  and  the  experiment  is  satisfactory.  Nevertheless, 
it  would  be  desirable  to  formulate  a  more  general  theory  based  on  as 
strict  approach  to  simulation  of  the  collision  terms  in  the  kinetic  equa¬ 
tions  of  a  binary  mixture. 

Brock  [239]  pointed  out  that  there  should  exist  yet  another  form  of 
motion  of  aerosol  particles,  "photodiffusiophoresis",  caused  by  a  change 
in  the  physical  or  chemical  equilibrium  between  particles  and  the  gas 
atmosphere  under  the  influence  of  electromagnetic  radiation.  For  ordinary 
photophoresis,  caused,  as  we  know,  by  uneven  heating  of  a  particle  by 
radiation.  Brock  proposed  the  term  "photothermophoresis". 
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§14.  Coagulation  of  HDA 


As  Smolukhovskiy  demonstrated,  the  problem  of  brown- in  coagulation 
is  mathematically  equivalent  ot  the  problem  of  condensation  of  a  vapor 
on  a  particle.  Therefore,  everything  that  was  stated  in  §9  concerning 
mass  transfer  is  also  applicable  to  coagulation.  The  only  difference  is 
that  the  Knudsen  number  in  the  case  of  coagulation  should  be  determined 
somewhat  differently--primarily,  the  role  of  the  mean  free  path  of  gas 
molecules  here  will  be  performed  by  the  apparent  (or  effective,  see  p.  5) 
free  path  of  aerosol  particles,  i.e.,  the  average  distance  at  which  cor¬ 
relation  with  the  initial  direction  of  motion  of_the  particlje  is  lost. 

This  distance  is  given  by  the  relationship  =  v^t,  where  Vp  a  /8kT/  nT 

is  the  thermal  velocity  of  the  particle,  i  =  m  ,  t  is  the  relaxation  time, 

B  is  mobility,  m  is  particle  mass.  Further,  in  the  case  of  coagulatior 

i*.  is  necessary  to  replace  l with  the  corresponding  value,  which  relates 

to  the  relative  motion  of  two  particles,  i.e.,  by  Z  /T”.  Finally,  in 

place  of  particle  radio  it  is  necessary  to  substitute  the  radio  of  the 
absorbing  medium  r  =  2a.  In  view  of  the  above  statement,  the  aerosol 

Knudsen  number  is  expressed  as  Kn  ~  l  '/ijlz..  The  values  of  l  and 

P  P  P 

Knp  for  spherical  particles  of  unit  density,  suspended  in  air  under 
normal  conditions,  are  presented  in  Table  4. 


TABLE  4.  VALUES  THAT  CHARACTERIZE  HDA  COAGULATION. 


Particle  Radius,  & 

10 

20 

50 

100 

200 

1 

500 

1000 

l  ,  A 

P 

059 

40* 

300 

.  ! 

!»4  1 

124  | 

113 

Kn  number 

40.4 

10,5 

42,4 

i  .  50 

0.5*  ! 

0.170; 

0,0* 

♦» 

Coagulation  constant, 

- 

. 

K  10^  cm/sec: 

According  to  Smoiukhovskiy 

42  i 

loir 

t,3,8 

54.0 

i*,0 

8,57 

oJ'‘) 

According  to  equation  (14.1) 

4,42 

4.59 

0.20 

0.14 

9.80 

9.07 

14.0 

12,0 

19,8 

11,0 

3!  ,3 
7.1  i 

i 

5,14 

According  to  Sakhni  theory 

At  large  Knudsenp  numbers,  obviously,  the  coagulation  constant  should 
be  expressed  by  the  formula  for  the  number  of  collisions  between  gas 
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molecules,  i.e.,  aerosol  particles  may  be  regarded  as  large  gas  molecules 
the  probability  of  collisions  among  which  is  not  influenced  by  the 
presence  of  another,  light  component  (Lorentz  gas): 


K  **  2  V2*a%  =  4  (3akT/fyi*, 


(14. 1) 


where  p  is  particle  density.  The  validity  of  this  expression  was  proved 
rigorously  by  Hidy  and  Brock  [242],  who  solved  the  kinetic  equation  for 
Lorentz' s  gas.  The  expression  derived  by  them  for  the  coagulation  con¬ 
stant  takes  into  account  the  electrostatic  interaction  among  particles, 
influence  of  external  fields  and  the  presence  of  velocity  gradients,  but 
unfortunately  does  not  include  terms  with  molecular  forces,  which  cannot 
but  play  some  part  in  the  coagulation  of  HDA,  since  their  particle 
dimension  is  comparable  to  the  background  wavelength.  The  intermolecular 
interaction  potential  for  two  spherical  particles  of  identical  radius  is 
described  by  the  Gamaker  equation: 

U  <#>-  5  [l? +  *."«.+  '»(>  +  £)]•  CM. 2) 


2 

where  R  is  the  intercentral  distance,-  A  =  irq  S/6,  q  is  the  number  of 
molecules  per  unit  volume  of  the  substance,  S  is  London's  constant.  From 
the  equation  of  moments  it  is  easy  to  obtain  the  expression 

(14  3) 

where  b  is  sighting  distance  (see  Figure  7),  is  the  minimum  distance 

between  particle  centers,  E  is  the  total  energy  of  the  particle.  Calculations 
by  the  successive  approximation  method  indicates  that  the  function  b(Rm)  in 

the  range  of  values  2a  <  R  <  «  should  have  a  unique  minimum  b  .  When 

m  m 

b  <  b  ,  R  becomes  imaginary,  i.e.,  the  minimum  distance  of  the  trajectory 
mm 

from  point  0  does  not  exist  and  the  particles  should  collide.  When  b  >  b^, 

R^  has  a  real  value  and  this  corresponds  to  the  case  shown  in  Figure  6-- 

absence  of  collision.  For  two  particles  in  motion  E  denotes  kinetic  energy 
of  their  relative  motion  for  V  =  0.  The  average  E  is  therefore  equal  to 

3  kT.  In  the  absence  of  molecular  forces  b  ■  2a,  i.e.,  the  molecular  forces 

m 

are  increased  by  a  factor  of  bm2/4a2  the  collision  cross-section  of  particles 
and,  consequently,  coagulation  constant. 


By  means  of  approximate  calculation  with  the  kinetic  equation  Brock 
[243]  determined  the  coaguatlion  constant  for  intermediate  Kn^  numbers. 


In  his  derivation  each  of  the  particles  was  regarded  as  a  source  0f  molecules 
with  a  perturbed  distribution  function.  By  colliding  with  the  surface  of 
another  particle  such  molecules  create  a  repelling  force  acting  along  the 


axis  of  the  center.  The  formula  derived  by  Brock  is  extremely  complex. 
When  Kn^  <  5  the  correction  factor  for  free-molecular  formula  (14.1), 


which  Brock  derived,  is  the  same  as  one,  but  as  Kn 


P 


is  reduced  to  about 


two,  this  factor  decreases  somewhat  and  becomes  a  fraction  of  a  percent 
less  than  one.  However,  at  **  2  it  begins  to  increase  rapidly  as 


Kn  decreases . 
P 


Thus,  for  Kn 


P 


85  ,1  the  Brock  formula  predicts  coagulation 


constants  greater  than  the  free-molecular  values,  i.e.,  produces  a  result 
that  contradicts  the  physical  concept  upon  which  it  is  based. 


Figure  7.  Trajectory  of  Particle 
in  Central  Field  of  Attraction. 


Obviously,  the  method  described  in  §9  using  formulas  (9.15)  or  (9.16) 
gives  more  accurate  data  concerning  the  coagulation  constant.  The  co¬ 
agulation  constants  calculated  by  this  method  are  listed  in  Tabic  4.  For 
comparison  also  included  are  values  calculated  according  to  equation  (14.1) 
and  also  according  to  Smolukhovskiy ' s  equation,  the  validity  of  which 
has  already  been  checked  by  several  researchers  who  are  working  on  HDA. 


Astakhov 
the  so-called 
derived  gives 


[231],  for  determination  of  the  coagulation  constant,  used 
thirteen-moment  approximation,  but  the  formula  which  he 

correct  values  only  in  the  range  Kn  <  1. 

F 
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The  question  concerning  the  influence  of  molecular  forces  on  coagula- 
tion  at  intermediate  Kn^  numbers  is  extremely  complex,  since  it  is  necessary 

to  take  into  account  the  effect  of  lag  of  molecular  forces  at  distances 
comparable  to  the  background  wavelength.  Up  to  now  no  theoretical  ex¬ 
pression  that  takes  into  account  molecular  forces  has  been  proposed  for 
%  *  !• 


Noteworthy  is  the  conclusion  reached  recently  by  Wang  and  Friedlander 
[244]  to  the  effect  that  during  coagulation  of  aerosols  under  conditions 
where  the  coagulation  constant  is  a  function  of  the  Knudsen  number,  no 
self-perpetuating  distribution  of  particles  with  respect  to  dimensions 
is  established. 

There  has  been  very  little  experimental  research  concerning  HDA 
coagulation.  The  chief  difficulty  in  such  studies  are  the  substantial 
particle  losses  on  the  walls  of  the  apparatus  and  connecting  lines, 
caused  by  the  high  diffusion  coefficient  of  the  particles.  Therefore, 
the  classical  method  of  studying  coagulation--periodic  determination  of 
the  concentration  of  particles  in  samples  taken  from  the  chamber --cannot 
be  applied  to  HDA.  It  is  not  surprising  that  O'Conner  [245]  and  Quon 
[246],  who  used  this  method,  obtained  inconclusive  results  from  the  point 
of  view  of  theory,  specifically,  that  the  coagulation  constant  increases 
monotonically  as  particle  dimension  decreases  in  the  range  of  intermediate 
Kn^  numbers,  and  even  more  rapidly  than  predicted  by  Smolukhovskiy's 

equation.  O'Conner  worked  with  room  condensation  nuclei,  the  composition 
and  shape  of  whose  particles  could  not  be  determined.  Quon  [246]  used 
aerosols  which  are  products  of  the  incomplete  combustion  of  hydrocarbons. 
Even  though,  due  to  the  low  count  concentration,  it  was  possible  to  assume 
that  the  aerosols  consisted  of  individual  particles  rather  than  aggregates, 
the  question  of  their  shape  and  degree  of  polydispersion  remained  un¬ 
answered.  Only  qualitative  information  is  available  in[[101]  concerning 
the  character  of  HDA  coagulation. 

Fuchs  and  Sutugin  [98]  used  the  flow  method  for  investigating  HDA 
coagulation.  An  aerosol  with  a  high,  a  priori  measured  count  concentration 
was  passed  in  laminar  flow  through  a  wide  tube  and  the  ratio  of  the  con¬ 
centrations  of  HDA  at  the  intake  and  outlet  was  measured  nephelometrically. 
Before  nephelometric  analysis  the  HDA  particles  were  enlarged  to  a 
dimension  of  38  0.5  u  by  means  of  condensation  of  dibutylphthalate  vapors 
on  them  in  a  KUST  apparatus  [97],  In  individual  tests  the  diffusion  losses 
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of  the  particles  in  the  tube  were  determined,  for  which  purpose  a  strongly 
diluted  aerosol  was  passed  through  the  tube.  The  initial  count  concentra¬ 
tion  was  determined  ultramicroscopically  after  dilution  of  the  HDA  in  » 
specially  designed  diluter  and  after  enlargement  in  the  KUST  instrument. 
The  coagulation  constants  determined  by  this  method  for0monodispersed 
NaCl  aerosols  with  average  particle  radii  of  25  and  45  A  were  13. 4* 10" 10 

and  15 . 1  *  10“ 1 0  cm3/sec.  The  measurement  error  was  11%.  The  magnitude  of 

correction  b  ^/4a^  to  the  free-molecular  expression  for  the  rate  of 

coagulation,  governed  by  the  influence  of  molecular  forces,  was  also 
calculated.  Using  the  London  constant  calculated  by  Mayer  [247]  from  the 
absorption  bands  of  NaCl  in  the  ultraviolet  range,  a  correction  factor 
of  2.19  was  obtained.  In  consideration  of  this  correction  the  calculated 
coagulation  constants  were  10.2*10'10  cm3/sec  for  a  =  25  X  and  13.4‘10~10 

cm3/sec  for  a  =  4S  X.  The  latter  of  these  values,  at  least,  is  in  satis¬ 
factory  agreement  with  the  experimental  data.  The  discrepancy  between 
experiment  and  theory  can  be  explained  by  some  polydispersion  of  the 
aerosol  and  inaccuracy  in  the  determination  of  the  London  constant. 


The  above-described  procedure  was  used  later  [248]  for  measurement 
of  the  coagulation  constant  of  monodispersed  aerosols  of  dioctylsebacinate 
° 

with  a  =  100  A.  The  number  Kn^  in  these  tests  was  1.54.  The  experimental 

coagulation  constant  was  13. a • io~ 1 0  cm3/sec,  i.e.,  somewhat  higher  than 
given  in  Table  4.  In  the  given  case  the  discrepancy  is  completely  explained 
by  experimental  error  and  the  influence  of  molecular  forces. 


Stockham  [13]  calculated  the  coagulation  growth  of  silver  particles 
and  found  that  his  experimental  data  can  approximate  the  theoretical 
values  if  values  2-10  times  greater  than  according  to  the  free-molecular 
expression  are  used  for  the  coagulation  constant.  The  discrepancy  between 
theory  and  experiment  increased  as  the  concentration  of  the  vapor  from 
which  the  aerosol  was  formed  increased.  This  is  obviously  related  to 
the  formation  of  branched  chain  aggregates  during  coagulation  of  aerosols 
with  a  high  weight  concentration.  At  low  concentrations,  as  follows  from 
electron  microscopic  photographs  made  by  Stockham,  the  aggregates  were 
compact,  and,  under  these  conditions  the  coagulation  constant  was  about 
two  times  greater  than  the  free-molecular  value,  which  is  in  good  agree¬ 
ment  with  the  factor  of  the  influence  of  molecular  forces  or.  silver  particle 
coagulation,  equal  to  2.21,  calculated  in  [«>]. 


From  the  above  results  it  can  be  concluded  that  there  is  satisfactory 
agreement  between  theory  and  experiment  in  connection  with  the  question 
of  HDA  coagulation. 


The  theory  of  coagulation  of  charged  aerosols  in  the  molecular  and 
transition  modes  is  mathematically  equivalent  to  the  theory  of  the 
charging  of  aerosol  particles.  In  the  free-molecular  mode  the  coagulation 
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constant  of  charged  aerosols  can  be  calculated  with  the  aid  of  an  expression 
derived  by  Hidy  and  Brock  [242].  For  intermediate  Kn^  the  corresponding 

calculations  were  made  only  on  the  basis  of  the  boundary  sphere  method 
[249]. 

5)5.  Capillary  Effects  and  Structure  of  HDA  Particles 

In  addition  to  the  features  inherent  to  HDA  in  relation  to  transfer 
phenomena,  optical  and  electrical  properties,  there  are  features  in  their 
thermodynamic  proper  vies,  most  of  important  of  which--high  va^or  pressure 
of  fine  particles--is  expressed  in  the  familiar  Kelvin  equation: 


In  (plpd  «*  2Ma/p^?r« 


(15.1) 


where  pQ  is  equilibrium  vapor  pressure  over  a  plane  surface,  p  is  the  vapoi 

pressure  over  a  droplet  of  radius  a,  o  is  the  surface  tension,  p  is  the 
density  of  the  material.  The  Kelvin  effect  becomes  notable  starting  with 
a  droplet  dimension  of  0.1  to  0.3  u,  and  increases  rapidly  as  a  decreases, 
which  is  one  of  the  reasons  for  instability,  and  even  impossibility  of 
producing  HDA  from  materials  possessing  a  vapor  pressure  that  is  not  too 
low  at  room  temperature.  Lamer  and  coworkers  [95,  250]  subjected  relation 
(15.1)  to  indirect  experimental  checking  on  droplets  of  sulphuric  acid  and 
dibutylphthalate.  Recently,  direct  observations  were  made  of  the  rate  of 
evaporation  of  liquid  particles  of  Pb  and  Bi  and  of  solid  Au  particles 
under  the  electron  microscope  on  a  carbon  substrate  [251].  These  tests, 
in  which  particle  dimensions  were  100  to  500  X,  also  verified  the  validity 
of  the  Kelvin  equation. 

Higher  pressure  over  a  surface  of  positive  curvature  is  sometimes 
referred  to  as  the  capillary  effect  of  the  first  kind.  Dependence  of 
surface  tension  on  surface  curvature  and  dependence  of  the  pressure  in  a 
droplet  and  of  its  density  on  dimension  are  ascribed  to  capillary  effects 
of  the  second  kind. 

Tolman  [252]  was  the  first  to  prove  the  existence  of  the  dependence 
’(a)  is  the  logical  outcome  of  Gibbs'  thermodynamics.  He  derived  the 
integral  equation: 


r  t»\,  u  \  Jj») 


In  u  m 
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where  6q  is  the  magnitude  of  the  interatomic  distance.  In  the  first 
approximation 
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Here  is  surface  tension  of  a  plane  surface.  The  theory  of  surface  tension 

of  curved  surfaces  underwent  further  development  (253,  254,  255],  in  which 
the  inaccuracy  of  Tolman's  conclusions,  related  to  equalization  of  pressure 
in  a  droplet  and  in  a  macroscopic  volume  of  liquid,  is  c'iminated.  In 
these  works  the  pressure  was  calculated  by  means  of  various  equations  of 
state.  For  6Q  the  following  integral  equation  was  derived: 

j  (R  +  30)  IPo  “  Px*(R)\  •»  0  (15.4) 

“00 

where  R  is  the  radial  coordinate;  p^  is  the  tensor  component  of  pressure 

describing  the  forces  acting  perpendicular  to  R  on  an  element  of  spherical 
surface.  The  function  o(a)  is  of  tremendous  importance  in  the  theory  of 
spontaneous  condensation,  although,  as  pointed  out  by  Kirkwood  and  Baff 
[256],  it  is  not  permissible  to  jnsider  this  function  alone,  ignoring 
the  reduction  of  density  of  fine  particles.  Radial  density  distribution 
and  surface  tension  of  particles  with  diameters  3-25  times  greater  than 
the  molecular  dimension  were  calculated  by  Flesner  [257]  using  the  equation 
of  state  for  solid  spheres  possessing  an  attraction  potential  inversely 
proportional  to  the  sixth  power  of  distance.  The  result  reached  in  this 
work  was  astonishing- -the  surface  tension  of  such  droplets  is  independent 
of  dimension  and  amounts  to  ~ ( 1/3)  This  means  that  in  some  range  of 

dimens ions --rather  narrow,  since  experimentally  the  function  o(a)  does  not 
hold  true  all  the  way  up  to  values  of  curvature  of  10'6  cm*1, --there  should 
be  a  sharp  reduction  in  the  value  of  c.  The  latter  seems  improbable. 
According  to  Plesner,  macroscopic  density  is  not  reached  in  extremely 
fine  particles,  even  at  their  center. 

Shcherbakov  [257,  258]  stated  that  for  extremely  fine  particles  the 
excess  free  energy  inherent  to  them  is  relate)  to  the  entire  volume  of 
the  droplet  as  a  whole,  and  not  simply  to  the  surface.  He  also  found 
that  the  simple  statement  of  the  function  a (a)  in  the  Hibbs -Kelvin  equation 
cannot  lead  to  the  correct  .esult  and  that  it  is  necessary  to  add  to  this 
equation  a  tern  that  contains  do/drr.  Hibbs  and  Kelvir  previously  expounded 
an  analogous  point  of  view,  although  their  opinion  was  based  only  on 


intuitive  premises.  Shcherbakov  examines  surface  tension  as  the  amount  of 
excess  free  energy  related  to  unit  area  of  a  particle.  The  function  o(a) 
according  to  Shcherbakov,  has  the  form: 


(15.5) 


where  1/k  *  2--5  A.  Due  to  the  increasing  influence  of  the  third  term  on 
the  reduction  of  a  the  function  o(a)  has  a  minimum  and  at  extremely  small 
a  increases  as  dimension  decreases.  This  result  is  in  disagreement  with 
the  results  published  in  [252-255]  and  with  a  conclusion  recently  reached 
by  Bcllemans  [2b0]  in  his  generalization  of  the  Mayer  theory  of  the  equation 
of  state  for  the  case  of  a  system  containing  an  interphase  boundary  of 
arbitrary  curvature. 

The  Cahn-llilliard  theory  [261,  262]  and  the  Hart  theory  [263]  of 
surface  tension,  based  on  the  consideration  that  the  density  of  free 
energy  is  a  function  of  concentration  and  concentration  gradient  of 
molecules,  also  predict  a  monotonically  increasing  function  a (a).  Accord¬ 
ing  to  this  theory,  at  high  supersaturations  in  dense  vapors,  the  particles 
of  the  liquid  phase  have  an  extremely  diffuse  surface,  and  their  density 
can  differ  only  slightly  from  the  density  of  the  vapor.  Developing  this 
theory,  Hart  [264]  arrived  at  an  extremely  important  conclusion.  The  fact 
is  that  in  Hibbs ’• theory  of  surface  tension  the  value  of  o  is  constant  in 
relation  to  the  choice  of  position  of  the  interphase  boundary  in  the  trans¬ 
ition  region  between  the  phases.  For  a  curved  surface,  however,  this 
invariance  does  not  hold  true  and  the  calculated  value  of  a  depends  on  the 
position  of  this  boundary.  In  this  case  the  position  of  the  boundary  and 
the  value  are  found  from  the  condition  do/dR  =  0,  where  R  is  the  radial 
coordinate.  Hart  established  that  under  certain  conditions  the  function 
c(R)  may  not  have  a  minimum,  i.e.,  the  concept  of  surface  tension  becomes 
meaningless.  Okuyama  and  Zung  [187]  analyzed  theoretically  the  dependence 
of  the  evaporation  coefficient  on  curvature;  these  authors  stated  that 
the  existence  of  this  dependence  is  the  consequence  of  the  need  to  take 
into  consideration  the  effort  expended  in  increasing  the  area  of  the  droplet 
after  the  combining  with  it  is  one  more  molecule.  Considering  this 
circumstance  they  found  that 


h 


(15.6) 


where  x  n  is  the  condensation  coefficient  on  a  plane  surface,  e_  is  the 
cU  0 

energy  activation  of  adherence  to  the  surface  at  0°  K.  In  the  opinion 

of  these  authors  •  ^  is  equal  to  the  surface  energy  of  the  droplet,  divided 

by  the  number  of  molecules  n  within  it.  Then 
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(15.7) 


o  —  ««,  exp  t—  Irtfo/nkT). 
c 

Substantial  reduction  of  "[subscript  illegible].  calculated  usln«  V 

O 

becomes  notable,  starting  with  a  <  100  A,  and  for  water  at  0°C  and 

O  O 

a  ~  1  A,  a£  decreases  by  a  factor  10°  in  comparison  with  a^.  Obviously, 

consideration  of  the  function  o(a)  should  greatly  weaken,  if  not  eliminate 
this  effect. 

Reduction  of  the  melting  point  of  fine  particles  is  described  by  the 
Thompson  formula: 


Vm  *  2oV/ja  (15*8) 

where  ATm  is  reduction  of  the  melting  point  compared  to  the  macroscopic 

value  T  ,  v7  is  the  mol  volume  of  the  liquid,  L  is  the  heat  of  fusion.  The 

capacity  of  HDA  particles  to  remain  for  a  long  period  of  time  in  the 
amorphous  state  and  their  capacfty  to  agglomerate  or  stick  together  on 
coming  into  contact  are  related  to  a  red  ction  in  the  heat  of  fusion.  The 
results  of  an  entire  series  of  studies  [32,  36,  46]  indicate  that  particles 
of  condensed  HDA  are  usually  amorphous.  This  fact  is  in  agreement  with  the 
Ostwald  law  of  stages  and  is  unique,  considering  that  the  free  energy  of 
formation  from  a  vapor  of  the  seed  of  the  liquid  phase  is  less  than  for 
the  formation  of  crystalline  seeds.  The  crystallization  of  aerosol  sediments 
of  aluminum  [32]  and  the  oxides  of  several  metals  [36]  requires  heating  or 
x-ray  irradiation.  It  was  shown  in  [45]  that  NaCl  particles  measuring 

~10  ^  cm  iast  for  a  long  time  in  the  amorphous  state  in  dry  air.  In  humid 
air,  however,  large  particles  recrystallize  after  a  few  hours  and  the 
fine  particles  immediately  .'ecome  crystals. 

The  capacity  of  highly  dispersed  particles  to  adhere  or  agglomerate 
on  contact  has  an  extremely  strong  effect  on  the  character  of  their  co¬ 
agulation.  As  we  know,  the  Tamman  temperature,  at  which  solid  bodies  begin 
to  agglomerate,  is  equal  in  the  absolute  scale  to  0.6,  the  melting  point 
for  many  bodies  [265],  The  existence  of  the  dependence  of  the  Tamman 
temperature  on  particle  dimension  was  proved  theoretically  by  Higuchi  [266]. 
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A  sharp  reduction  in  the  specific  area  of  the  sediment  of  nickel 
aerosol  with  a  particle  dimension  of  ~30  X  at.  room  temperature  after 
several  hours  was  described  in  [38].  It  is  interesting  to  note  that  the 
influence  of  the  composition  of  the  atmosphere,  in  which  the  sediment 
was  located,  on  this  process  could  not  be  determined.  Coagulation  of 
silver  aerosol  with  a  particle  dimension  of  100  X  is  accompanied  by 
aglutination  and  adherence  of  particles  even  at  120-150°  [6],  i.e., 
much  below  the  normal  Tamman  temperature.  Reduction  of  this  temperature 
by  several  hundred  degrees  for  the  sediments  of  the  HDA  of  a  number  of 
refractory  materials  was  noted  in  [37].  It  is  precisely  the  capacity 
of  particles  to  adhere  during  coagulation  that  limits  the  production  of 
highly  dispersed  powders  through  the  aerosol  state  [6]. 

The  authors  express  their  appreciation  to  I.  B.  Stechkina  and 
V.  I.  Smirnov  for  discussion  of  certain  questions. 
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